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Science and Society” 


By DAVID SARNOFF 


President, Radio Corporation of America, New York 


I' is a high honor to be invited to address this 
gathering of men who search for the kind of 
knowledge on which our present industrial civi- 
lization is based. When we analyze the economic 
and social structure of the world in which we 
live, we find that it rests upon a vast foundation 
of scientific discovery. The steam engine, the 
development of electric power and its applica- 
tion, the internal combustion engine, the funda- 
mental steps in chemistry and metallurgy which 
made large-scale steel manufacturing possible, 
were vital technological inven- 
tions. They brought in their train 
such social institutions as the 
factory, the large city, the trade 
union, advertising, instalment 
selling, and mass production. 
Underneath these technolog- 
ical inventions there lies a pain- 
fully accumulated store of fun- 
damental knowledge concerning 
the properties of matter. That 
is the domain of the physicist. 
Without the method of 
scientific thought, and the phy- 
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sical laws discovered and developed by the great 
physicists—Aristotle, Bacon, Newton, Faraday, 
Maxwell, and many others who lived before the 
present century—there would have been none 
of these inventions, or they would have been 
evolved at a much slower rate, and certainly 
not to the same degree of perfection. 

The 


century 


industrial revolution of the nineteenth 


the ‘‘machine age’’—had its inception 
in the work of the physicist. Since the earliest 
days, it was he who led the way in the search 
for undiscovered facts and prin- 
ciples of nature. Upon these, 
piece by piece, we have reared 
our present economic structure. 

Few of us fully realize the 
long, continuous chain of evo- 
lution that has led from the 
purely scientific and theoretical 
work of men like Newton and 
Maxwell to the technology which 
we see all about us. Few men 
who are not professional physi- 
cists realize the implications for 
the future of humanity that lie 
in the theoretical and experi- 
mental work being done by 
the leaders in physical research 





today. The relationship of the physicist to the 
far-reaching social changes which follow his dis- 
coveries is one of great importance to all who 
are concerned with human progress, whether 
they be scientists, business men, or government 
officials. 


Il. 
One could scarcely find a better example of 
the evolution of a scientific idea 
pure mathematical 


starting in 
theory and achieving a 
practical result affecting the daily lives of all 
civilized people—than the one provided by the 
industry with which I happen to be associated, 
that of radio. It is one of the most recent fields 
of applied physics. But it is also one in which 
the debt to pure physics can be most clearly 
traced. 

Radio had its origin in the purely theoretical 
reasoning of James Clerk Maxwell, a professor 
at King’s College, London, when in 1865 he 
advanced reasons for the existence of electro- 
magnetic waves. Twenty-two years later, Pro- 
fessor Heinrich Hertz at Bonn University, in- 
spired by Maxwell's theoretical work, proved by 
experiment the actual existence of these waves 
and their through space. 
Following his publications other university sci- 
entists began experimenting with electric waves. 
Work went on in the laboratories, and crude 
transmitting and receiving arrangements were 
devised. At this point, Marconi, in Italy, be- 
came interested. In 1895 he invented the 
elevated radiator, or antenna as we now call it. 


ability to travel 


By means of improved transmitting and re- 
ceiving equipment, a ground connection, and 
a telegraph key, he put together the first com- 
mercially successful method of transmitting 
electric waves through the air over consider- 
able distances, and gave to the world a practical 
system of wireless telegraphy. From that point, 
development of the new invention progressed 
rapidly. 

The social results of radio have been far- 
reaching, and the end is not yet in sight. Dr. 
William F. Ogburn, professor of sociology at the 
University of Chicago, not long ago compiled a 
list of 150 social effects directly traceable to 
radio. From Maxwell's theory of electromag- 
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netic waves to 150 social results, the links in 
the chain are continuous. 

Since the beginning of the twentieth century, 
many men in your profession have explored the 
physical phenomena which have given us the 
basis of the greatest developments in radio. 
You are all aware of the tremendously important 
part the vacuum tube plays in our industry. 
This device is the direct result of early work 
by physicists in the study of electron emission 
from heated filaments. Their work stimulated 
further study in the application of high vacuum 
technique to the development of modern tubes. 
Many men who became famous in other fields 
of physics contributed substantially to the early 
growth of radio. 

Today the radio industry includes not only 
countless radio devices as such, but sound- 
amplifying systems, sound-motion-picture ma- 
chines, photoelectric apparatus, and all sorts of 
vacuum tube applications. Every branch of 
physics has contributed to the creation of these 
things: mechanics, heat, light, sound, atomic 
physics, electricity and magnetism. Mechanics is 
applied to obtain uniformity of movement in 
sound-motion-picture machines and in phono- 
graphs; heat, in the development of high power 
vacuum tubes and loudspeakers; light, in the 
optical systems of photographic sound recording 
and reproducing equipment; sound, in loud- 
speakers and other vibrating mechanisms, and 
in the design of broadcasting studios and 
auditoriums; atomic physics, in vacuum and 
cathode-ray tubes; and electricity and mag- 
netism, in all circuit design, as well as in con- 
struction of all magnetically operated devices. 

In our utilization of particles of the atomic 
system we stand today where the early astron- 
omers stood in their exploration of the heavens. 
They studied the planets and the stars which 
seemed to be the heavenly bodies nearest them. 
In our knowledge of the atom,. we first dis- 
covered and utilized the negative electron— 
the outermost and most easily accessible struc- 
tural element, and the one which, in a sense, is 
nearest to us. Some day we shall know more 
about and doubtless utilize some of the other 
elementary nuclear particles which have been 
discovered in recent years—protons, neutrons, 
positrons, deutrons, and their various combina- 
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tions. These new discoveries in turn may give us 
new vacuum tubes, new sources of power, new 
modes of travel and communication, new manu- 
facturing processes, new forms of illumination, 
new cures of dreaded diseases, new highways to 
health. 

Even as the astronomers penetrate farther 
and farther into the depths of space, using ever 
more powerful telescopes, so does the physicist, 
with his bombardment apparatus, penetrate 
deeper and deeper into the atom. The fact that 
one of his most powerful tools, the cyclotron, 
utilizes an ultra-short wave radio transmitter as 
one of its basic elements, illustrates the rela- 
tionship which exists among all the physical 
sciences and their applications. 


II. 


Radio—which grew from the seed planted by 
physicists to the point where it affects the life of 
nations—has “‘arrived,”’ but only at an early 
station on its journey. We are just beginning to 
enter, in any practical way, the fascinating 
domain of ultra-high frequencies—in which 
radio sight will be added to radio sound. 

Short wave of pictures and 
printed or written material has been an accom- 
plished fact for several years. It is now in daily 
service between Europe and America. The 
broadcasting of a facsimile newspaper into every 
business office and home—in half-hourly instal- 
ments if desired—is perfectly feasible. The 
establishment of such a service is now an 
economic rather than a technical problem. 

The new art of television is also making 
progress. For the past ten months we have been 
transmitting experimental television programs 
from the Empire State tower in New York 
City to receivers in the hands of engineers at 
observation points throughout the metropolitan 
area. We believe that acceptable standards of 
picture definition, to which transmitters and re- 
ceivers will be synchronized, have been reached. 

No field of applied science leans more heavily 
on all branches of physics than does television. 
This is particularly true with respect to the 
recent work in atomic physics. Unexplained 
electronic phenomena occur in “‘iconoscopes,” 
the devices which convert light into electrical 
currents, and in the cathode-ray 


transmission 


tubes, or 
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‘‘kinescopes,”’ which convert electricity back into 
light. Here is an absorbing and fruitful field of 
research for the modern physicist. 

The major obstacles to the public introduc- 
tion of television are no longer in the field of 
research and engineering. They lie in a new 
domain. Television now demands the creation of 
a new art form, allied with, yet distinctive 
from, the arts of the stage, of the motion 
picture, and of sound broadcasting. It requires 
new talent, new techniques of writing, direction, 
and studio control. It must set in motion an 
ascending spiral whereby good programs create 
a demand for receiving sets, thus creating a 
growing audience, which in turn will make 
possible better programs. Television must build 
networks, and justify an economic base capable 
of supporting an expensive program service. 
These are some of the problems of television, 
solution of which will one day make it a major 
industry. They are the kind of problems that 
sprout from a seed planted in the soil of pure 
physics. 

An evolutionary sequence of thought, experi- 
ment, production, and use, similar to that which 
constitutes the history of radio, might be traced 
in almost every product or service at the com- 
mand of mankind today. The food we eat, the 
clothes we wear, the roofs that shelter us—as 
well as our services of communication, trans- 
port, health, education, and recreation—all bear 
witness to the impact of scientific thought and 
processes upon human activities that are older 
than recorded history. 


IV. 
Yet that impact had scarcely begun to make 
itself felt a century ago. Prior to the 19th 


century, technological advances were largely the 
result of ingenuity or guesswork or trial-and- 
error. The brave soul who dared to dream 
dreams and see visions in the realms of pure 
science was regarded as a sort of witch-doctor— 
possibly harmless, or perhaps a suitable subject 
for an inquisition and burning at the stake. 

I congratulate you, gentlemen, upon having 
been born to your distinguished calling in a 
more enlightened age. Society today honors you 
and acknowledges its debt of gratitude to you. 
Our radio, the newspapers and the magazines 


375 











reflect a new public taste for scientific subjects. 
They dramatize your discoveries. Every man 
who reads or listens can acquire a smattering of 
scientific knowledge, and converse about elec- 
trons and relativity and light-years. A light- 
year is something like a billion dollars: we have 
learned to talk about it even if we can’t quite 
comprehend it. 

The scientific approach to the solution of 
human problems, the cumulative effect of our 
ever-increasing store of knowledge, the speeding- 
up of technological improvements through the 
coordination of discoveries in diverse fields 
these are assets of such recent origin that the 
ink with which mankind has entered them upon 
its balance sheet is scarcely dry. 

V. 

It is one thing, however, to have assets, and 
another to know how to use them. The inductive 
processes of thought by which science lays brick 
on brick—testing each as it goes along to make 
eventually an impregnable structure—are not the 
processes we use in our attempt to construct a 
better social order. Whereas the scientific world 
glorifies truth, the social world still operates 
largely with prejudice. Thus it comes about that 
we have assets which are products of the scientific 
type of mind, but we do not yet know how 
adequately to use those assets to make them pay 
dividends to humanity. 

Humanity’s balance sheet is all right. It 
shows resources that are priceless, and a net 
worth that is beyond computation. But the value 
of assets without yields is purely theoretical. 
The document we need to be concerned about 
is the human profit-and-loss statement. 

The scientist and the sociologist have view- 
points which are perhaps too widely separated. 
The scientist is engaged in the pursuit of truth, 
the knowledge of which is the most valuable of 
our acquired assets. The sociologist is con- 
cerned with the ultimate effects of truth upon 
human behavior. To the scientist the discovery 
of truth is an end in itself. To the sociologist it 
is only a means to an end. The scientist is 
interested in our balance sheet. The sociologist 
is interested in our profit-and-loss statement. 

In an industrial age, the difference between 
profit and loss for humanity is the difference 
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between use and misuse of the products of 
industry. 

In bygone days the principal products of the 
industrial anvil were simple plough-shares, but, 
as science has given us infinitely richer and more 
varied materials with which to work, our modern 
industrial output has become more and more 
complicated. Today the sociologist rightly may 
claim that many of the gifts of science and in- 
dustry are in the nature of a two-edged sword. 
It is a sword which, like the Nothung of Siegfried, 
can be used to slay the dragons of ignorance, 
intolerance, and greed; but there is always a 
chance that it will turn out to be a weapon with 
which civilization may destroy itself. 

One does not have to go far to find illustra- 
tions of the blessings and dangers that go side 
by side in the discoveries of science. The chem- 
ical that safeguards the work of the surgeon can 
poison the city’s water supply. The airplane, 
that speeds transportation and commerce, can 
drop bombs from the air to blow women and 
children to atoms. 

There are other dangers, less obvious and 
more subtle. Radio, for example, can be used 
for propaganda and regimentation, as well as 
for education and entertainment. Science laid 
the same gift of radio at the feet of society in 
Europe as it did in America. It is true that in the 
United States there is room for improvement in 
some of the programs broadcast on the air, and 
that we still have to learn how to derive the 
greatest social benefit from radio. But no one 
raised in the tradition of liberty and democracy 
can doubt that our use of it is in the direction 
of social betterment, and that in certain parts 
of Europe, where radio has been commandeered 
by the forces of regimentation, its misuse points 
toward social degradation. 

The new art of television has similar poten- 
tialities to build up or tear down social values. 
Like sound broadcasting, it can make friendly 
neighbors of people who differ in race, creed, 
politics, and language; while at the same time 
it offers a powerful weapon to the war-maker, 
and a medium of propaganda for the autocrat. 

Fortunately, man’s deliberate abuses of the 
gifts of science are the exception rather than the 
rule. Otherwise we would have to agree with 
the ancient and pessimistic philosopher who 
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asserted that, ‘““He who increases knowledge in- 
creases sorrow.’’ Such a belief would warrant 
converting our universities into trade schools, 
and our doctors of science into expert mechanics. 

On the other hand, it is too optimistic to 
assume that the mere translation of a scientific 
discovery into a usable commodity or instrument 
always advances civilization; that just because 
humanity can travel faster, communicate more 
freely, cook, wash, iron, and gather ice cubes 
with less effort than ever before, it has reached 
the all-time peak of civilization. Giving a man 
a hoe or a microscope does not make him a 
farmer or a scientist, and giving him a radio or 
an automobile does not make him civilized. It 
is the use that society makes of the products of 
science and industry that determines whether 
civilization is advanced or retarded. 

VI. 

Civilization depends for its advance upon our 
expanding knowledge of the social as well as the 
physical sciences, for society, no matter how 
benevolent its intentions, cannot solve its 
problems by intuition or rule-of-thumb. It must 
develop its own standards and technique. Yet it 
is only natural that the growth of social science 
should lag behind that of physical science, for 
the social scientist is dealing with human rather 
than inanimate materials; with moral rather 
than mathematical equations; with experiments 
which yield their final results, not in the course 
of days or years, but only in generations and 
centuries. 

Experimentation in any field is an undertaking 
of the pioneering spirit—never that of the 
reactionary and the standpatter. The advance 
of social science, no less than that of physical 
science, calls for the creative imagination of a 
Newton and a Maxwell, an Edison and a Mar- 
coni. Obsolescence is a factor in social as well 
as in industrial machines. In society as in in- 
dustry, to stand still is to go backward. 

The social experimenter is important to so- 
ciety. But if he is ruled by his emotions, and 
abandons the spirit of research for that of 
crusading, there is danger for society. The spirit 
of the social crusader is not always the spirit 
which seeks for truth. To attempt to cure the 
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ills of society with untested remedies is not the 
method of science. 

The social experimenter must be one who 
approaches the problems of society in the clear 
light of an unbiased mind. He must collect and 
analyze facts, seek to fathom causes and estab- 
lish principles. He must always be willing to 
reconsider those principles in the light of their 
results. 

It is, of course, impossible to draw a com- 
plete analogy between the world of science and 
society itself. The engineer deals with things 
that have no will of their own. His is the will, 
which can command them to take a certain 
shape and serve a certain purpose. He who 
would be an engineer of society must deal with 
human beings, who are flesh and blood, emo- 
tion and prejudice, dreams and aspiration. They 
each have a will, and in a democratic society 
people must be persuaded and not coerced. But 
the scientific spirit can direct the means of their 
persuasion. A surgeon does not try to persuade 
a patient to undergo an operation until he has 
made a careful diagnosis, and can convince the 
patient that the operation will help him. He 
appeals not to his prejudices or to his fears, but 
to his sense. 

If the cause of civilization is to be advanced 
on a permanent basis, we must learn to follow a 
procedure similar to that of the experienced sur- 
geon. We must not be afraid to operate; but it 
is vitally important to know when to operate 
and when to advocate a less spectacular and 
safer treatment. One decision may require as 
much courage as the other. We must investigate 
and verify social and economic facts; analyze 
them in order to arrive at logical conclusions; 
and obtain widespread acceptance of those con- 
clusions through popular education and debate. 

All these are slow, laborious processes, and 
may prove irksome in a fast-moving age. But 
they are necessary if we are to preserve a free 
society. The enlightened citizen accepts the 
fact that the momentum of society, like the 
momentum of physical bodies, is the product of 
mass times velocity. The velocity of the few 
often has to wait upon the inertia of the many, 
and it is only by overcoming that inertia that 
genuine social progress among a free people is 
achieved. 
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VII. 


Industry today is following the vanguard of 
science into new and infinite realms of knowledge. 
It would be a rash astronomer who said that he 
had calculated the outermost limits of space, 
beyond which there is nothing. It would be a 
rash physicist who claimed that he had dis- 
sected the atom into its ultimate, indivisible 
fragments. Science and knowledge have no 
boundaries. 

So it would be a rash economist who pre- 
dicted any limit to the tangible results of scien- 
tific thought in the form of new goods and 
services placed at the disposal of mankind. In 
fact, it is only by a constant development of new 
goods and services that we may expect to re- 
engage the man-power released by technological 
improvements in established industries. The 
market for every new commodity eventually 
reaches a saturation point and becomes primarily 
a replacement market, so that a more efficient 
technology reduces the number of workers 
needed in that field. 

But science is simultaneously creating new 
employment, both by the modernization of 
established industries and by the creation of 
new ones. In our own generation we have seen 
the automobile, the airplane, the motion picture, 
and the radio provide totally fresh fields of 
activity for millions of men and women. Many 
of our older industries have engaged scientists, 
with notable success, to develop new and re- 
model old products to meet the needs of a 
modern era. 

The industry which has not learned how to 
employ scientists to make it new, and keep it 
new, is doomed. Few industries are so stagnant 
as not to be aware of this; but there are some so 
conservative that the scientist is called upon to 
turn salesman and show them how modern 
science can rejuvenate them to meet present-day 
realities and survive. 

In all respects, I hope we can bring about a 
closer understanding and cooperation between 
you, who seek new truths in the universities, and 
leaders of industry, who seek to make the truth 
you discover serve society. Not only research 
staffs but industrial managers should at all times 
be kept informed of your new discoveries. ‘With 
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such knowledge, promptly obtained, I am cer- 
tain we can shorten the time-gap which now 
separates technological unemployment and use- 
ful re-employment. 

Any measure of unemployment relief ob- 
tained by placing a check-rein upon technology, 
or by arbitrarily hampering men’s efficiency, is 
unsound, uneconomic, and cannot endure. 


VIIL. 


The problems created by technical science 
must be solved by increasing and applying our 
knowledge. of social science. Society cannot ex- 
pect to see every detail of its future plotted in 
advance, any more than the scientist in his 
laboratory knows beforehand where each new 
experiment will lead him. But a study of the 
social strains and stresses in the world today 
points toward a clearly defined objective which 
should be the immediate goal of social science. 
That goal should be to achieve economic justice, 
peace and prosperity in a free democracy. 

The autocracies of Europe have come into 
power by exploiting fears of insecurity and by 
promising impossible Utopias. Fear is an enemy 
of progressive development, whether of indi- 
viduals or of nations: fear of economic and 
social forces too powerful to be struggled 
against; fear of insecurity, of injustice, of un- 
predictable economic changes. 

The American pioneer in the pre-machine age 
did not know this kind of fear. His enemies 
were the forces of Nature from which he could 
wrest a living with the strength of his two hands. 

That strength is no longer enough. The indi- 
vidual is no longer self-sufficient. The fears of 
today must be banished by the organized effort 
of an enlightened society, capable of bringing 
into equilibrium the economic forces that huge 
populations and the industrial age have set in 
motion. 

Because a collective effort is called for, how- 
ever, it is a false and dangerous assumption to 
think that the freedom of the individual must be 
turned over to the state. The governments which 
have demanded such a price have not been able 
to deliver the economic security they promised 
in exchange. But even if they could make good 
their promises, the belief that any material gain 
is worth the loss of individual liberty denies the 
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spirit and tradition of America, and ignores the 
bitter lessons of history. 


IX. 


Freedom of the individual is essential to the 
full expression of his creative faculties—whether 
in art, in science, or in technology. Only out of 
free minds has modern civilization been created; 
only out of free minds can it be advanced. 

To reach our goal of economic and _ social 
freedom, it is not sufficient merely to have an 
idealistic viewpoint. Idealism is only wishful 
thinking until it is given shape and direction by 
fact-finding and logic. We must learn to apply 
science to the business of government and to the 
government of business. 

At a time when the people of many countries 
are being encouraged by their leaders to believe 
that salvation lies in autarchy and self-suffi- 
ciency, it is well to remember that the minds of 
men in all nations have collaborated across all 
borders to advance the service of science to 
society. Without such collaboration, progress 


will be slowed up, and may even stop. Regard- 
less of how the human emotions may operate, 
the human mind knows no national boundaries. 

Science and Society depend upon each other 
for their welfare and advance. Science teaches 
a basic lesson—that knowledge of the truth, 
without fear or prejudice, is indispensable to 
progress. Where that lesson is rejected, and 
scientists are ordered to shape their reasoning to 
fit the purposes of an autocratic government, 
science goes backward instead of forward. The 
tragic position of a scientist in a regimented 
society was recently revealed by a professor in 
one of Europe’s oldest universities, when he 
made this public statement: ‘‘We do not know or 
recognize truth for truth’s sake or science for the 
sake of science.”’ 

Under such a doctrine, Science and Society 
can only serve each other badly. Democracy 
believes in a different philosophy, written in 
words of wisdom nearly two thousand years 
ago: ‘Ye shall know the truth, and the truth shall 
make you free.” 





As a result of all my reading and observation it seemed to me that improvement in 
the standard of living of the human race could come about in the future only by 
reason of new discoveries and inventions, just as, in the past, the steam engine and 
other inventions had been responsible for many improvements in the standard of 
living enjoyed by the average man today. It was these things, and not governmental 
or political action, that had increased production, lowered costs, raised wages, 


elevated the standard of living and so had brought about a greater participation of 


the human race in these benefits. 


ANDREW W. MELLON—Address at the Dedication 
of the New Building of Mellon Institute, Pitts- 
burgh, May 6, 1937 
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Alloys 


for 


Tools and Wear Resistance 


By ZAY JEFFRIES 


Incandescent Lamp Department, General Electric Company 


OME twenty years ago there appeared in one 

of the French periodicals a statement to the 
effect that when metals are subjected to great 
hydrostatic pressure the hardness is substantially 
increased. That didn’t seem to me to be a logical 
thing. I didn’t see how a metal could be hard- 
ened without being plastically deformed—that 
is, a metal which was not susceptible to heat 
treatment. Dr. Bridgman collaborated with me 
at the time in some experiments, and we demon- 
strated quite conclusively that when the pressure 
is truly hydrostatic there is no plastic deforma- 
tion and there is no change in hardness that 
could be measured. 

Dr. Bridgman’s diagnosis of the experiments 
made in France was that the fluid which had been 
used was one which itself solidified under high 
pressure, and in that way the supposed hydro- 
static pressure was nonhydrostatic and made 
possible a certain amount of plastic deformation 
which, in turn, gave rise to the hardening. 

In connection with tools, I want to say in 
passing that when one gives a title, “‘Alloys for 
Tools and Wear Resistance,” he is not neces- 
sarily limited to the hard materials. Many 
hammers used industrially are made of lead. 
The old time hammers, of course, were hard, and 
most of the present day ones are relatively hard. 

There are many uses for soft hammers. One 
of them is in the molding art, and especially in 
connection with permanent molds. These molds 


* Presented at the Symposium on Metals, Cambridge, 
Massachusetts, January 28-30, 1937. 
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must be hammered from time to time to separate 
one part from another, and it is essential that 
the mold part remain unimpaired. A lead or 
aluminum hammer may be used for that purpose. 
Of course, the hammers easily deform and quickly 
wear out, but their replacement is not expensive. 
It is well to remember that tool materials are 
not necessarily hard. 

However, it is my intention to discuss in a 
broad way certain of the harder tool materials. 
We have all read of the profound effect that tools 
have had on the progress of civilization. It has 
been stated that the progress of peoples can be 
traced through their development of tools and 
tool materials. That is quite obvious from his- 
torical studies. 

The application of tools to our modern indus- 
trial development, however, is not quite so 
obvious. What is happening in the industrial 
world is the transfer of hand tools to machines. 
Consider, for example, the anvil, the hammer, 
shears, punches, knives, saws, chisels, tongs and 
the like. These came into existence originally as 
hand tools. For centuries, much of the progress 
of the civilized world depended on, the develop- 
ment and use of these as hand tools. 

Our oldest tool material of any great moment 
is carbon steel. The carbon steel tool is still used 
as a typical tool material in tremendous quanti- 
ties. It is used, however, for machine parts in 
even greater quantities. Take, for example, 
lathe chucks. These involve some combination 
of the tongs and the vise, being a mechanical 
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grip to serve the same purpose as the hand 
grip. 

The use of tool steels for such purposes is very 
much greater than it is for tools themselves. 
We can visualize literally thousands of types of 
industrial machines today whose function is 
essentially that of holding tools which simulate 
the old time hand tool and even the old time 
hand operation. In one machine, many combina- 
tions such as the grip, the file or its equivalent, 
the saw, the shear, and the punch, operate 
simultaneously or in quick succession and pro- 
duce parts all the way from hundreds per minute, 
down. 

The evolution, then, of a tool material is from 
its use as a tool, and usually in the old days as a 
hand tool, to a machine part. Let us consider 
one use of what used to be a typical tool steel 
for ball bearings and roller bearings. The steels 
used for these parts are tool steels. The tonnage 
of ball bearings and ‘roller bearings today is 
greater than that of the same materials used for 
the simple tools. Other examples will occur to 
you, illustrating the tremendous use of these 
old time tool materials in machines and as part 
of the equipment of our modern industrial 
civilization. 

When high speed steel came into use in the 
early part of this century, many people thought 
it would replace carbon steel in tools. It did re- 
place many carbon steel tools, but there are more 
carbon steel tools made today than there were in 
the first part,of the century. What high speed 
steel has done is to make possible a greater 
machine age, and it has taken up part of the 
growth of the tool industry, but it has not 
decreased the output of carbon steels. 

High speed steel first came into use as a tool 
material, and now it has many other uses. It is 
destined to have still more uses outside of the 
typical tool field. The same is true of stellite, 
which is essentially an alloy of cobalt, chromium 
and tungsten. When it is cast it has full hardness. 
It is not machineable, like high speed steel, 
because it has no soft state. High speed steel 
has the virtue of being soft enough in one state 
to be machined. After it is machined it can be 
hardened to any desired degree up to about 
65 Rockwell ‘‘C.”’ Stellite must be cast to shape, 
or nearly to shape, and formed finally by grind- 
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ing. It is true that it can be machined to an 
extent with diamond-pointed tools or with 
some of the newer cemented carbide tools, but 
the forming of stellite tools by machining is not 
a large factor at the present time. 

Siellite first came into use as a tool material, 
and now it is being used in greater quantities for 
machine parts, not only for hard surface facing, 
but also as cast stellite parts themselves, for 
wear resistance. 

The newer tool materials are the cemented 
carbides, representing a new art. The utilization 
of powder metallurgy in the production of these 
materials represents an advanced step and a new 
step in the tool art. This industry at the present 
time is well launched the world over. The ma- 
terial handles more like diamond than like a 
metal at the present time. On the average, in the 
United States, a pound of the material is di- 
vided into ninety pieces before it goes into use. 
Even so, in the year 1936 upwards of fifty tons 
of material were consumed throughout the 
world. Even the cemented carbide has already 
started in the direction of use for machine parts 
as distinguished from the tools themselves. I 
have in mind one machine which uses 25,000 
cemented carbide wearing points. Within a 
period of two or three months there was supplied 
the material to equip this particular machine 
and also one die for drawing steel tubing. The 
die in one piece weighed five and a half pounds. 
The 25,000 pieces for the one machine weighed 
five and a half pounds. The application to the 
machine was equivalent to making 25,000 small 
tools, because each piece of cemented carbide 
was brazed or welded onto the end of a small 
steel shank. 

The above is mentioned only as an example. 
The cemented carbides are being used success- 
fully in many wear resistance applications in 
many types of machines, as, for example, in 
lamp-making machinery. One of the parts in one 
of the machines is a clamp which grips a piece 
of tungsten wire and holds it tightly while it 
pushes one end into hot glass: that is for the 
purpose of making the filament supports. Until 
the cemented carbides were available, these 
clamps had to be changed sometimes three times 
a day. Before a change, as a rule, there would be 
some defective mounts. With cemented carbide 
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Turning and facing a cast iron textile machinery part 
with carbide tools at 1170 r.p.m. Production per hour: 150 
pieces. Number of pieces between tool grinds: 12,000 to 
1400. 


clamps the machines can run along for weeks 
without any attention. 

Another example is in the application of 
cemented carbides to cutters and anvils for 
flattening wire. The difference in life is at least 
three or four hundred to one over the best 
materials available prior to the use of cemented 
carbides. 

The advantage, of course, is not only in the 
longer life itself and in the fewer defective parts. 
Some of the advantage derives from lessening 
the attention required for the operation of the 
machine. If one has to change a part two or 
three times daily, the machine needs almost 
constant watching during operation. If the part 
has to be changed only once in six months, the 
operator need give little attention to that par- 
ticular part. 

The cemented carbides are used for valves and 
valve seats. The control valves in certain ice- 
making machines are made of cemented carbide 

not only the seats, but the needle valves as 
well. There are other types of valves which are 
relatively inaccessible, such as those in mine 
pumps and pumps which must handle not only a 
clear liquid but slimes and materials containing 
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grit. It has been found that the lives of valves 
made of cemented carbide under such circum- 
stances may be hundreds of times those made of 
the ordinary valve materials. 

Although the cemented carbide industry is 
young, it has within a few years changed from 
use as a tool material solely to application in 
machine construction. Hence it is undergoing in 
the course of a few years the type of evolution 
through which the ordinary carbon tool steel 
took thousands of years to pass. 

I wish to say something about the intermediate 
material, the alloy of iron, cobalt and chromium. 
By substituting iron for the chromium of stellite 
and changing the cobalt and tungsten content 
materially, an alloy can be made which is not 
hard when it is first cast. In fact, it can be made 
of a hardness which is machineable. This alloy 
can also be forged and rolled, and then it is 
susceptible to heat treatment. There is no carbon 

or at least not enough to produce carbon 
hardening. The carbon content may be up to a 
quarter of a percent without changing the hard- 
ness much and without interfering with the 
heat treatment. After this alloy is forged and 
rolled, it can be heated to a high temperature, 
say 1250° centigrade, and cooled fairly rapidly, 
after which it is relatively soft and machineable. 
The hardness at that stage is about 45 Rockwell 
“C.”’ It can then be reheated to 550 or 650° 
centigrade and the hardness can be brought up 
to any value desired between 45 and about 
70 Rockwell “C.”’ 

The alloy has some very strange characteris- 
tics. In certain conditions it is so brittle that 
when dropped on the floor it will break. Also, 
if the high temperature treatment is not high 
enough, for instance 1150 instead of 1250°, the 
hardness can be developed with ease but the 
material has poor cutting properties. This is an 
unexplained phenomenon for some of the 
younger metallurgists and physicists to work on 
in the future. Why do we get the requisite hard- 
ness and yet poor cutting properties? 

The alloy can be given good cutting qualities 
by suitable heat treatment, and it has some 
rather peculiar properties still not explainable 
when it is produced in the preferred state. It 
has a peculiar behavior in machining the austen- 
itic steels. It machines Hadfield’s manganese 
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steel better than anything that we have ever 
seen, the cemented carbides included. It ma- 
chines the stainless steels better than any other 
tool material. Why that is I do not know, but I 
am offering it as one of the unexplained facts. 
In the course of time we shall know more about it. 

I may say, too, that we do not yet know why 
this material is so weak after a certain heat 
treatment. That is another unsolved problem, 
so far as we are concerned. 

Even this material, new as it is, and used to 
only a very limited extent for tools at the present 
time, is already finding use outside the tool field. 
It is newer material than cemented 
carbide and yet, within two or three years, it is 


an even 


going through the same type of evolution which 
took thousands of years for the ordinary tool 
material. 

Finally, I wish to mention some things about 
the carbide materials which may be of interest. 
Of all the foundation carbides, tungsten carbide 
is, so far, the best. Why that is we do not know. 
It seems to have what Dr. Hoyt stressed a 
number of years ago in connection with tool 
steels—‘‘body.” It is relatively insensitive to 
temperature changes, and has a toughness out 
of proportion to its strength as determined by 
transverse tests. It is a fact of experience that 
tungsten seems to be the best material. Tungsten 
is the foundation of nearly all the commercial 
grades of cemented carbide. The most important 
of these materials is straight tungsten carbide 
and cobalt, the cobalt ranging, according to the 
toughness desired, from about three percent to, 
say, twenty percent. 

There are other materials, however, which 
have found definite places in the art. One is 
the titanium 
carbide. Some material is made containing no 


tantalum carbide, and other is 
tungsten carbide, but it is very small in pound- 
age—that is, tantalum carbide cemented with 
nickel or cobalt, or both. Most of the tantalum 
carbide is used in conjunction with tungsten 
carbide. Cobalt is used as the cementing metal, 
although nickel can be used, and both nickel 
and cobalt in the same mixture can be used. 
The tantalum materials do not crater in steel 
turning or steel cutting as badly as do the 
straight tungsten carbides. That phenomenon 
is not understood. If twenty percent, say, of 
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Carbide shearing blade for shearing steel tag wire. Carbide 
tools shear 1,000,000 pieces between grinds as compared 
with 40,000 to 60,000 with high speed steel tools. 


tantalum carbide is added to tungsten carbide, 
the steel cutting qualities of the material are 
greatly enhanced. 

Titanium carbide acts in somewhat the same 
way, but different in magnitude. The materials 
which can be used at the highest temperatures 
are those which contain large amounts of ti- 
tanium carbide. The cutting speeds attainable 
in steel are bevond the capacity of most existing 
lathes and machine tools at the present time. 
Soft steel can be cut up to fourteen hundred 
surface feet a minute. Very few commercial 
machines will stand that speed. Why the ti- 
tanium behaves as it does is just as much of a 
mystery as that relating to tantalum carbide. 

The surface of this field has barely been 
scratched. While these three materials are all 
used industrially and are important at the 
present time, that doesn’t mean that the end 
with these 


reached in connection 


cemented tool materials. 


has been 

There is one recent industrial development 
which I think is very important. Even though 
the tantalum and titanium carbide 
materials will cut steel, in this country the art of 
cutting steel with the carbides has advanced 
slowly. In Europe about thirty-five percent of 
the cemented carbide is used for steel cutting. 
This is easy to explain. They have a large number 
of single tool Conditions 
fixed for favorable operation with the carbides. 
Here, to a large extent, a machine must perform 


carbide 


machines. can be 
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many operations, some at low and some at high 
speeds. High speed steel is well adapted to 
handle such operations because it cuts well at 
low speeds. The carbides do not cut steel so 
well at low speeds. It has been a problem to 
out material and 


tool design to cut at low speeds as well as at 


work combinations of tool 
high speeds. 

Within the last year some of these complicated 
machine set-ups have been operated entirely 
with cemented carbide tools. By selecting tool 
materials and angles—in other words, by de- 
signing tool materials for the work—and by 
using them at the proper speeds and feeds, 
steel jobs are now being run entirely with ce- 
mented carbides. During the past year much 
effort has been expended on this problem and 
much information has been forthcoming. Inas- 
much as the steel cutting field is much larger 
than the cast iron or nonferrous field, I think 
this development has great potentiality. 





From these remarks you might say that the 
cemented carbides are on the way to replace 
high speed steels. That is no more true than it 
was true in 1900 that high speed steel was 
destined to replace carbon steel. The cemented 
carbides will take up a proportion of the growth 
which we hope is on the way in this country. 
High speed steel, fifteen or twenty years from 
now, will be used in greater tonnage than it is 
today. In other words, it will be a bigger industry 
as the years go by, just as the carbon steel tool 
industry, or its equivalent, is a greater industry 
today than it ever was before. 

I look forward 


to the growth—the rapid 


growth—of the cemented carbide industry and 
the growth of these newer alloys which, inci- 
dentally, will provide wonderful materials for 
structural application in connection with high 
temperature equipment, but I anticipate also an 
increased use of ordinary tool steel, stellite and 
high speed steel. 





T here need be no fear of science becoming the prostitute of commercialism, as some 
seem to think. On the contrary, the widespread and growing invasion of industry by 
scientists is a wholesome sign. The scientists of our industrial research laboratories 
are the same men who, twenty years or so ago, found opportunities for work only in 
our universities and schools, or perhaps in the medical profession, or in some 
relatively obscure scientific bureau of the government. Indeed, many of them have 
come to industry from these sources, where they were honored and revered for their 
records of able and unselfish service to mankind. 

Industry is changing these men but little, if at all. The mere act of moving from 
one city to another, from one laboratory to another, from an old allegiance to a new 
one does not change human character. These are the men who have given ideals and 
altruism to our uniwersities, and industry is now giving them opportunity for 


expression in a new field. It is asking but one thing of its scientists 


the practical benefit of their fellow men. 
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URING the last five years, the electron 
theory of metals has come from a point 
where the mere explaining of electrical conduc- 
tivity was something of a feat, to a stage where 
its interpretation of alloys should be a definite 
help to practical metallurgists. Further progress 
in the theory is likely to come more and more 
from cooperation between theoretical physicists 
and metallurgists, and the advance of physical 
metallurgy in the next few years is likely to be 
almost as dependent on the electron theory of 
metals as the advances of the last generation 
were dependent on thermodynamics and the 
theory of solutions. For the first time we can 
make a clear picture of what metallic solid solu- 
tions and intermetallic compounds really are, 
and as the theory progresses their properties are 
bound to be more and more clearly understood. 
The ordinary concepts of chemical valence, of 
course, are not applicable to the metallic state. 
Alloys are compounds of variable composition, 
and the occurrence of intermetallic compounds of 
practically fixed composition is determined more 
by peculiarities of the crystal structure than by 
valence laws. The forces between atoms in a 
metal, then, are not of either of the two familiar 
chemical types. They not electrostatic 
forces between ions of opposite sign, as in the 


are 


ionic crystals of inorganic chemistry; neither are 
they homopolar valence forces, as found in the 
familiar inorganic gaseous molecules, and in 
organic compounds. The first thing which the 
modern electron theory has done is to provide 
us with a definite picture of the nature of the 
forces which actually operate within a metal. 
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In the last analysis, of course, the forces are 
electrical ; the quantum theory has not changed 
the fundamental fact that all interatomic and 
intermolecular the 
in some 


forces of significance in 
structure of matter are electrical, or 
minor cases magnetic. An atom of an electro- 
positive element, or metal, contains several 
loosely bound electrons, surrounding a much 
more stable positive core, consisting of the 
positive nucleus and a number of negative elec- 
trons. The loosely bound electrons are the ones 
which can be removed when the atom in question 
enters an ionic crystal. But in a metallic lattice 
these electrons are not removed, for there is no 
electronegative element to attract them. They 
merely become free to wander through the crys- 
tal, thereby carrying electric current. The picture 
of a metal, then, is roughly a sea of free electrons, 
in which are floating enough positively charged 
ions to make the whole thing electrically neutral. 
Different metals vary greatly in the fraction of 
the volume taken up by the positive ions. In an 
alkali metal, like sodium, the fraction may be 
only ten percent, the other ninety percent being 
taken up by quite free electrons, while*in the 
transition groups of elements the may 
occupy fifty percent or more of the space. Now 


ions 


electrostatics alone is able to explain the co- 
hesion of such a structure. In the metal, the 
volume of the whole structure is less than the 
sum of the volumes of the isolated atoms; that 
means that the free electrons on the whole are 
crowded closer to the positive ions than in the 
separated atoms, resulting in a decrease of 
energy as the crystal is formed, or an increase of 
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Fic. 1. Energy distribution of electrons in crystals. 
Abscissa, energy; ordinate, number of electrons per unit 
energy range. Shaded regions occupied by electrons. (a) 
Free electrons in metal. (b) Separated bands in an insu- 
lator. (c) Partly separated bands in a poor conductor. 


energy as it is broken up, which would show 
itself as the heat of sublimation of the metal. 
This explanation shows that the energy of the 
metal does not depend in any striking way on 
small motions of the ions. In an organic molecule, 
for instance, or in an inorganic structure like a 
silicate, the valence forces act in definite direc- 
tions, like the tetrahedral bonds of carbon and 
silicon, and any displacement which changes 
interatomic bonds, either in direction or in 
interatomic distance, meets with large resistance. 
In a metal, on the other hand, an individual ion 
does not greatly care just how many neighbors 
it has, or in what directions they are located, 
with the result that it is easily deformed, and that 
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the same metal often exists in several crystal 
structures, with different numbers of neighbors, 
or different coordination numbers, in the various 
phases. Furthermore, the metals with small 
ions, like the alkalis, are softer and more easily 
deformed than those with large ions, like the 
ferrous metals. The size of the ions is not the only 
factor affecting ductility, however. It is a 
familiar fact that intermetallic compounds of 
definite composition, and some few metals, 
show a surprising brittleness as compared with 
most metals, and at the same time show un- 
usually small electrical conductivity. Thus in the 
copper-zinc alloys, the gamma-phase, a peculiar 
crystal structure occurring only for a narrow 
range of compositions, is much more brittle 
than the other brasses, and a much poorer con- 
ductor. One of the most conspicuously brittle 
metals is bismuth, which at the same time is a 
poor conductor. Many other examples could be 
given. A more detailed study of the free electrons 
gives the interpretation of these less ductile 
and metallic structures. 

To understand these points, we must consider 
the kinetic energy of the free electrons. One 
might suppose at first sight that at low tempera- 
tures the free electrons of a metal would have 
no kinetic energy, and that at higher tempera- 
tures the mean kinetic energy would be 3/2 kT, 
as with a perfect gas. The situation is quite 
different, however, on account of the principle 
called Pauli’s exclusion principle, whose counter- 
part in statistical mechanics is the Fermi sta- 
tistics. According to this principle, no two 
electrons in the metal can have the same kinetic 
energy, but instead the electrons must be spaced 
out in energy according to a definite law. The 
result is that while the slowest free electron has 
no kinetic energy, it is the only one of this sort. 
There is then one with a small energy, one with 
a larger value, and so on, until the fastest of the 
free electrons has quite a large kinetic energy, 
of a magnitude equal to what it would get by 
falling through a difference of potential of several 
volts. This holds even at the absolute zero; 
increase of temperature, in fact, makes very little 
change in the energy distribution, so that the 
free electrons contribute very little to the specific 
heat. The law of spacing of the energies is 
simple for free electrons: the number of electrons 
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per unit energy range is proportional to the 
square root of the kinetic energy. This is shown 
in Fig. la, where the abscissa is the total energy 
of the electron, the ordinate the number of 
electrons, and where the energies actually oc- 
cupied by electrons are shaded. We also show the 
zero of energy, drawn to be correct for sodium. 
All the electrons, as we can see, have negative 
energy. The total energy is the potential energy, 
which is negative, on account of the attractive 
forces between the positive ions and electrons 
which we have already mentioned, plus the 
kinetic energy, and to get the total energy of the 
metal we should have to take the average energy, 
or the center of gravity of the blackened part of 
the figure. This corresponds to the heat of sub- 
limation of sodium, and the theoretical calcula- 
tions check experiment fairly accurately. 

Now the alkali metals are almost the only 
ones where the electrons act as if they were 
really free. The more difficult an electron is to 
detach from an atom, the less it acts like a free 
electron in the metal. An electron fairly tightly 
bound to an atom can only occasionally jump 
over to a neighboring atom, with the result that 
its average kinetic energy is very small. As a 
result of this, the energy distribution curves 
generally are higher, meaning that more elec- 
trons can be crowded into a given energy range, 
so that the mean kinetic energy is lower. There 
is also another change: the distribution curves 
are no longer smooth and regular, but show a 
structure such as is shown in Figs. 1b and Ic. 
These show two cases: first, in 1b, an extreme one 
in which the energies are separated into distinct 
bands; second, in 1c, an intermediate case where 
there is almost complete separation. In the first 
case, if one band is completely filled, the other 
completely empty, the crystal is an insulator. 
The reason is that an electric field in a conductor 
accelerates the electrons, giving them larger 
kinetic energy, in the process of building up an 
electric current. But in a completely filled band, 
there are no states of higher energy to which 
these electrons could go by the action of an 
electric field, so that no current can be set up. 
Such a pattern, in other words, is characteristic 
of an insulator; that is, in general, of a chemical 
compound held together by ionic or homopolar 
valences. The intermediate case, where the bands 
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are almost completely separated, corresponds to 
a poor conductor, which at the same time ap- 
proaches an intermetallic compound in proper- 
ties. It is this intermediate case which we have 
mentioned above, as resulting in structures of 
low conductivity and low ductility, sort of a 
compromise between a metal and a chemical 
compound, of which two examples are the 
gamma-phase of brass, and bismuth, both 
previously mentioned. 

The particular points which I wish to bring 
out in the present paper, however, do not relate 
to problems of poor conductors. I wish instead 
to see how this theory of energy bands applies 
to the iron group of elements, particularly to 
the ferromagnetic elements, iron, cobalt, and 
nickel, and to their alloys. Fig. 2 shows the en- 
ergy distribution curve, which applies schemati- 
cally for this whole series of elements. The typical 
feature is that there are two distinct bands of 
electrons, overlying each other in energy. The 
low, broad band, much like the distribution 
curve for free electrons, though on this scale it 
looks much lower, comes from the few easily 
removed electrons of the atom, the 4s electrons 
in the language of atomic structure, of which 
most of the isolated atoms of the iron group 
possess either one or two. The sharp, high band 
comes from the next inner shell of electrons, the 
3d shell, which begins to be filled up with scan- 
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Fic. 2. Energy distribution of electrons in the iron 
group. Numbers 1-12 represent the energies to which the 
bands are filled in elements with 1—12 outer electrons. 


dium, which normally has one such electron, 
and which is completed with nickel, which 
normally has ten electrons, the number which 
this band can accommodate. In any one of these 
elements, the band is filled up to a definite 
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point, such that the combined area of both 3d 
and 4s bands up to this energy is enough to 
the total 
In general, this results in a fractional number of 


accommodate number of electrons. 
electrons in each band. Thus in nickel there are 
ten outer electrons, of which about 9.4 electrons 
per atom on the average are in the 3d band, 0.6 
in the 4s; in cobalt, the next preceding element, 
there are 9, about 8.3 in the 3d, 0.7 in the 4s; in 
iron, there are 8, about 7.3 in the 3d, 0.7 in 4s; 
etc. These figures are not very accurate, but 
are approximately correct. Of course, they do not 
mean that atoms really have fractional numbers 
of electrons ; they represent averages over all the 
atoms of the metal, each atom having a whole 
number of electrons. Now the 4s electrons are the 
ones which behave in a typically metallic way, 
forming a sea of charge in which the positive ions 
float. The 3d electrons, on the other hand, are 
attached rather tightly to the ions, as the high 
narrow form of the distribution curve shows, so 
that they are not very active in conduction, 
One might think that they would give rise to 
binding forces like chemical valence, but they 
are rather too far apart to have a very large 
effect of this sort, though they do have some 
effect on binding. As the 3d shell fills up, the first 
electrons entering it of course have the lowest 
possible energy, and later ones have higher 
energies. At first, then, the 3d electrons tend to 
lower the average energy of the electrons, result- 
ing in a progressively tighter and tighter binding 
as more and more are added until finally they 
begin having such high energy that they tend to 
increase rather than decrease the average 
energy, and the binding gets weak again. This 
fits in with the fact that criteria such as heats of 
sublimation, compressibility, melting point, and 
other properties show the binding of the elements 
in the middle of the transition group to be greater 
than the binding at either end of the group. 
The electrical conductivity of the elements 
results largely from the 4s electrons, the ones 
which act mostly like free electrons. Since the 
number of these is less on the average than unity, 
the conductivity is poorer than for alkali metals ; 
other features in addition decrease the conduc- 
tivity, which there is not time to discuss. And 
the 4s electrons are also largely responsible for 
the binding and for the lattice spacing, which is 
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remarkably constant through this group of ele- 
ments, though we have just seen that the 3d 
electrons have some effect here. The most 
characteristic property depending on the 3d 
electrons is ferromagnetism. Consider iron, for 
instance, which has about 7.3 electrons per atom 
in the 3d band, as against a possible 10. To dis- 
cuss magnetism, it is essential to know that the 
individual electrons are magnets, and that it is 
their magnetism which results in ferromagnetism. 
Now in the presence of a magnetic field there is a 
principle which limits these electron magnets 
to two orientations only: their magnetic mo- 
ments can point parallel to, or opposite to, the 
field, but not at any other angle. Furthermore, 
the band of 3d electrons, holding 10 electrons as 
a maximum, is subdivided into two identical 
bands, one for electrons of each orientation, each 
holding 5 electrons. Let us now consider two ex- 
treme cases of magnetism of our iron crystal. 
First, it could be completely unmagnetized, even 
in its smallest elements. Then there will be 3.65 
electrons per atom of each orientation in the 3d 
shell. On the other hand, consider the completely 
magnetized state. Then the best we can do is to 
have 5 electrons per atom parallel to the field, 
for that fills up the band, leaving 2.3 electrons 
pointing opposite to the field. The net magnetic 
moment is then 5— 2.3=2.7 electronic moments 
per atom. A similar calculation gives 1.7 elec- 
trons per atom for cobalt, 0.6 for nickel, and 
shows that copper, the next element beyond 
nickel, for which the 3d shell is completely filled, 





could not possibly be ferromagnetic. Now when 
we convert into practical units, these net mag- 
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Fic. 3. Saturation magnetization, in electronic units, as 
function of number of vacant spaces in shell of 3d electrons, 
for binary alloys Ni-Cu, Ni-Zn, Ni-Co, Fe-Ni, Fe-Co, 
Fe-V, Fe-Cr. 
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netic moments actually agree with the observed 
values of the saturation magnetization for these 
three atoms. If we look back at the figures, we 
see that the general rule is that the saturation 
magnetic moment equals the number of empty 
spaces in the 3d shell, times the moment for one 
electron per atom. Of course, for the elements 
before iron, which are not ferromagnetic, this 
rule cannot be expected to hold. 

Fig. 
saturation magnetization holds for the whole 


3 shows how well this rule regarding 


series of ferromagnetic binary alloys of iron, 
cobalt, and nickel, and some neighboring ele- 
ments. Here we plot the magnetic moment (in 
units of the electronic moment) as a function of 
the number of empty spaces in the 3d shell. The 
significant feature of this plot is the way we have 
included alloys as well as the pure metals. 
Consider for instance an alloy of iron and cobalt. 
Each of these elements would have about 0.7 
electron per atom in the 4s state, so that we can 
assume that the alloy has the same. Then the 
average number of 3d electrons per atom will be 
ta and 8.3 
(as in cobalt), the exact amount depending on 
composition. It is this number which we have 


(as in iron) 


somewhere between 


used as the abscissa in the plot. The graph shows 
that the magnetic moment is almost precisely 
what we should get for a pure metal (if it ex- 
isted) with the number of 3d electrons per atom 
which the alloy has. In other words, the only 
effect of going from a pure metal to an alloy, as 
far as this property is concerned, is to change the 
average number of electrons per atom. We see 
that this result holds, not merely for the iron- 
cobalt alloys, but for Fe-Ni, Ni-Co, Ni-Cu, 
Ni-Zn, Fe-\V, Fe-Cr, all of which lie practically 
on a single line. This line corresponds to our rule 
that the number of empty states in the 3d shell 
gives the magnetic saturation, up to a point 
between iron and cobalt; above this value, the 
saturation moment falls off, as we should expect. 

The method which we have used here, of 
plotting all the alloys as a function of the num- 
ber of empty spaces in the 3d shell, might be a 
useful concept in other applications of the theory 
to alloys. Of course, many properties, particu- 
larly those depending on the entropy change in 
mixing unlike atoms, depend on other things in 
addition to the number of electrons per atom, 
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but still it should be a useful method for cor- 
relating experimental facts. For instance, we can 
consider the phase diagrams of the binary alloys 
of Fe, Co, Ni and neighboring elements from this 
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Fic. 4. Phase equilibrium diagram, as function of num- 
ber of vacant spaces in shell of 3d electrons, for binary 
alloys Cu-Ni, Ni-Co, Ni-Fe, Co-Fe, Fe-Cr. 


point of view. Fig. 4 shows the more important 
for Cu-Ni, Ni-Co, Ni-Fe, 
Co-Fe and Fe-Cr, plotted against this same 
parameter, the number of empty spaces in the 


equilibrium lines 


3d shell. Certain general features become clearer 
in this way than by considering the separate 
phase diagrams. Thus the liquidus is a fairly 
smooth curve, rising toward the middle of the 
transition group at Cr (according to our previous 
statement that the most tightly bound elements 
tend to come toward the middle of such a group). 
In a broad way, the alloys with less than about 
two empty spaces in the 3d shell are face cen- 
tered, those with more are body centered. The 
magnetic transformation, or Curie point (shown 
by a dotted line), starts at —273°C for zero 
empty spaces in the 3d shell, the point at which 
ferromagnetism starts, then rises to a maximum 
at something like two empty spaces; falling 
again in the region of the body centered lattice 
in much the same way that the magnetic satura- 
tion value falls. But the details are very compli- 
cated. In particular, the properties of cobalt are 
quite different from those of an iron-nickel alloy 
with the same number of electrons per atom. The 
transition lines are all much higher for cobalt 
than for the Fe-Ni alloy, and the hexagonal 
phase of cobalt appears, analogous perhaps to a 
hypothetical hexagonal phase of Fe-Ni which 
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might appear at temperatures too low to reach. 
The reason for this is not perfectly clear, but 
presumably is connected with the entropy 
difference between the pure cobalt, where there 
is no entropy of mixing, and the Fe-Ni alloys, 
which has such an entropy effect. This method 
of considering the whole series of alloys as far 
as possible as a one-parameter family has not 
been used much up to the present, but may be 
valuable in visualizing regularities in the series 
which might otherwise escape attention. 

Of course, the methods just described can be 
used only for substitution solid solutions, in 
which the substituted atoms are almost of the 
same size. This is more nearly true in the iron 
group, perhaps, than anywhere else in the 
periodic table. But we must not lose sight of the 
fact that a great deal of the variety in the struc- 
ture of alloys arises on account of the combina- 
tion of atoms of quite different sizes. The theory 
of energy bands suggests, as we should expect, 
that structures consisting of two or more types of 
atoms of different sizes would tend to form com- 
pounds rather than purely metallic solid solu- 
tions. The reason is essentially that valence 
electrons in one type of atom do not have ordi- 
narily the right energy to act at the same time as 
valence electrons in the other type. In a pure 
metal, or an alloy of very similar atoms, of 
course this difficulty does not appear. But if 
the electrons in one atom are not free to wander 
at will to their neighbors, the average kinetic 
energy is decreased, there is less freedom of 
movement and poorer conductivity, and the 
energy bands tend to narrow down, approaching 
the situation which we have seen is to be identi- 





fied with the compound. At the same time the 
interatomic forces become more like valence 
forces, more dependent on the exact arrangement 
and spacing of the atoms, and the structure tends 
to become hard and relatively brittle, rather than 
soft and ductile like the more typical metals. 
Such behavior, of course, is typical of metals 
of the iron group alloyed with quite different 
elements, as carbon, silicon, and nitrogen. The 
theory of such alloys has not been touched as yet 
but we may look forward to important results 
when it can be worked on. 

From what we have seen, the theory of metals 
and alloys is still in an early stage, where it can 
suggest rather than work out the answers to 
problems. Its promise, nevertheless, seems to be 
great. Thermodynamics has been the great guide 
of metallurgy in the last generation. But thermo- 
dynamics has great weaknesses. As far as it goes 
it is correct, but it can give us very little insight 
into the mechanism underlying the properties 
of a physical or chemical system. To get further, 
we must have a model of the system, whose 
properties can be worked out from first princi- 
ples, and from which, by use of the methods of 
statistical mechanics, we can derive values of the 
quantities which thermodynamics must postu- 
late as empirical facts. For the first time, the 
theory of metals is able to see its way clear to 
provide such models, accurate enough so that 
we can really hope to get quantitative informa- 
tion from them without a prohibitive amount of 
work. It does not seem fanciful to hope that this 
theory will bring about as great advances in 
theoretical metallurgy in the next generation as 
thermodynamics has provided in the last one. 





Research is the reconnaissance party of industry, roving the unknown territories 
ahead independently, yet not without purpose, seeing for the first time things that 
all the following world will see a few years hence.—S. M. KINTNER 
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New York, New York 


ET me begin by quoting one or two sentences 

from Professor Van Wert’s new and ad- 
mirable Introduction to Physical Metallurgy. 
Professor Van Wert says that ‘‘Progress in this 
search”’ [the search ‘‘to discover the essential 
nature and meaning of metal behavior” | ‘‘has 
been beset by many obstacles. Experimental 
inquiry has often been hampered by the lack of 
adequate research tools, for in metal research as 
elsewhere investigative enterprises have, on 
occasion, had to wait on the invention or de- 
velopment of more precise analytic or measuring 
appliances or methods.” 

I quote this statement of Professor Van Wert’s 
because it is of a new research tool that I am to 
speak, a tool the lack of which has hampered 
certain types of experimental inquiry in metal- 
lurgy up to the present time. I speak of electron 
crystal analysis as a new research tool because, 
although 10 years of age, it seems not to be well 
known to metallurgists. It is true that metal- 
lurgists here and there have made use of it. 
But it is true also that no mention is to be found 
of either electron diffraction or of electron 
crystal analysis in any book on metallurgy; at 
least there is no mention of either of these 
subjects in any one of four books on metallurgy 
reviewed in Nature only last month. It is fair to 
say, I think, that electron crystal analysis is not 
as yet a recognized tool of metallurgical research. 

This simplifies matters for me, for it allows me 
to treat the subject in as elementary a fashion 
as I please. What I have thought of doing is 
first to make some general statements about 


* Presented at the Symposium on Metals, Cambridge, 
Massachusetts, January 28-30, 1927. 
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What Electrons Can Tell Us 
About Metals’ 


By C. J. DAVISSON 
Bell Teiephone Laboratories 






und 2 


electron waves and electron diffraction, then to 
describe a few typical investigations in which 
electron diffraction has been employed, and 
finally to say something about the technique of 
this new type of crystal analysis. 

The most general statement one can make in 
regard to electron diffraction is perhaps that it is 
very similar to x-ray diffraction. A beam of 
electrons in which the electrons all travel with 
the same speed has a definite wave-length just 
as has a beam of characteristic x-rays. The 
wave-length of a beam of electrons is inversely 
proportional to the momentum of the individual 
electrons which compose it. The present practice 
in electron crystal analysis is to use electrons 
which have been accelerated through from 30 to 
70 kilovolts. For these the wave-length range is 
from 0.07 to 0.05 angstrom units. The wave- 
length of the Ka, radiation of tungsten is greater 
than 0.2 angstroms and this is the shortest wave 
regularly used in x-ray crystal analysis, so that 
one difference between electron crystal analysis 
and x-ray crystal analysis is that in the former 
we work with shorter waves. 

But the really great difference between x-rays 
and electrons is in the extent to which they are 
scattered and absorbed by matter. The char- 
acteristic for which x-rays are. most celebrated 
is their power to penetrate thick slabs of material 
entirely opaque to ordinary light. Associated 
with this transparency which all materials ex- 
hibit to a greater or less extent, is a correspond- 
ingly low scattering power. The specimens 
examined by x-ray diffraction must not be too 
thin or the radiation scattered into the pattern 
will be negligible compared with that remaining 
in the primary beam. 
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a, as prepared. 
b, heated 10 hours at 690°C, 
c, plus 10 hours at 680°C. 


Fic. 1. Electron diffraction patterns of oxide film from aluminum—by 


With electrons the case is entirely different. 
All substances are more opaque to electrons of 
the wave-lengths used for crystal analysis than 
are metals to visible light. Electrons are no more 
suitable for examining sheets of metal by trans- 
mission than metal sheets are suitable for 
replacing glass in windows. To be suitable for 
examination by electrons by transmission, a 
specimen must be no more than a few hundred 
angstroms in thickness It must be just the sort 
of specimen which cannot be examined by x-rays. 

Massive specimens can be examined by elec- 
trons by reflection. The beam is directed onto 
the surface at near grazing incidence, and the 
half-pattern the 


crystalline state of a surface layer of excessive 


which is produéed reveals 
thinness or, more often, the crystalline state of 
tiny points and edges which protrude from the 
surface and transmit parts of the incident beam. 
Invisible films of material, different chemically 
from the bulk of the specimen, are frequently 
discovered by this method. 

The reflection method has, of course, long been 
used in examining metals by x-rays. Here too 
the patterns are determined by material lying 
near the surface. But the layer contributing to 
the pattern is never so slight for x-rays as it is 
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d, plus 50 hours at 680°C. 
e, plus 2 hours at 800°C. 
f, plus 1 hour at 1000°C. 


Preston and Bircumshaw. 


for electrons. An invisible surface film would 
never be found by x-rays; the pattern it produced 
would be swamped by the much stronger pattern 
due to the underlying metal. Thus x-rays and 
electrons are complementary in their uses, rather 
than coincident. Gossamer films is the province 
of electrons—relatively thick and bulky speci- 
mens is the province of x-rays. 

I will now describe very briefly three out of 
some 700 investigations which have been made 
by electron diffraction during the last ten years. 
Of the three I have chosen, one only has to do 
with ferrous metals, but all will illustrate the 
kind of information a ferrous metallurgist might 
expect to obtain concerning his materials. 

The first is an investigation made by Preston 
the National 
Laboratory. Preston and Bircumshaw have been 


and Bircumshaw! in Physical 
using electron diffraction for some time to study 


the oxidation of metals. In their most recent 
paper they describe their investigation of the 
oxide film which forms on aluminum at room 
first the 


method using as specimens small plates cut from 


temperature. They tried reflection 
cold-rolled sheets. The patterns they obtained 
! Phil. Mag. 22, 654-665 (1936). 
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consisted of portions of two broad fuzzy rings 
such as are produced when x-rays are diffracted 
by a liquid. This might mean that the oxide 
formed on aluminum at room temperature is 
amorphous, or it might mean that the surface 
metal had been rendered amorphous by the 
rolling and that the pattern was due to the metal 
rather than to the oxide. To dispose of this 
uncertainty Preston and Bircumshaw dissolved 
away thin sheets of aluminum by gaseous HCl, 
leaving only the oxide films. These they examined 
by transmission and again they found the broad 
fuzzy rings. The oxide film formed on cold-rolled 
aluminum at room temperature is undoubtedly 
amorphous. 

On heating, the film remains amorphous to 
about 700°C at which temperature a pattern of 
sharp rings begins forming. When fully formed 
the pattern is that of y—Al,O . It is identical 
with the pattern produced by oxide films fished 
from aluminum. Patterns obtained at 
various stages of the transition from amorphous 


molten 


to crystalline oxide are reproduced in Fig. 1. 
There are certain peculiarities about the early 
stages of the crystallization, as revealed by the 
patterns, which are not yet understood. 


face after abrasion—more specifically, after the 
face had been smoothed by filing. The results 
which Dr. Germer obtained are perhaps of some 
interest to metallurgists because of a similarity 
between what happens when a galena crystal is 
filed, and what happens when a metal is hardened 
by cold working. 

To understand the patterns produced by the 
abraded crystal one must understand first the 
pattern produced by the normal or undisturbed 
crystal. Galena forms crystals of the rocksalt 
type—a checkerboard of lead and sulphur atoms 
extended into three dimensions. The crystal 
cleaves on its cube face. The reflection pattern 
produced by this face when well etched is under 
most circumstances an array of spots. If the 
primary beam, always at grazing incidence, lies 
near [100] direction in the surface, the array is 
square, Fig. 2. A score or more of diffraction 
beams issue from the crystal with the crystal 
and incident beam in fixed relation to one 
another, which may seem rather surprising to 
those familiar only with x-ray diffraction. 

There are two circumstances which contribute 
to this production of multiple beams; the first is 
that the wave-length of the electrons is quite 





Fic. 2. Cleavage face thoroughly 


etched—beam parallel to [100] di- 


rection. 


Fics. 2, 3, and 4. Electron diffraction patterns produced by a galena crystal—L. H 


Next I will describe an investigation made by 
Dr. L. H. Germer? of the Bell Telephone Labo- 
ratories. The material investigated was galena 
in the form of single crystals. One of the matters 
studied was the crystalline state of the cleavage 


2 Phys. Rev. 50, 659-671 (1936). 
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Fic. 3. After abrasion (filing at 
right angles to primary beam). 


Fic. 4. Same after etching (pat- 
tern due to undisturbed crystal 
showing through that due to rolled 
crystallites). 7 


. Germer. 


small compared to the constant of the crystal, 
and the second is that there is an appreciable 
misalignment among the crystallites which form 
the galena crystal. The electron wave-length is 
less than 1/100 of the linear dimension of the 
crystal cell. The Bragg angles are correspondingly 
small. The smallest of them are considerably 
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smaller than the largest angles by which the 
crystallites fail of perfect alignment. In these 
circumstances Bragg reflections of differing in- 
dices can occur from different crystallites for one 
and the same setting of beam and crystal. This 
is an explanation, though not the only possible 
explanation, of the spot pattern observed in this 
case. For present purposes we need know only 
that this square array of spots is the pattern 
produced by the cleavage face of the undisturbed 
crystal when the incident beam is parallel to a 
cube edge direction, or nearly so. 

face was then abraded—filed 
parallel to a [100] direction. In this condition 
it produced a pattern of rings, Fig. 3. The 
surface material had been broken up and skewed 


The crystal 


° 
a + 10° 





DIRECTION 
OF FILING 


BROKEN UP 
MATERIAL 


__ 
——" RING PATTERN 


/ \ 


UNDISTURBED CRYSTAL ROLLED CRYSTALLITES 
SPOT PATTERN ARC PATTERN 


Fic. 5. Schematic cross section of a galena crystal illus- 
trating the deformation produced by abrasion. 


What 


observed was the powder pattern of galena. 


about into random orientations. was 

The surface of the crystal was then etched 
away a little at a time, and the pattern recorded 
at each stage of the etching. What was observed 
depended upon the relation between the direction 
of the primary beam and the direction of filing. 
When the direction of the primary beam was at 
right angles to the direction of filing the ring 
pattern weakened, some of the rings vanished 
completely and others were reduced to short arcs. 
This pattern of arcs, which exhibited no recog- 
nizable symmetry, persisted without alteration 
through some stages of etching. Then the square 
array of spots began showing through, indicating 
that undisturbed crystal was beginning to out- 
crop, Fig. 4. When this occurred it was dis- 
covered that each of the arcs was associated 
with one of the spots, and that the arc pattern 
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could be described as a rotation and streaking 
out of the spot pattern. It was as if the crystal- 
lites in a layer immediately below the surface 
had been rolled by the file through angles 
ranging from 10 to 30 degrees, and this un- 
doubtedly is what had happened. 

The behavior is illustrated in Fig. 5. The 
drawing represents a cross section of the crystal 
after filing, the direction of filing being down- 
ward. The surface layer is badly broken up and 
disorganized and produces the ring pattern. 
Below this is a relatively thick layer in which 
the crystallites, represented as blocks, have been 
rolled by the file. They roll from 10 to 30 degrees 
and then jam. The material in this layer produces 
the arc pattern. Below this is the undisturbed 
material which produces the spot pattern. This 
exemplifies rather well one of the actions which 
is known from x-ray analysis to occur in work 
hardening. 

The last experiment I will describe is by Dr. 
H. R. Nelson* of the Battelle Institute. Dr. 
Nelson has vaporized iron onto flat surfaces in 
vacuum and has studied the patterns produced 
by these deposits before and after exposure to 
air. One of these is reproduced in Fig. 6. Most 
of the spots or short arcs of this pattern fall on 
equally spaced lines parallel to the surface of 
the specimen. These alone are due to the iron 
and what they reveal is that the iron forms 
crystals as it is deposited, and that a [111] 
direction of each of the crystals stands normal 
to the general surface of the specimen, or nearly 
so. The orientation is that preferred by iron 
when deposited electrolytically. 

The spots and short ares which do not fall on 
the equally spaced lines are due to oxide formed 
on the iron crystals during an accidental 5 
minute exposure of the deposit to air at room 
temperature and at a pressure of about 10-* mm 
of mercury. The fact that the oxide forms a 
symmetrical pattern of spots rather than a 
pattern of rings means that like the iron crystals 
the oxide crystals have a strongly preferred 
orientation. 

With further exposure to dry air, Fig. 7, the 
oxide pattern became stronger, the iron pattern 


' weaker and eventually, Fig. 8, only the oxide 


3]. Chem. Phys. 5, 252 (1937). 
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Fic. 6. Oriented film slightly 
exposed to air. 


Fic. 7. After 2 minutes exposure at 
room temperature. 200°C, 


Fic. 8. After 15 minutes exposure at 


Fics. 6, 7, and 8. Electron diffraction patterns produced by a vaporized deposit of iron during oxidation—H. R. Nelson. 


pattern remained. Dr. Nelson has identified the 
structure of the oxide as that of magnetite and 
the preferred orientations as that in which a 
[ v2, 1, 0] direction in the crystal stands normal 
to the surface. 

It seems odd, unlikely even, that there should 
be anything unique about the [ v2, 1, 0 ] direction, 
for it is not, so to speak, a crystallographic 
direction—these have in all cases rational indices. 
The explanation is, however, quite simple. 
Mehl, McCandless and Rhines* have shown that 
when iron is oxidized at high temperature, the 
orientation of the oxide crystals is related in a 
simple way to the orientation of the iron crystals 
upon which they are formed. We stand the 
oxide crystal cube squarely upon the iron cube, 
rotate the oxide cube 45 degrees about a vertical 
axis, and we have the relation between the two 
structures. One feature of this relationship is 
that every [111] direction in the iron structure 
is parallel to a ['v2, 1, 0] direction in the oxide 
structure. It follows that if the oxide crystals 
formed at room temperature build onto the 
oriented iron crystals in the same way as those 
formed at high temperature, the oxide crystals 
will indeed stand with [v2, 1, 0] directions 
pointing directly outward from the surface. 

These experiments illustrate some of the 
things electrons can tell us about metals. They 
can give us information about the chemical and 
crystalline constitutions of metallic surfaces and 
of thin films of metals and their compounds. 
These, I believe, are not matters to which 
metallurgists are entirely indifferent. 


4‘ Nature 134, 1009 (1934). 
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Let me say a few words now about the means 
employed in obtaining electron diffraction pat- 
terns such as we have seen. The material require- 
ments are a source of electrons, a supply of 
constant high potential, a set of apertured 
electrodes for use in accelerating the electrons 
and forming them into a narrow beam, an 
adjustable mounting for holding the specimen, 
a fluorescent screen for use in adjusting beam 
and specimen, and a photographic plate for 
recording the pattern. These parts except the 
voltage supply are sealed into a _ gas-tight 
chamber, which is exhausted and held by con- 
tinuous pumping at a pressure of 10-* mm of 
mercury or less. 

These features are found in all electron 
diffraction equipments, but the form they take 
varies considerably from one equipment to 
another. Almost every experimenter works with 
a camera of his own design which he regards as 
superior in some respects to all others. 

A diagram of the diffraction camera, used by 
Finch, Quarrell and Wilman® at the Imperial 
College of Science and Technology, is shown in 
Fig. 9. The source of the electrons is a low 
pressure gas discharge which is maintained in 
the chamber at the top. Cathode rays from this 
discharge enter the main chamber through a 
narrow tube and form the electron beam. 
Halfway down the tube is the specimen holder, 
and at the bottom the fluorescent screen and 
photograph plate, which are manipulated by 
mechanical controls passing through the base. 
An important feature of this equipment is a coil 


5 Trans. Faraday Soc. 31, 1051-1080 (1935). 
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magnet which surrounds the tube between source 
and specimen and is used to focus the beam on 
the photographic plate. It is used also, by 
tilting it, to align the beam. The position and 
orientation of the specimen are adjustable from 
outside through gas-tight bearings and a metallic 
bellows. The parts shown near the top of the 
main tube are used for deflecting the beam on 
and off the photographic plate, and for measuring 
the beam current. The beam current used in 
this and other cameras is usually a small fraction 
of a microampere. The exposure times range 
from a tenth of a second or less to several 
seconds, depending on the beam power and the 
nature of the specimen. 

The main chamber is pumped continuously 
and necessarily rapidly since gas enters it 
continuously from the discharge chamber. This 
camera or one of the same type is now manu- 
factured for sale by the Cambridge Instrument 
Company. It is capable of giving very fine results. 

Another embodiment of the same essential 
features is illustrated in Fig. 10. This is a 
diagram of a camera developed at the Bell 
Telephone Laboratories. The electrons come 
from a tungsten filament at A, and are acceler- 
ated from A to B. The beam is defined by small 
apertures in parts B and B’. The alignment of 
the beam is achieved by adjustment of the 
positions of the apertures through sylphon 
bellows and by the use of crossed electrostatic 
fields in the region between B and B’. The 
specimen, in a mounting adjustable through 
bellows, is placed just beyond the second 
aperture. The fluorescent screen closes the end 
of the main tube. The photographic plate is 
held until wanted in the upper part of the 
vertical cylinder shown in the drawing. It can 
be lowered into position for exposure or with- 
drawn by gas-tight mechanical couplings to the 
outside. Adjustable stops are set to limit the 
amount of plate exposed. With exposures limited 
to strips one or two centimeters wide, portions 
of 10 or 15 patterns can be recorded on a single 
plate. This camera is in daily use on corrosion 
and other problems arising in the telephone 
plant. 

‘hardly need point out that what with having 
Fic. 9, Gas discharge type of electron diffraction camera to open and reseal the egeree ” change er 
Finch and Quarrell. mens and plates, and what with having to worry 
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Fic. 10. Filament type of electron diffraction camera—Bell Telephone Laboratories. 


about pressures and pumps, the recording of 
electron diffraction patterns is more bothersome 
than the recording of x-ray patterns. And yet 
when the equipment is properly installed and 
in continuous service the technique is much 
less difficult than might be imagined. Dr. Nelson 
reports that with his apparatus he goes from 
vacuum to vacuum in from 4 to 7 minutes. 

It is rather pointless, however, to make 
comparison of the two techniques, for there are 
few investigations in which one may use x-rays 


or electrons as he chooses. Either the investiga- 
tion is one to which x-rays alone are applicable, 
or it is one to which only electrons are applicable. 
If the latter, one has no choice but to accept the 
experimental difficulties, such as they are, in 
the usual dauntless spirit. 

Metallurgists have found x-ray diffraction 
applicable to so many of their problems it will 
be surprising, indeed, if, in the near future, 
they do not find electron diffraction applicable 
to many others. 





The application and development of scientific ideas in industry depend upon a 
complete understanding of the way in which industry can make use of science and 
scientific method, and this can only be achieved if the problem of cooperation 1s 
studied by the man of science as well as by the industrialist. The differences of r 
experience, training and outlook between men of scientific education and many 
others engaged in productive industry are not always considered by those who are 
anxious to secure the potential benefits which science, wisely used, can bring to 

the community. The scientific worker has to meet the industrialist half way. 


VOLUME 8, JUNE, 1937 


Science, Industry and Society 


(Editorial in Nature, March 27, 1937) 
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i Kaen chemical reactions which occur in the 

corrosion of metals in the ordinary atmos- 
phere are often complex and varied. They depend 
upon the composition, physical state and surface 
condition of the metallic material, as well as 
upon the chemical components of the environ- 
ment,—their phases and concentrations. These 
reactions are affected by temperature and 
temperature fluctuation, by the movement or 
circulation of the air (or electrolyte in case of 
submerged corrosion), by the nature and the 
solubility of the corrosion products and_ par- 
ticularly by the position in which solid corrosion 
products are precipitated with reference to the 
surface of the corroding metal. Moreover, in 
actual cases of corrosion, marked changes may 
occur in the nature and concentrations of the 
surrounding nonmetallic reactants during the 
process of corrosion. It is apparent that a 
rather detailed knowledge of both the metal and 
of the environment in which it is to be used are 
essential in the chemical investigation of corro- 
sion problems. 

The engineering need of information of the 
relative durability of metallic materials has led 
to the widespread use of corrosion tests. These 
vary widely in character from exposure for 
prolonged periods to the actual environments in 
which the metals are to be used to laboratory 
tests of great variety. Such tests, while of con- 
siderable assistance in the selection of suitable 
materials for given purposes, have not generally 
furnished sufficient information as to. the 


* Presented at the Symposium on Metals, Cambridge, 
Massachusetts, January 28-30, 1937. 
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processes of corrosion to permit prediction as to 
what will happen under slightly different condi- 
tions—the problem which usually presents 
itself. 

The relatively long period of time required to 
obtain measurable results on metal specimens 
exposed to natural environments led to an 
interest in accelerated corrosion tests in which 
the rates of reaction are accelerated by increasing 
the temperature, the humidity, or by otherwise 
modifying the natural environment. It soon 
became apparent, however, that the results of 
accelerated tests could not be correlated usually 
with actual rates of corrosion under the condi- 
tions in which the metals are used, and more- 
over, that the quality ratings established by 
such tests were often erroneous. Except for 
certain limited uses, the idea of accelerated 
corrosion testing has been largely abandoned. 

The more recent trend in corrosion investiga- 
tions has been in the direction of the study of 


‘corrosion processes with experimental techniques 


of sufficient sensitivity to measure the actual 
rates of reaction occurring in natural environ- 
ments. In this way significant information can be 
obtained in a reasonable time without the danger 
of .distortion introduced by accelerating the 
corrosion process. The present paper describes 
the application of some of these methods and 
what they reveal concerning the mechanisms 
of corrosion processes. 

In recent years considerable progress has been 
made in the elucidation of corrosion reactions. 
It has been well established, for example, that 
corrosion in the presence of moisture is an 
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electrolytic process in which the metal dissolves 
at certain areas or points—the anodes of small 
corrosion cells, the cathodes of which are the 
adjacent areas on the metal surface at which 
hydrogen is deposited. The driving force of these 
cells arises either from some chemical or physical 
inhomogeneity of the metal or some inhomo- 
geneity of the environment. Their electrolytic 
operation is influenced by the composition, size 
and distribution of the anodic and cathodic areas, 
by the character of the corrosion products and 
by the chemical nature and the conductance of 
the surrounding environment. This statement 
of the replacement of hydrogen from the con- 
tiguous electrolyte by the dissolution of the 
corroding metal does not present a detailed view 
of the actual mechanism which must involve 
the action of the solvent as well as the exchange 
of electrons at electrode areas. The essential 
fact is that metal ions replace hydrogen ions in 
the electrolytic environment. A second step of 
equal importance in the process lies in what 
happens to the electrodeposited hydrogen since 
the continuation of the attack depends upon its 
disposal. It may be evolved as molecular hydro- 
gen by metallic cathodes of low overvoltage, or 
it may be removed by combination with oxygen 
or other oxidants. In either case, when the rate 
of corrosion depends upon the rate of hydrogen 
disposal, the process is said to be cathodically 
controlled. Similarly, when the anodic reaction 
controls the rate of corrosion, as, for example, in 
the corrosion of magnesium in water where the 
rate of corrosion depends upon the rate of dis- 
charge of hydroxyl ions,' the process is said to 
be anodically controlled. In other words, the 
progress of corrosion may be controlled by the 
extent of polarization at either cathode or anode 
or at both. Either process may reduce the effec- 
tive potentials of the corrosion cells virtually to 
zero. 

It is of interest to consider the manner in 
which various environmental constituents in- 
fluence the current voltage or polarization char- 
‘acteristics of corrosion cells. Reference has al- 
ready been made to the depolarizing influence of 
oxygen at cathode areas. In this function it is the 
most important promoter of corrosion. While 
some hydrogen is usually evolved in the corrosion 

1 Whitby, Trans. Faraday Soc. 29, 1318 (1933). 
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of most metals (owing often to the presence of 
impurities of low overvoltage) the extent of 
attack depends mainly in most cases upon oxy- 
gen availability. Nitrates and other oxidants 
play a role similar to that of oxygen. The rate of 
cathode depolarization and therefore the rate of 
corrosion in these cases depends upon the rate 
at which oxygen or oxidants reach the cathode 
areas. This is usually a function of the concentra- 
tion of these substances and of temperature and 
is largely influenced by convection currents. 
There are certain substances, however, which if 
present in the environment may hinder cathodic 
depolarization. Gelatine, diethylamine, hydro- 
quinone, and various colloidal or reducing sub- 
stances inhibit acids 
because of their interference with the cathodic 


the corrosion of iron in 


reaction. Colloidal soils, experience shows, may 
retard corrosive action by assisting in the main- 
tenance of hydrogen films at cathodic areas,’ 
or by otherwise reducing the accessibility of 
oxygen to these regions. 

Anodic processes in corrosion cells consist in 
either the dissolution of metal ions or the pre- 
cipitation of anions at the metal surface. An 
increase in metal the 
accumulation of products is 
reflected in anodic polarization. The solubility of 
compounds of the metal and anions of the 
solution largely determines the course of the 
reaction. Compounds of low solubility tend to 


ion concentration or 


solid corrosion 


form films over anodic areas. The hydroxides 
(or hydrated oxides) of most metals fall in this 
class as do salts such as the silicate or sulfate of 
lead, and phosphate of iron. When, however, the 
corrosion products are more soluble, other con- 
stituents of the environment may, by means of 
secondary reactions, precipitate a solid phase at 
an appreciable distance from the metal surface. 
Thus, in the corrosion of iron, the moderate 
solubility of ferrous hydroxide permits migration 
of ferrous ions from the surface to an adjacent 
zone where, owing to the action of dissolved 
oxygen, the less soluble ferric hydroxide may be 
precipitated as a rust scale. Similarly, lead 
carbonate may be precipitated in the body of the 
solution when lead corrodes in distilled water 
exposed to the air. In neither of these cases does 


the precipitation of the solid phase interfere 


? Burns and Salley, Ind. Eng. Chem. 22, 293 (1930). 
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Fic. 1. Potentials of iron and stainless steel in water and 
certain solutions. 


with the anodic reaction—i.e., induce anodic 
polarization—or interfere with corrosion 

As indicated above, anodic polarization is 
most readily observed in solutions containing 
strong oxidants or anions which form relatively 
insoluble salts with the metal. It is conceivable 
that products comprising these surface films 
have a lesser solubility than the analogous com- 
pounds in bulk.* The process of anodic polariza- 
tion of a metallic surface is a progressive action 
proceeding from the more anodic to the less 
anodic regions until the surface becomes sub- 
stantially film-coated or passivated. 

From the foregoing discussion it will be ap- 
parent that a knowledge of the electrochemical 
polarization characteristics of corrosion cells 
would indicate the nature of the process and 
make possible the prediction of corrosion be- 
havior. It is usually impossible, however, to 
measure the current density potential relation- 
ships of individual corrosion cells owing to a 
lack of knowledge of the electrode areas involved. 
Actually, these areas are probably of a wide 
range of sizes and change in size during the 
progress of corrosion. Even the nature of the 
cathodes is often unknown. Practically, however, 
it is a simple matter to determine the composite 
of the resulting potentials and their change with 
time. If the observed potential of the metal 
surface as a whole tends to become more electro- 
negative with time, it is indicative of increased 
cathodic polarization and the process is probably 





’ Whitby, Trans. Faraday Soc. 29, 415 (1933). 
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cathodically controlled. If, on the other hand, 
the trend of potential with time is in the electro- 
positive or more noble direction, anodic polariza- 
tion and anodic control of the process is apparent. 
These time-potential measurements may be 
readily made with a recording potentiometer. 

In Fig. 1 there is shown the actual time-po- 
tential behavior of specimens of iron and 18-8 
nickel-chromium steel immersed in water and 
certain solutions. It will be observed that the 
potential of iron in normal sulfuric acid rapidly 
becomes electronegative. The initial trend of a 
similar specimen immersed in distilled water is 
toward the noble side, owing probably to the 
passivity of the air-formed film. Within a few 
minutes, this trend is reversed, the specimen 
becoming more electronegative probably because 
of the destruction of this film. In the presence 
of an oxidizing solution, potassium dichromate, 
it will be seen that the potential of iron (repre- 
sented by curve 4) becomes markedly electro- 
positive. A similar passive behavior is shown 
by stainless steel in tap water (curve 5). The 
iron specimens, the potentials of which are 
represented by curves 1 and 3, displayed active 
corrosion, while iron in potassium dichromate 
solution (curve 4) and stainless steel in water 
did not corrode. In other words, a trend in the 
electronegative direction, which, as before men- 
tioned, is characteristic of cathodic polarization, 
is associated with active corrosion, while the 
electropositive trend of potential is indicative of 
passivation. It is of interest to note in this con- 
nection that certain organic compounds used as 
corrosion inhibitors which, as before mentioned, 
function by retardation of the cathodic reaction 
appear to have little or no influence on the 
characteristic of time-potential curves. The ex- 
periment represented by curve 3, for example, 
differs from that of curve 2 only by the presence 
of 0.2 percent diethylamine, a typical corrosion 
inhibitor. Corrosion has been greatly retarded in 
this case but the iron is still in the active or 
corroding class. 

Another instance of the usefulness of time- 
potential studies in revealing information as to 
corrodibility is illustrated in Fig. 2.4 An alloy of 
copper containing 25 percent iron had in previous 


* Schumacher and Souden, Met. and Alloys 7, 95 (1936). 
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tests in another locality shown satisfactory re- 
sistance to atmospheric corrosion, but upon ex- 
posure to New York City air it was severely 
corroded. The time-potential curves show a 
tendency toward passivity which is accentuated 
by providing greater accessibility to oxygen 
(compare the experiment in which specimen was 
immersed in a body of solution, represented by 
curve 3, with that in which the specimen was 
merely in contact with a film of solution, 
curve 2). Passivity is readily destroyed, as 
shown by the abrupt break in curve 2 upon the 
addition of a drop of dilute sodium chloride. 
From this evidence it seems likely that the less 
favorable behavior of the copper-iron alloy in 
New York may be due to the small chloride 
component of the atmosphere. 

A rather unique method of determining the 
relative rates at which similar alloys tend to cor- 
rode is to compare the rates of passivation in the 
presence of film-forming constituents.’ For ex- 
ample, when specimens of lead and certain lead 
alloys are immersed in solutions of 7 normal 
sulfuric acid exposed to the air the initial po- 
tential of lead breaks abruptly after varying 
periods of time to the potential of the cathodic 
areas. A comparison of the rates of sulfation of 
several grades of lead and certain lead alloys, 
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Fic. 2. Potentials of iron copper and an iron-copper alloy in 
dilute sodium sulfate solution. 


determined in this manner, is given in Fig. 3. 
It should be emphasized that the relative surface 
activities of metals determined in this manner 
cannot be correlated directly to rates of corro- 


> Burns, Bell Sys. Tech. J. 15, 20, 603 (1936). 
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sion in most actual environments because of the 
interference of film-forming substances which 
are usually present and which tend to block the 
primary reaction. 

In the experiments described above it was 
found that a potentiometer of ordinary sensi- 


’ 
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RELATIVE RATES OF SULFATION 




















Fic. 3. Relative rates of sulfation of lead and certain lead 
alloys in 7 normal sulfuric acid. 


tivity was entirely adequate and that the record- 
ing feature was of considerable advantage in 
obtaining a complete record of potential with 
time. When, however, it became desirable to 
determine the corrosion behavior of an iron 
surface beneath a paint coating, a vacuum tube 
electrometer circuit proved essential for the 
measurement of the potential of the metal and 
its trend with time. By this means, it has been 
possible to determine whether the inhibitive 
action of primer pigments is due to chemical 
action or to their influence on the physical 
properties of paint films.’ This work is being 
extended to include a wide variety of pigments. 

Still another instrument which promises to be 
of considerable assistance in the study of passive 
film stability is the oscillograph. Fig. 4 (records 
1 to 4, respectively) shows the behavior of the 
potential of passive iron as additional quantities 
(beyond 5 percent) of concentrated hydrochloric 
acid are added to the nitric acid solution re- 
sponsible for the induction of passivity. It will 
be observed that the magnitude of the potential 
break increased with the concentration of the 
film-destroying chloride, and that the time 
required for repair of the film increased with the 
extent of the film failure. The instrument em- 


6 Burns and Haring, Trans. Electrochem. Soc. 69, 169 
(1936). 
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Fic. 4. Time-potential curve of iron in 16m HNO, upon 
adding 5—6 percent 12m HCl. 


ployed in these measurements has a flat response 
characteristic up to 10,000 cycles, and would 
appear to be suitable for the study and compari- 
son of protective films on various metals. Another 
investigator has recently reported results of an 
oscillographic study of the behavior of iron in 
chromate-sulfate solutions.’ 

These potentiometric studies show whether a 
metal is in a corroding or corrosion-resistant 
state, and this is the most important phase of the 
corrosion problem which presents itself in the 
commercial use of metals. As has been seen, 
such studies, while furnishing some indication of 
rates of corrosion, do not provide quantitative 
information on this point. The usual methods of 
determining corrosion rates—changes in weight 
or physical properties or measurement of the 
consumption of reactants—are usually laborious, 
and sometimes of low accuracy. Recently a new 
device of analytical chemistry, known as the 
“polarograph,” has become available and pre- 
liminary experiments indicate that it should be 
useful in the determination of corrosion rates.* 
Essentially, it consists of a cell comprised of a 
capillary dropping mercury cathode, a mercury 
anode of relatively large area, and the electro- 
lyte to be analyzed. In operation, an electro- 
motive force is applied to the cell at a constant 
rate and the resulting current is automatically 
measured and recorded as a function of the 
polarizing voltage. With increasing voltage the 
current rises very slowly until the cathode po- 


tential attains a value at which a given positive 


’ Karschulin, Zeits. f. Elektrochem. 42, 722 (1936). 
* Hevrovsky, Mikrochem. 12, 25 (1932): Kemula, Zeits. 
f. Elektrochem. 37, 779 (1931). 
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ion is discharged (or reduced) at the mercury 
surface. At this point a rapid rise in current oc- 
curs, the magnitude of which is determined by 
the mobility and concentration of the ion. By 
providing a suitable (and constant) ionic en- 
vironment (with lithium chloride in the experi- 
ments to be described) and employing a con- 
stant mercury dropping rate and capillary tip, 
it becomes possible to make rapid determinations 
of certain metals by comparison with solutions of 
known concentrations under the same conditions. 

In Fig. 5 there is shown by means of polaro- 
graphic measurements a comparison of the rates 
of corrosion in distilled water exposed to the 
air of one specimen each of chemical lead and 
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Fic. 5. Polarographic comparison of rates of solution of 
two grades of lead in distilled water at room temperature. 


corroding lead. It will be seen that initially 
chemical lead corrodes more rapidly than cor- 
roding lead, which finding is in harmony with the 
results on their respective rates of sulfation 
contained in Fig. 3. It is of interest to note that 
after 360 minutes, the amount of corrosion indi- 
cated by the polarograms is nearly the same for 
both lead specimens, indicating that the reaction 
after its initial stages must be controlled by 
external factors rather than by the composition 
of the metal. 

The polarographic method is capable of meas- 
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uring one microgram of lead in a volume of one 
milliliter with a precision exceeding +5 percent. 

The high sensitivity of the method permits the 
measurement of the reaction rate in the initial 
stages of a corrosion process, thus contributing 
to the elucidation of the mechanism. In addition 
to speed and high sensitivity, another advantage 
of the method lies in the fact that two or more 
metal ions may be determined in the same 
measurement provided their deposition poten- 
tials (on mercury) differ by as much as 0.2 volt. 
It should prove to be particularly useful in 
investigations of the corrosion of alloys in which 
more than one component dissolves and where 
the nature of the process is obscure. 

From the foregoing account it is seen that 
significant information concerning the nature of 
corrosion processes in the presence of moisture 
may be obtained by utilizing certain techniques 
of high sensitivity in place of the common and 





more cumbersome corrosion-life tests. Moreover, 
it is no longer necessary to resort to the practice 
of accelerating corrosion reactions in order to 
measure them. Observations of the potentials of 
materials in known environments and changes of 
these potential values with time which may be 
determined by means of a potentiometer or a 
vacuum tube electrometer show whether the 
metal surface is in a corroding or passive state. 
The existece and stability of passive surface 
films may be investigated advantageously by 
oscillographic measurements. Rates of corrosion 
may be measured readily by the use of the 
polarograph, the instrument being particularly 
suitable for determining the initial rates of 
reaction. These and other techniques of equal 
promise, such as electron beam diffraction and 
the newer methods of microchemistry, will be the 
means of important contributions to our knowl- 
edge of the corrosion of metals. 





[t seems to me a hopeless task to provide a complete and finished liberal education 
suitable to this century by four years of college work. The only worthwhile liberal 
education today is one which is a continuing process going on throughout life. 
It is from this premise, I suggest, that this question of ‘‘knowledge of wide surface 
and little depth’’ must be judged. I should reject all informational tests applied to 
recent graduates as indicative of the effectiveness of our general education. Whether 
a liberal education has been a success or failure should be measured by the student's 
breadth of vision fifteen or twenty years after graduation. Has the smattering 
acquired in college worn thinner with each succeeding year? If so, it has been of 


little value. Or has it provided a basis for continued intellectual and spiritual 


growth? In this case it has been the most significant part of the college training. 
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SERIES of articles by the author dealing 

with the static strains and stresses in a 
drawn bow have been published in Ye Sylvan 
Archer.' These articles have shown the effect of 
the shape and form of bending of the bow upon 
these strains and stresses. Some interesting and 
valuable information was obtained from that 
work which has materially changed the design 
of modern bows. 

Considerable time has been spent by the 
author in studying, by experimental methods, 
the behavior of a bow while the arrow is being 
discharged. Some of the results of these investi- 
gations have been published in the magazine 
previously mentioned and in The Journal of the 
Franklin Institute (October 1929). He has also 
had high speed pictures (one to three thousand 
frames per second) taken by Electric Research 
Products which, in addition to putting an end 
to all controversies regarding the ‘‘archers para- 
dos,”’ furnish an excellent means of studying the 
motions of the arrow and bow limbs. 

While these experimental methods are inter- 
esting and valuable, they are necessarily slow 
and expensive. An analytical method was there- 
fore developed for determining the dynamical 
forces, the velocities and accelerations of the 
arrow and bow limbs. This method has been 
useful in improving the design and cast of bows. 
The method is briefly set forth in this paper. 
Since the dynamical treatment depends on the 
static conditions, a few of the static equations 
will first be given. 

Let us assume that the bow is made up of a 
short rigid middle section and two limbs, which, 
for all positions of the draw, bend in arcs of 
circles. This assumption is reasonable because 


' Nov. 1930 to Aug. 1932 (published at Albany, Oregon). 
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practically all bowyers construct their bows so 
that they bend in this form. In the articles 
referred to above, it has been shown that a very 
simple design having this form of bending will 
stress all sections of the bow equally. It is 
believed that for most bows this form of bending 
is desirable. If for some particular bow another 
form of bending is desired, the method used 
here may be applied to it. 

Referring to Fig. 1, let B equal one-half the 
length of the bow, L equal one-half the length 
of the rigid handle section, B};= B—L equal the 
length of the active bending portion of each 
limb, S equal one-half the length of the string, 
IT equal the distance from middle of bow to 
middle of line connecting bow tips, P equal 
distance from arrow nock to middle of line 
connecting bow tips, D=H+P equal length of 
draw, Y equal one-half length of line connecting 
bow tips, A equal angle between the line con- 
necting bow tip and point O and the line repre- 
senting the position of the undeflected limb. 
(The point O is located at a distance of 3B,/4 
from the tip of the bow.) Let E equal angle 
made by string with the line connecting the bow 
tips, f equal the static force at each bow tip in a 
direction tangent to its path, 7 equal static 
tension in the string, F equal static drawing 
force, ¥}= Y—B,/4—L and N equal the distance 
along the path made by the bow tip during 
the draw. 

It may be shown that the path made by the 
bow tip is part of a cardioid. The portion of the 
cardioid traversed by the bow tip is almost a 
perfect arc of a circle whose radius is 3B,/4 and 
whose center is located at a distance of 3B,/4 
from the tip of the undeflected bow. Since the 
equations which are to be obtained are based on 
the assumption that the bow tips travel along 
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arcs of circles having radii of 3B,/4, the accuracy 
of this assumption will now be investigated. 

In Fig. 2, let B, represent the length of the 
bending portion of one of the bow limbs. Since 
the limb bends in the arc of a circle, B, is the 
arc of a circle having a radius of r and the angle @ 
corresponds to this arc. The chord K =2r sin 6/2 
and r=B,/@ where @ is measured in radians. 
Therefore K=(2B,/@) sin 6/2. The point Q is 
selected so that QY=3B,/4 and QX=B,/4. 
X Y=B, and represents the length and position 
of the active bending portion of the limb before 
the draw. The angle QXZ=6/2. Let QZ=R. 
Then: 


R=[(B,?/16+ K?—(KB,/2) cos 6/2 }}. 


Substituting the above value of K and making 
use of well-known trigonometric relations it may 
be shown that: 


R= (3B,/4)(1—1604/9|6+---)}. 


Since the angle @ for a bent bow never exceeds 
unity (i.e. 57.3°), all the terms under the radical 
except the first one are negligible. Therefore 
R=3B,/4. (The maximum value of @ for the 
average bow is equal to about %. For this value 
of 6, the second term under the radical would 
equal 16/32,805. Neglecting this term results in 
an error of less than two parts in 1000.) 

Referring again to Fig. 1, the force f is pro- 
portional to the displacement N and is therefore 
proportional to the angle A which corresponds 
to the are N. 

Let f=CA where C is a constant depending 
on the dimensions and material of the bow. 

A subscript 0 will be used to indicate the 
value of any of the above variables when the 
bow is in its braced but undrawn position. 

The following trignometric and algebraic re- 
lationships may now be obtained: 


Y, = (3B,/4) cos A, Y= /Y,+B,/4+L, 
S= Y)=(3B,/4) cos An+ B,/44+L 

=} length of string, 
H=(3B,/4) sin A, TIT) = (3.B,/4) sin Ao, 
E=cos" Y/S, P=S sin E=(S— Y?)!, 
D=H+P, T=f/sin (A+E)=CA/sin (A+), 
To=f/sin Ap=CAo/sin Ao, 
F=2T sin E=(2AC sin E)/sin (A+ 8), 
N=3B,A /4, No=3B,Ao/4. 
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Dynamical Treatment 


Let x represent the distance traveled by the 
arrow in ¢ seconds. Then: dx/dt= V where V is 
the instantaneous velocity of the arrow. But 
x= D’'—D where D’ is the value of D for the 
fully drawn position. (Prime letters will be used 
to represent the value of the variables for the 
fully drawn position. Prime letters will therefore 
be constants. ) 


dx/dt=d(D’ — D)/dt= —dD/dt= V. 


In like manner: —dN/dt=v equal velocity of 
bow tip. But: dD/dN=(dD/dt)/(dN/dt) = V/v 
=R where R is the ratio of the two velocities. 
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Also, 
dD/dN=(dD/dA)(dA/dN). 


But D=//+P. Therefore, substituting the values 
given for D=//+P and differentiating 
dD/dA =(3B,/4) cos A+(3B,Y/4P) sin A. 
But A =4N/3B, and dA/dN=4/3B,. Therefore: 
dD/dN=V/v=R=cos A+(Y/P) sin A. 


Let / equal mass of arrow, m/2 equal effective 
mass of each limb located at tip and W equal 
potential energy of both limbs. Since N is the 
displacement of each bow tip: W’—W=C(A‘N’ 
—AN) is the loss in potential energy of both 
limbs from position of full draw to any other 
position. 

Neglecting the dissipation in heat and sound, 
the gain in kinetic energy is equal to the loss in 
potential energy. Therefore : 


M V?/2+mv?/2=C(A'N’—AN). 
But V = Rv. Therefore: 
v=[2C(A’N’—AN)/(MR*+m) }}. 
But N=3B,A /4. Therefore: 


v=[3B,C(A” —A?)/2(MR*?+m) }}. 
And 
V = Rv=R[3B,C(A” —A?)/2(MR?+m) }}. 


Let a equal the acceleration of the bow tip 
(Bold face 


type will be used to indicate dynamic values of 


and A the acceleration of the arrow. 
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variables. ) 
a=dv/dt=(dv/dA)(dA /dt) 
=(dv/dA)(dA dN)(dN dt). 


But dA /dN=4/3B, and dN/dt= —v. Therefore: 
a= — (4v/3B,)(dv/dA). Therefore by differentia- 
tion: 


a=2CA/(MR?+m) 
—2MRC(A” —A?)Z/(MR?+m)?, 


where 
Z=sin A — ( Y P) cos A +(3B,S° 4P?*) sin? A. 
Also A=dV/dt=d(Rv) /dt=vdR/di+ Rdv/dt. 


dR/dt=(dR/dA)(dA /dt) 
=(dR/dA)(dA/dN)(dN/dt) = —(4v/3B,)dR/dt. 


Therefore 


A=2CAR/(MR?*+m) 
+2mC(A” —A?)Z/(MR?+m)?. 


The dynamical force f on the tip is: f=ma/2. 
The dynamical force F on the arrow is: 


F= MA. 





The author with a shooting machine which he de- 


signed for use in testing arrows. This machine was used 


in conjunction with a chronograph for measuring the 


velocity and acceleration of arrows. It is equipped with 


pneumatic release. 
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The dynamical tension in the string is: 


T=F/2 sin E=SF/2P=SMA/2P 
=(f—f)/sin (A + E). 


If the angle A is expressed in radians, the 
masses in pounds, the lengths in feet and the 
forces in poundals, the velocity and acceleration 
will be obtained in feet per second and feet per 
second per second, respectively. 

It is more convenient, however, to use degrees 
for the angle A, inches for the lengths, pounds 
for the forces and grains for the masses. It is 
also desirable to have the constant C of such a 
value that f=CA where f is expressed in pounds 
and A is expressed in degrees. In order to use 
the units as just explained and still obtain the 
velocities and accelerations in feet per second 
and feet per second per second, respectively, 
the above formulas become : 


v=[492B,C(A” — A?) /(MR?+m) }}, Vince 


a= 450,800CA /(.R?+m) 
—7,868MRC(A” —A?*)Z/(MR?+m)’, 


f/ = STATIC FORCE ON BOW TIP 
#= DYNAMIC FORCE ON BOW TIP 


FORCE IN POUNDS 





° < 8 12 6 20 24 26 
ORAW IN INCHES 
Fic. 3. 


A=450,800CA R/(.\R?+m) 
+7,868mC(A” —A?)Z/(MR?+m)?, 


f=2.22maxX10-*, F=4.44VAX10~°, 
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T=F/2 sin E=SF/2P=2.22SMAxX10-*/P 
= (f—f), sin (A+E), 
R=V/v=cosA+(Y/P) sin A. 


In order to show how these formulas may be 
used, they are now applied to a six-foot bow 
which has an eight-inch rigid middle section 
(L=4 inches), and which has a six-inch bracing 
height (J7),=6 inches), and a draw D=27.5 
inches. 


F = STATIC FORCE ON ARROW 
F = DYNAMIC FORCE ON ARROW 


FORCE IN POUNDS 





° a 8 12 '6 20 24 26 
DRAW IN INCHES 


Fic. 4, 


From the equation //)>=(3B,/4) sin Ao, Ao 
= 14°—29’, 

From the equation S=(3B,/4) cos Aj+B,/4 
+L, S=35.2373 inches. 

The active bending portion of each limb, 
B,=32 inches. 

Let each limb have a constant thickness, 
‘=0.6 inch. Let the width at the dip equal 
W =1.5 inches. Let the width at any other point 
equal w= Wx/B,, where x is the distance from 
the nock. This type bow is discussed in the 
January, 1932 issue of Ye Sylvan Archer. All 
sections of such a bow are equally stressed and 
the limbs bend in true arcs of circles. 

It may be shown that the kinetic energy of a 
limb for this type bow K=tWdB,V?/60 where 
d is the density of the bow wood and V is the 
velocity of the nock as it passes through its 
neutral position. 

If a bow limb having the same kinetic energy 
has all its mass m/2 concentrated at the nock, 
then: m V?/4=tWdB,V?/60 
and m=tWdB,/15. 


For m in grains, t, W, and B, in inches and d 
in pounds per cubic foot, this equation becomes : 


m=0.270tWdB,. 
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C STRING TENSION 
C STRING TENSION 





+ a 2 6 20 24 


ORAW IN INCHES 


Fic. 5. 
For yew wood, d=43. Therefore: 
m=0.270X .6X1.5 X43 XK 32 = 334.4 grains. 


If the arrow weighs 350 grains and the string 
weighs 105 grains, M will equal 385 grains, 
because it has been shown in previous publica- 
tions that approximately 4} the weight of the 
string should be added to the weight of the arrow 
in order to obtain its effective mass M. 

It is now necessary to determine the value of 
the constant C. 

In the January, 1932 issue of Ye Sylvan 
Archer, it was shown that the deflection at the 


nock of this type bbw N=6fB,*/WtY where Y 


is Young’s modulus for bending. But N= B,rA 
240 where A is expressed in degrees and 


f=CA. 
Therefore: B,rA/240=6CAB,*/WPhY 
C= Wt*x Y/1440B,’. 
Y=1.46X 10°. 
C= 1.005. 


from which 
For yew, 
Therefore 
The formulas developed in this paper are now 
applied and the results plotted in graph form. 


Fig. 3 shows the static and dynamic forces on 
the bow tip or nock, as a function of the draw. 
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The curve marked f represents the force tangent 
to the path required to deflect the tip for all 
positions of draw. When the string is released 
with the arrow in place, most of the potential 
energy of the bow is used in accelerating the 
arrow so that there is initially only about five 
pounds of force effective in accelerating the bow 
tip. This dynamical force f becomes smaller as 
the arrow is being discharged. It reaches zero 
for a position corresponding to a draw of about 
15.5 inches. This force then becomes negative. 
This means that the static force tending to 
move the tips forward is not as great as the 
retarding force resulting from the string tension. 
This difference is made up by the kinetic energy 
that has been acquired by the tips. The force 
which is tending to stop the tips is now being 
used to accelerate the arrow. 

Fig. 4 shows the static and dynamic forces on 
the arrow. Here we see that at the moment of 
loose, the dynamic force F on the arrow is 
smaller than the static force F which was on 


N FEET PER SECOND 


V=VELOCITY OF APROW 
v @VELOC/TY OF BOW TIPS 





a 8 te 6 20 24 26 
ORAW IN INCHES 


Fic. 6. 


the holding fingers. This is due to the fact that 
some force is required to start the limbs in 
motion. As the arrow moves forward, the force 
on it decreases as shown by the curve F. The 
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dynamical force F does not decrease as rapidly 
as the static force F. This is due to the kinetic 
energy of the tips. This explains why massive 
bow tips do not reduce the velocity of the arrow 
by a very great amount. The energy required to 
put the tips in motion is given to the arrow at a 
later time. However, this increase in the force 
F due to the kinetic energy of the limbs may not 
be desirable from the standpoint of accuracy. 

The static string tension T is given in Fig. 5. 
In the same graph we find that as the arrow is 
loosed, the dynamical tension T drops slightly 
but soon begins to increase. Were it not for the 
fact that the string may stretch slightly, this 
tension would break the string. Heavy tips 
increase this tension. 

Fig. 6 gives the velocity of the arrow and the 
velocity of the bow tips as a function of the draw. 
We see that the velocity of the bow tips becomes 
a Maximum corresponding to a draw of about 15 
inches. It then becomes smaller and smaller, 
reaching zero at the time the. string reaches its 
neutral position. Because of stretching of the 
string this zero velocity is reached slightly later 
than shown. However, in this case the tips 
actually reverse their direction and move back 
again. A heavy string accentuates this effect. 
The velocity of the arrow increases throughout 
the time of discharge. However, for light arrows 
and heavy tips, the rate of increase may become 
much less during the middle portion of the 
discharge and then increase quite sharply near 
the end. 

Fig. 7 shows the accelerations of the arrow and 
of the bow tips as functions of the draw. The 
acceleration of the bow tips becomes smaller and 
smaller until it reaches zero at a position corre- 
sponding to a draw of about 15.5 inches. The 
acceleration then becomes negative. This means 
that the tips are slowing up. The acceleration of 
the arrow decreases for a short distance and then 
the decrease becomes less rapid. For very light 
arrows it may increase again. The shorter the 
bow, the more the reflex, the heavier the tips 
and the lower the stringing height, the more 
pronounced is this effect. That is why so many 
archers condemn all these conditions where 
accuracy is desired. A bow having a permanent 
set will reduce this effect. Heavy arrows will also 
prevent the acceleration from increasing near 
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the close of the discharge. A whip-ended bow 
will also reduce the acceleration near the end of 
the discharge. This explains why a whip-ended 
bow does not kick. 

The computations given here are for a bow 
having an efficiency of 100 percent. If the string 


A= ACCELERATION OF APROW 
@ = ACCELERATION OF BOW TIPS 


ACCELERATION IN FEET PER SECOND PER SECOND 





° 7 8 12 16 20 24 26 
DRAW IN INCHES 
Fic. 7. 


is taken into consideration as in the above 
computations, bows may be built that have an 
efficiency of 93 percent. An Osage bow of the 
type described above was built by Dr. P. E. 
Klopsteg which, after taking the string into 
consideration, had an efficiency of 93 percent. 
Most bows have much lower efficiencies. If the 
bow wood is good and no part is stressed beyond 
its elastic limit, the efficiency will be, lowered 
only by the fact that the kinetic energy of the 
limbs is not all transferred to the arrow. The 
loss within the wood and in air resistance is 
usually negligible. If, in making the computa- 
tions, the final velocity of the arrow is taken for 
a position about 2 inches in advance of the 
bracing height, a value will be obtained that 
compares well with experimental results. In 
other words, in practice all of the kinetic energy 
of the tips is not converted into arrow velocity 
as assumed in the computations. 
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National Research Council Fellowships 


National Research Council Fellowships in physics for 
1937 have been awarded to the following: 
PROPOSED INSTITUTE FOR STUDY 


Cornell University 
Harvard University 


Emit JOHN KONOPINSKI 
ALFRED Otto CARL NIER 
HupertT Jos& YEARIAN Harvard University 

ERWIN RuDOLF GAERTTNER California Insiitute of Technology 
Jesse LEONARD GREENSTEIN Yerkes Observatory 

WILLIAM CONYERS HERRING Massachusetts Institute of Technology 
WILLIAM RUDOLPH KANN! University of Wisconsin 

ALEXANDER S. LANGSDORF, Jr. University of California 

JOHN REGINALD RICHARDSON Princeton University 

LEONARD ISAAC SCHIFF University of California 
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$1000 American Chemical Society Prize 


@ The $1000 prize of the American Chemical Society 
for 1937 has been awarded to Dr. E. Bright Wilson, Jr., 
assistant professor of chemistry in Harvard University. 
The award, bestowed annually upon a scientist under 
thirty-one years of age and of unusual promise, goes to 
Dr. Wilson, who is twenty-nine, for his experimental 
work in physical chemistry. 

The prize will be formally presented at the ninety-fourth 
meeting of the society to be held in Rochester, N. Y., 
September 6 to 10. Maintained by Dr. A. C. Langmuir of 
Hastings-on-Hudson, N. Y., and his brother, Dr. Irving 
Langmuir, Nobel prize winner in chemistry, it provides 
for ‘recognition of the accomplishment in North America 
of outstanding research in pure chemistry by a young man 
or woman, preferably working in a college or university.” 


* 


Science on the Radio 


@ A description of the grinding of the world’s largest 
telescopic mirror and a demonstration of the operation of 
gigantic dynamos which create artificial lightning were 
highlights on a broadcast from the California Institute of 
Technology, at Pasadena, on Thursday, April 29, from 
8:00 to 8:30 p.m., EDST, over the Blue Network of the 
National Broadcasting Company. 

Dr. Carl D. Anderson, winner of the Nobel prize, was 
among the experts who gave listeners actual demonstra- 
tions of tests conducted in wind tunnels, high tension 
laboratories, and the optical laboratory. From the optical 
laboratory, under the supervision of Dr. Max Mason, 
listeners heard the actual grinding of the world’s largest 
telescopic mirror, a 200-inch glass. Professor Clark Milli- 
kan, head of the aeronautical division, described the 
aeronautical laboratory. 
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Physics Teachers of Western Pennsylvania and Environs 


@ The Association of Physics Teachers of Western 
Pennsylvania and Environs met May 15, at St. Vincent 
College near Latrobe, Pa. The program was as follows: 
Welcome. Rr. Rev. ALrrep Kocn, O.S.B., Archabbot. 

What is considered an adequate lecture in physics? W. N. Sr. Peter, 
University of Pittsburgh. 

Apparatus for automatically demonstrating the appearance of various 
transparent substances in plane polarized light (with demonstration). 
J. B. NATHANSON. AND C. W. PRINE, Carnegie Institute of Technology. 

Supply and maintenance of equipment. FATHER BERNARD, O.S.B., St. 
Vincent College. 

The reflecting telescope as an undergraduate project. W. P. CUNNING- 
HAM, State Teachers College, California, Pa. 

Some characteristics of good definitions. A. G. WortnHinc, University 
of Pittsburgh. om 

Physics in industry. Davip J. Gites, Works Manager, Latrobe Electric 
Steel Company. 

Physics and steel. LEoNarD GRIMSHAW, Latrobe Electric Steel Company. 

Energy loss in hysteresis, in terms of the recent ferro-magnetic theory. 
G. A. Scort, University of Pittsburgh. 

Laboratory experiments. F. A. Mo.sy, West Virginia Universily. 


* 
Award of the Acheson Medal 


@ The Electrochemical Society has just awarded the 
Acheson Medal and $1000 prize to Dr. Frederick Mark 
Becket, president of Union Carbide and Carbon Research 
Laboratories, Inc., and vice president of the Electro 
Metallurgical Company and the Union Carbide Company, 
for his outstanding contributions to electrometallurgy. 
The highly developed alloy steel industry of today, with 
its stainless alloys, its high speed tools, and high tempera- 
ture resistant products used in the electrical industry, owes 
much of its rapid spread to Becket’s inventions. 


* 


Fellowship in Electrochemistry 


g The Electrochemical Society has just awarded the 
ninth Weston Fellowship of $1000.00 to Mr. Garth L. 
Putnam, of Seattle, Washington. Mr. Putnam will con- 
tinue his research at Columbia University, investigating 
methods of efficiently depositing alloys of active metals 
from aqueous solutions. The work will be done under the 
direction of Professor Colin G. Fink. 

The Fellowship was founded by Dr. Edward Weston, 
a pioneer in electrochemistry. 


* 


Behind the Scenes in Industrial Europe 


@ The European laboratory tour party of the National 
Research Council sailed May 19 on the S.S. Normandy. 
The party numbered about forty. They will spend two 
weeks visiting laboratories in England and about two 
weeks each in France and Germany. 
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University Department Changes 


Q In recognition cf steadily expanding scientific knowl- 
edge and the prospect of increasingly important develop- 
ments in industry, the Massachusetts Institute of Tech- 
nology has announced the division of its Department of 
Mining and Metallurgy into two distinct fields of pro- 
fessional education. 

Beginning next autumn the Department of Mining 
Engineering will be administered as a separate department 
under the direction of Professor W. Spencer Hutchinson, 
who has long been head of the courses in mining and 
metallurgy. His department will also include the course in 
petroleum production. The new Department of Metallurgy 
will be directed by Dr. Robert S. Williams, for many years 
professor of physical metallurgy, and will include the 
course in ceramics. 

The new Department of Metallurgy will give special 
attention to physical and process metallurgy, as well as 
metallurgical preduction, the latter being concerned pri- 
marily with the economical and statistical functions of the 
industry. The two former involve application of basic and 
scientific principles to the production and treatment of 
metals. 


* 


@ At the College of Engineering of the Carnegie Institute 
of Technology, courses in chemical engineering will be 
administered next autumn under two separate departments 
instead of under one as at present. Dr. Thomas R. Alex- 
ander, associate professor of inorganic chemistry, will 
become acting head of the department of chemistry, and 
Dr. Warren L. McCabe, professor of chemical engineering, 
will be head of the department of chemical engineering. 


* 


The Homes of Tomorrow 


@ A conference on housing was held at the Massa- 
chusetts Institute of Technology in Cambridge, June 7. 
The program arranged for this meeting was as follows: 


Vannevar Bush, Vice President and Dean of Engineering, M.I.T., 
presiding and summing up at the end. 

Sir Raymond Unwin of London, England, world-famous architect and 
town planner, on “‘Better Homes and Neighborhoods for All—The 
Social and Economic Aspects of Shelter.”’ 

John Ely Burchard, Vice President of Bemis Industries, housing re- 
searcher, on “‘How Better Homes Will be Built—The Question Mark 
of Prefabrication: New Materials, New Methods of Design.”’ 

Ernest J. Bohn, sometime President, National Association of Housing 
Officials, consultant on low-cost housing, on ‘‘Government Inter- 
vention in Housing—What Governmental Agencies Have Done and 
Might Do.” 

Robert D. Kohn, distinguished New York architect, former Director, 
Housing Division of Public Works Administration, on ‘‘The Future of 
Housing—Some Possible Ways Leading to Better Housing.”’ 


* 


Q Dr. William H. Doherty, of the Bell Telephone 
Laboratories, has been awarded by the Institute of Radio 
Engineers the Morris Liebman Memorial Prize for 1937, in 
recognition of his work in the field of radiofrequency power 
amplifiers. The presentation was made at the annual dinner 
of the institute in New York on May 12. 


Science, May 21, 1937 
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New Appointments and Promotions 

@ At Massachusetts Institute of Technology: 

John Wulff from assistant professor to associate professor. 

Wayne B. Nottingham from assistant professor to associate professor. 

Nathaniel H. Frank from assistant professor to associate professor. 
At Princeton University: 

Gaylord P. Harnwell from assistant professor to associate professor. 
At Cornell University: 


Hans A. Bethe from assistant professor to professor. 
H. E. Howe from assistant professor to professor. 
R. F. Bacher from instructor to assistant professor. 


* 


Q The degree of doctor of laws will be conferred in June 
by the University of Glasgow on Dr. Max Planck, emeritus 
professor of theoretical physics at the University of Berlin. 


* 


Q Professor Linus Pauling has been appointed to 
succeed the late Professor Arthur A. Noyes as director of 
the Gates Chemical Laboratories of the California Institute 
of Technology. He has also been appointed George Fisher 
Baker lecturer in chemistry at Cornell University for the 
first term of the year 1937-38. 


* 


Q Professor W. L. Bragg, Langworthy professor of 
physics in the University of Manchester, has been ap- 
pointed director of the British National Physical Labo- 
ratory. 


* 


Q Dr. Victor Hicks has joined the staff of the Westing- 
house X-Ray Company, Inc., at Long Island City, New 
York. 

* 





Journal of Applied Physics 
For July 


Science, Invention and Society, by WALDE- 
MAR KAEMPFFERT 


Physics in the Paper Industry, by E. W. 
SAMSON 


Physics of Air Hygiene, by H. B. MELLER 
AND F. F. RUPERT 


Physics of Subsidence and Ground Movement 
in Coal Mines, by H. P. GREENWALD 


Trilinear Coordinates, by HOWARD AIKEN 


Impregnation Studies with Colloidal Graphite, 
by BERNARD H. PORTER 


Contributed Research Papers on Mecha- 
nisms for the Starting Potentials of Corona 
Discharge, Penetration Equipment for Esti- 
mating Viscosities, and Elastic Stretch and 
the Infrared Spectrum of Rubber. 
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q In response to a request for an opinion on “The Place 
of Science in National Planning,”’ by Senator Robert J. 
Bulkley, chairman of the United States Senate Committee 
on Manufactures, the American Institute sponsored a 
discussion at the Chemists’ Club, New York City, on 
April 22. Dr. Willis R. Whitney, of the General Electric 
Company, Waldemar Kaempffert, science editor of The 
New York Times, Dr. George Baehr, trustee of the New 
York Academy of Medicine, Julius H. Barnes, formerly 
chairman of the U. S. Chamber of Commerce, and others 
were present and led the discussion. 


Science, April 23, 1937 
* 


Announcements 


Q Sir William Bragg will deliver the forty-third James 
Forrest lecture before the Institution of Civil Engineers on 
April 27, at 6 p.m. The title of Sir William’s lecture will be 
“The Crystal and the Engineer.” 

Q Professor J. Chadwick will deliver the twenty-eighth 
Kelvin lecture before the Institution of Electrical Engineers 
on April 22 at 6 p.m. The subject of the lecture will be 
“Elementary Particles of Matter.” 


Nature, April 17, 1937 
* 


Electrochemical Society Elects New Officers 


@ At the annual meeting of the Electrochemical Society 
held at Philadelphia, Pennsylvania, the following new 
officers were elected: 


President; WiLtiAM G. Harvey, Aluminum Company of America, 
Cleveland, Ohio 
Vice Presidents; L. D. Vorce, New York City 
R. L. BALDWIN, Niagara Falls, New York 
O. W. Storey, Chicago, Illinois 
Managers; W. W. Winsuip, New York City 
E. M. Baker, Ann Arbor, Mich. 
S. SWANN, Jr., Urbana, Illinois 
Treasurer; ROBERT M. BuRNs, 463 West Street, New York City. 
Secretary; Coutin G. Fink, Columbia University, New York City. 


* 


A. S. H. V. E. Semi-Annual Meeting 


@ The American Society of Heating and Ventilating 
Engineers will hold their semi-annual meeting at the New 
Ocean House, Swampscott, Mass., June 24-26, 1937. The 
following technical papers are scheduled: 


Adaptability of pre-cooling coils to air washer systems. JouN Everetts, 
Jr. 

Degree days for summer air-conditioning. A. D.- Marstons 

Study of summer cooling in the research residence using water at 
temperatures of 52° and 46°F. A. P. Kratz, S. Konzo anv E. L. 
BRODERICK 

Cooling load analysis of a bank building. Joun James 

Ventilation requirements, Part II. C. P. YaGLou anp W..N. WitHe- 
RIDGE 

Performance of fin tube units for air cooling and dehumidifying. 
G. L. Tuve Ano C, A. MCKEEMAN 

Cooling and heating rates of a room with different types of steam 
radiators and convectors. A. P. Kratz, M. K. FAHNESTOCK AND E. L. 
BRODERICK 

Report of A.S.H.V.E.-A.S.K.E. Committee on National Standards for 

ir-Conditioning Applications. L. A. HARDING, Chairman 

«oe method of testing and rating air cleaning devices. ARTHUR 
NUTTING 

Direct reading effective temperature indicator. J. R. PARSONS 

Infiltration effects on heating Grant building. F. C. HouGuTEN 
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@ At the annual meeting of the American Academy of 
Arts and Sciences, held on May 12 at its house, 28 Newbury 
Street, Boston, Massachusetts, it was voted to award the 
Rumford Medal to William Weber Coblentz, physicist of 
the National Bureau of Standards, Washington, D. C., in 
recognition of his investigations in heat and light, including 
the physical study of the firefly, photoelectric properties of 
materials, the radiation of stars and pioneer work on the 
temperatures of the planets. 


Science, May 21, 1937 


@ Dr. Thomas C. Poulter, director of the Armour 
Institute of Technology, Chicago, who was second in 
command of the Byrd Antarctic expedition of 1933-35, 
was presented on April 27 with a special gold medal by the 
National Geographic Society. Dr. Gilbert Grosvenor, 
president, made the presentation in the presence of Rear 
Admiral Richard E. Byrd. On the Byrd expedition Dr. 
Poulter’s work included geophysical investigations, studies 
of ice conditions, magnetic soundings and observation of 
meteors and auroral phenomena. 


* 


@ Science Service announces that the Massachusetts 
Institute of Technology will build a ten-foot diameter 
tunnel that will simulate conditions at 35,000 feet altitudes 
and speeds of 450 miles an hour. Similar equipment is 
already being built at Langley Field by the National 
Advisory Committee for Aeronautics. Both wind pressure 
tunnels follow the general plan and design of the wind 
tunnel of the Reich Air Ministry at the University of 
Goettingen at Hannover, Germany. The rarefied atmos- 
phere and accompanying low air pressures encountered 
in the sub-stratosphere, where aviation now would like 
to do its flying, render inaccurate some experimental 
findings obtained in previous wind tunnels at ordinary 
atmospheric pressure. 


* 


@ The Research Foundation of Armour Institute of 
Technology announces the establishment of the Eli Wish- 
nick Fellowship to promote fundamental research in the 
field of extremely high pressure: This fellowship carries a 
stipend of Seven Hundred Fifty Dollars a year. The 
Fellow may work for either an M.S. or a Ph.D., although 
preference will be given to applicants having a master’s 
degree, and intending to work for a doctorate. The appoint- 
ment will be made for a two-year period beginning Septem- 
ber, 1937. Application should be made to Thomas C. 
Poulter, executive director of the Research Foundation of 
Armour Institute of Technology, Thirty-three Hundred 
Federal Street, Chicago, Iilinois. 


* 
Necrology 


@ Alexander Howard McConnell, founder, and chairman 
of the board of directors of Central Scientific Company, 
died at Albuquerque, New Mexico on May 4, 1937, in his 
sixty-first year, while en route to Chicago from California. 
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Electrical Engineering in Radiology. L. G. H. SARSFIELD. 
Pp. 284+-xiv, Figs. 199. Instruments Publishing Company, 
Pittsburgh, 1936. Price $6.00. 

This is a fairly rapid survey of x-ray apparatus from 
the point of view of the design engineer. The contents are 
a mixture of design data, circuit information, practical 
techniques in operation and control, and general descrip- 
tions of commercial apparatus which should prove useful 
to those wishing to be oriented in some of the problems, 
methods, and apparatus of the conventional x-ray field. 
There are chapters on the design of x-ray potential and 
filament transformers, on mechanical rectifiers and elec- 
tronic-tube rectifying circuits, on x-ray tube current and 
voltage measurements, high voltage cables, voltage and 
current control apparatus, and on equipment layout. 
The chapter surveying standard x-ray tubes includes 
mainly general descriptions and illustrations of commercial 
tubes. The book concludes with brief chapters on main- 
tenance of x-ray equipment, special equipments, and rules 
of electrical safety. 

Electrical Engineering in Radiology may prove of interest 
to roentgenologists and technicians in charge of x-ray 
equipment, to prospective purchasers of equipment, and 
to students. It is hardly sufficiently detailed or technical 
for reference use by engineers and physicists already in 
the field. The author has been connected with the electrical 
development of radiological equipment at the Woolwich 
Arsenal in England, and it is therefore quite natural that 
the emphasis should be on English and continental appa- 
ratus and procedures. It will perhaps be regretted by many 
that the subject matter has been restricted to the power 
supply, circuit, and practical engineering aspects of x-ray 
production, and that such important matters as x-ray 
phenomena, intensity and quality determinations, filtra- 
tion, electronics of x-ray tubes, and x-ray vacuum tech- 
nique, which might quite properly come under engineering 
in radiology, are not considered. The material covered, 
however, is well written with many illustrations, and 
supplies much practical information. 

Joun G. Trump 
Massachusetts Institute of Technology 


Dynamics, Part II. A. S. Ramsey. Pp. 344+xi, illus- 
trated. Cambridge University 1937. 
Price $4.50. 

This book is a continuation of an earlier volume under 
the same title by the same author and deals with the more 
advanced parts of dynamics. Both are written in true 
Cambridge style, the first volume being intended for 
“first year honours’ students” and the second part hav- 
ing for its object a three years’ “honours course’’ and the 
famous “Tripos examinations.” 


Press, England, 


The author is president 
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of Magdalene College in Cambridge, and in past years has 
prepared for publication a Treatise on Dynamics by Dr. 
Besant, from which the present book has drawn some 
material. 

Coming from these surroundings it is considerably more 
difficult and thorough, and contains much more subject 
matter than the usual textbook of American colleges. 
Comparing its scope and difficulty with other English 
books on the subject, it seems to stand somewhere between 
the works of Lamb and Routh, and, like those classics, 
contains a large number of problems and 
examples. 

The book opens with some special topics in particle 
dynamics, including theory of orbits and rotating tubes, 
followed by a short chapter on motions of chains. After a 
brief exposition of kinematics various motions of particles 
in three dimensions are taken up. This is followed by two 
chapters on 
Next including impulsive 
motions, are treated extensively by means of Newton's 
equations. 

The next 
brings Lagrange’s and Hamilton’s theorems with many 
examples. These methods are applied in the following two 
chapters to small oscillations and gyroscopic motions. An 


worked-out 


two-dimensional motion of rigid bodies. 


three-dimensional motions, 


large chapter on generalized coordinates 


appendix deals with vector analysis and its applications, 
of which the author states that “‘it is hardly possible to 
ignore the present day fashion of using vectors in dynamics, 
I doubt whether it is a fashion which has come to 
stay and I am unwilling to require all readers of this book 
to use vectors.” 

Summarizing, it may be said that the book is clearly 
written, contains a large number of good problems, and, 
with its first part, will be useful for those who want to 
study the subject thoroughly, but who hesitate to take up 
the more extended works of Routh or Rayleigh. 

J. P. Den HartoG 
Harvard University 


but 


° 

IX* Congrés International de Photographie. Pp. 852. 
Published by L. P. Clerc, La Revue d’Optique Théorique 
et Instrumentale, Paris, 1936. 

Representing a cross section of the photographic art at 
the time (1935), the Proceedings of the Ninth International 
Congress, held in Paris, are worthy of perusal by all who 
use photography as a scientific implement. American 
readers will be interested in this record of the first public 
European exhibition of Kodachrome film and of photo- 
graphs taken with the aid of Pola screens. Three papers 
are concerned with the sensitometry of bichromated 
gelatin, a subject on which little has been published. 
Theorists will be pleased to note the revival of interest in 
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the theory of photography that has taken place since the 
Dresden Congress in 1931. 

Of the papers having the most practical implications, it 
was invidious to select. As an indication of the scope of 
these papers, there might be mentioned the regeneration 
of M-Q developers (p. 219), 
plates totally devoid of gelatin (p. 309), a sensitometric 


the making of Schumann 
study of reducers (p. 365), the use of ethyl bromacetate 
in the developer to produce bluish tones on chloride papers 
(p. 388), a pessimistic study of fine-grain developers 
(p. 396), the application of the Herschel effect to the 
printing of negatives having excessive contrast (p. 438), 
an inertialess light valve for sound recording based on the 
Kerr effect (p. 717), a photographic method of reading a 
legend that has been obliterated with ink of a different 
color (p. 727), and a survey of the problem of photograph- 
ing paintings in monochrome (p. 768). Those who wish to 
amuse themselves by sensitizing their own paper will be 
interested in formulas given on pages 372-386. Unusual 
applications of photography in paleontology are described 
(p. 737), as well as x-ray motion pictures of the vocal 
cords, with sound accompaniment (p. 753), ‘qui provoqua 

. . les rires de l'assemblée.” 

Unlike the Dresden Congress, which standardized the 
type of light source to be used in sensitometry, the Paris 
Congress effected no standardization. For speed measure- 
DIN 


system (unfortunately not described in the Proceedings) 


ment, the German (Deutsche Industrie Norm) 
was the leading contender, but it was rejected, pending 
further investigation, because the required development 
conditions result in too high a gamma. The counter 
proposals of the other delegations and the arguments for 
and against the DIN system give the reader an idea of 
the extent to which a single quantity will indicate the 
sensitivity of a material in practical use. 
FRED H, PERRIN 


Eastman Kodak Company 


The Analysis of Engineering Structures. A. J. S. 
PipPARD AND J. F. BAKER. Pp. 553+ix, illustrated. 
Longmans Green, New York, 1937. Price $9.00. 

Both authors are professors of English engineering 
schools and both have been taking an active part in the 
work of the Steel Structures Research Committee of the 
Department of Scientific and Industrial Research. The 
book contains not only material which is usually given to 
students in courses on theory of structures, but also some 
new material representing the research work of the authors. 
The content of the book is briefly as follows: 

The first three introductory chapters contain a review 
of methods of analysis of simple trusses and methods of 
calculating stresses and deflections in straight and curved 
beams. The 4th chapter deals with general theorems 
concerning the analysis of statically indeterminate struc- 
tures. Chapters 5 and 6 give applications of these theorems 
in calculating deflections and stresses in redundant frames. 
Chapters 7 and 8 deal with struts and beam submitted 
and lateral loads. 


to the simultaneous action of axial 


Various column formulae are also these 


chapters. In chapters 9, 10, and 11 systems with rigid 


discussed in 
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joints are considered and the methods of calculating 
secondary stresses in trusses are discussed. Chapters 12 
and 13 deal with elastic arches and with reinforced concrete 
In chapters 14 and 15 the 


determinate 


beams. influence lines for 


statically and for redundant frames are 
shown, and the impact effect on bridges is discussed. 
From the remaining portion of the book we will mention 
chapter 19, treating the design of steel framed structures, 
in. which Steel 


Structures Research Committee are discussed. At the end 


the preliminary investigations by the 


there are six more chapters briefly discussing suspension 
bridges, curved beams loaded normally to the plane of 
curvature, flat slabs, gravity dams, earth pressure, and 
the mechanical methods of stress analysis. 

In several places new material can be found which 
usually is not given in textbooks. In chapter 4, for instance, 
a method is developed for design of redundant frames 
which makes possible the determination of cross-sectional 
area of all members so as to make the stresses equal to the 
allowable stress. This method is more direct than the usual 
trial and error method. In chapters 9 and 10 the methods 
of successive approximations are discussed at some length. 
Essentially, the work of American engineers is given in 
the book and the reader may obtain an impression that 
these methods represent very recent developments. In fact, 
the originator of these methods was O. Mohr, who used 
such a method for calculating secondary stresses.! The 
method to numerical examples in 


same was applied 


Waddell’s book? and in this way great simplification in 


calculating secondary stresses was obtained. The adapta- 


tion of the method to the analysis of rigid frame structures 
such as we encounter in the design of office buildings was 
made by K. A. CaliSev.* He developed, also, a very useful 
method of analysis of Vierendeel trusses.‘ The recent 
development of the method of successive approximation is 
similar in principle to the solution used by CaliSev. In 
chapter 14 new material regarding impact effects on 
bridges is given and vibration of bridges produced by 
balance-weights is discussed. The new material in chapter 
19 has already been mentioned. 

The authors refer principally to the literature in the 
English language and do not always give proper reference 
to the work done on the Continent. In discussing bending 
of beams, for instance, they refer to W. H. Macaulay’s 
paper, although the method was originated with Saint 
Venant. In presenting the generalized theorem of three 
moments, they mention the paper by Arthur Berry, in 
1916, although similar equations were obtained by several 
authors; indeed, in 1916 these equations were already in 
textbooks. Their discussion of beams curved in plane 
gives reference to the work by A. H. Gibson and E. G. 
Ritchie, although problems of this kind were discussed a 
Saint Venant. In fact, 


broad literature on this subject, principally in the German 


long time ago by there exists 
language. 

It is the opinion of the reviewer that the book is not 
complete enough to be used as a text in American engi- 
neering schools. There is very little given about methods of 
analysis of trusses and hinged three-dimensional systems. 
Practically nothing is said about the design of Voussoir 
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arches. All recent work on the choice of the center line of 
arches, the effect of compressive forces on bending of 
arches, the elastic stability of arches are not mentioned. 
The thin walled structures such as are 
encountered in airplane design is also omitted. At the same 


discussion of 


time, the book contains some new material and may be of 
interest to engineers who would like to know about the 
new developments in the theory of structures as accom- 
plished in England. 
S. TIMOSHENKO 
Stanford University 
1 See Mohr's book, Abhandlungen (1906), p. 429. 
2J. A. L. Waddell, Bridge Engineering (1916). : 
’ Techni¢ki List, Zagreb, 1922 and 1923. See also the recent paper of 


the same author in Memoirs of the International Association for Bridge 
and Structural Engineering, Vol. 4 (1936). 

See K. A. CaliSev, “Die Berechnung der Rahmentrager,”’ Der 
Bauingeniem, 1922. See also L. C. Maugh, Memoirs of the International 
Association for Bridge and Structural Engineering, Vol. 3 (1935). 


Katalytische Umsetzungen. W. FRANKENBURGER. Pp. 
444+ xi. Akademische Verlagsgesellschaft, Leipzig, 1937. 
Price, bound, RM 36. 


Lehrbuch der Hockfrequenztechnik. F. Vi_Bic. Pp. 
775+xx, Figs. 801. Akademische Verlagsgesellschaft, Leip- 
zig, 1937. Price, bound, RM 32.80. 


The Logical Structure of Science. A. CoRNELIUS BEN- 
JAMIN. Pp. 344. Kegan Paul, Trench, Trubner & Co., Ltd., 
Broadway House, London, 1936. Price 10s. 6d. 


Dielectric Phenomena in High Voltage Cables. D. M. 
Roptnson. Pp. 173+-xii, Figs. 96. Instruments Publishing 
Company, Pittsburgh, Pa., 1936. Price $5.00. 
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was Professor of 
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Contributed Original Research 


Double Refraction in Colloidal Solutions 


H. W. FARWELL 
Columbia University, New York, N. Y. 
(Received March 31, 1937) 


HE phenomenon of double refraction during 

the flow of certain colloidal solutions is well 
known.' In these solutions it is fairly clear that 
the colloidal particles must not be considered as 
amorphous spheres but as minute crystals which 
are themselves doubly refracting, and that the 
orientation in the stream is evidence of the more 
or less needle-like form of these crystals. The 
first evidence of the sort was obtained by Kundt? 
who refraction in 
the rotating cylinder method of 
Maxwell.‘ Kundt’s results indicated a differentia- 
tion between double refraction which is pro- 
portional to rate of slip and double refraction 
which is dependent upon orientation of long 
colloidal particles. Whether or not the two are 
produced by the same cause, an actual shearing 


was investigating double 


liquids by 


stress, it has seemed to be desirable to ascertain 
whether in reality there were two effects or 
only one. 

If the effect is single, then double refraction is 
definitely related to rate of slip, and becomes 
zero when the rate of slip is zero. As a first 
attempt to secure the proper conditions a thin 
stream of a vanadium pentoxide solution was 
forced through a narrow slot in the bottom of a 
reservoir over which air was admitted under 
pressure. The stream passed vertically downward 
with its plane parallel to two crossed Polaroid 
plates about 6.5 cm in diameter set with principal 
directions at 45° to the direction of motion of 
the stream. 

E.g., Freundlich, Elements of Colloidal Chemistry, p. 134. 
A general summary may be found in Coker and Filon, 
Photoelasticity. 


§ Kundt, Ann. d. Physik 13, 110-133 (1881). 
* Maxwell, Collected Papers, Vol. I1, p. 379. 


1 
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While there is shearing during the passage in 
the slot, once the liquid has passed beyond the 
orifice, there can be left only such small shears 
as are produced by surface tension and by 
contact with the air. The effect of contact with 
the air is to produce a possible slip in the direc- 
tion of motion. If the plane of slip is parallel to 
the direction of propagation of light, then the 
principal stresses are normal to the direction of 
propagation of light and at 45° to the vertical, 
that is, parallel to the directions of polarizer and 
analyzer. Hence this cannot change the intensity 
of the light transmitted. If the plane of slip is 
perpendicular to the direction of propagation, 
the principal stresses are in the vertical plane 
parallel to the direction of the light. Hence a 
small effect may be due to this slip. The effect of 
surface tension is to contract the stream and is 
therefore tending to produce slip in the hori- 
zontal direction, the effect of which can be 
analyzed in a similar way. 

Actually the entire stream in the field of the 
Polaroid plates showed a brilliant fairly uniform 
color. The intensity was much greater than that 
observed when the plates were set with polarizing 
directions vertical and horizontal. The conclusion 
is that in that stream the rod-shaped particles 
were largely oriented with the long dimension 
along the direction of the stream, and held that 
position at least during the passage through the 
field. However, there is not a complete elimina- 
tion of the effect of slip. 

A second attempt to find double refraction in 
a fluid with rate of slip zero was made by study- 
ing the effect of introducing a cell containing a 
V0; sol or a benzopurpurin sol between the 
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crossed Polaroid plates. 
(The cell was about 5 





cm high, 3 cm_ wide, 
and 1 cm thick.) The 
double refraction in the 
cell is 


very plain if 


there is any flow in 


the cell, and some of 
the patterns obtainable 
quite remarkable. 
(Several are 
indicate the type of 
phenomenon observed.) 


are 
shown to 


The patterns were often 


changing rapidly, and 


even small temperature 


differences were quite 


sufficient to start con- 


vection currents in 





which long particles were 











given other than a ran- 5 
dom orientation. 

With the observation 
cell surrounded by a 
water jacket the tem- 
perature gradients could 
be eliminated. Various patterns could then be 
produced by jarring the cell slightly. Yet what- 
ever patterns were formed, all tended to persist 
long after the jarring had ceased. However, the 
patterns eventually grew less sharp, changed 
form, and almost disappeared, leaving the 
general appearance of the cell not completely 
dark but nearly uniform. For instance in a typi- 
cal case the pattern formed after the cell had been 
slightly shaken from side to side, not stirred 
internally, could be followed for over an hour. 

It is then certain that double refraction occurs 
in this case because the particles have been 
suitably oriented, and that having been so 
arranged they tend to maintain the configuration 
though ultimately the ordinary molecular dis- 
turbances will again result in bringing about a 
random orientation. 

In both cases the shearing stresses involved 
under the conditions are very small, far below 
values required to produce the double refraction 
of the type discussed by Maxwell. Hence the 
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Fic. 1. Cell after standing overnight. No water bath. Fic. 2. Lower right side 
of cell is being warmed by a jet of air on the outside. Fic. 3. Air jet discontinued. 
Fic. 4. Several minutes after No. 3. Fic. 5. Cell oscillated slightly for a few seconds. 
Fic. 6. Another view after oscillation. Fic. 7. Cell contents after stirring with a rod. 
Fic. 8. After a different stirring. 


conclusion that the double refraction observed 
by Kundt is of two kinds. 

The patterns observed in the cells are in 
themselves of no particular importance, because 
it is impossible to deduce therefrom anything 
definite as to the orientation of particles in a 
particular region, the phase difference being an 
integrated value along the light path through a 
heterogeneous medium. However, it is possible 
that the results may be of interest in studying 
some aspects of disturbances in fluids, and 
several photographs of observed patterns are 
appended. For these photographs I am indebted 
to Dr. G. N. Glasoe, who has succeeded in 
making satisfactory exposures, although both 
the colloids used, vanadium pentoxide and 
benzopurpurin, transmit only narrow bands in 
the red. 

My thanks are also due to Professor A. W. 
Thomas of the chemistry department for the 
preparation of the colloidal solutions, and for 
discussions of the problem. 
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On the Creep of Solids at Elevated Temperature* 


A. NADAI 
Westinghouse Electric and Manufacturing Company, East Pittsburgh, Pa. 


(Received March 26, 1937) 


In this paper equations are set up which describe the 
behavior of a solid under the action of combined stresses 
when a steady state of creep is assumed to exist. A material 
is said to be in a steady state of creep if it yields with a 
constant speed under a given constant stress. In an 
ordinary creep test the steady state exists when the 
deformation has reached the point where the ‘“‘creep 
curve” becomes a straight line. In mathematical language 
this can be expressed by saying that the stress is a function 
of the rate of strain and is independent of the amount of 
strain. Cases are worked out in which this stress vs. strain 


INTRODUCTION 


N account of the great variety in the 

mechanical behavior of solids and of the 
metals as engineering materials in particular the 
impossibility of ascribing to them simple general 
mechanical ideal properties such as can be 
attributed in fluid mechanics to liquids has been 
recognized. It seems that progress has rather 
been made in the opposite direction, namely by 
confining the treatment to smaller groups of 
solids, by limiting the ranges of the variables, 
such as the strains, the strain velocities and the 
temperatures and by describing the plastic 
deformations independently from the strength 
properties (for example the breaking conditions, 
the fatigue phenomena, etc., need separate 
treatments). 

It was mentioned elsewhere that the treatment 
of the plastic deformations of solids may be 
restrained either to cases of plastic flow in which 
the yield stresses remain independent of the 
speed of deformation or to cases in which they do 
not depend on the strains, but only on the strain 
rates, i.e., on the speeds with which the strains 
change with time. While the former kind of flow 
is encountered in the polycrystalline ductile 
metals at comparatively low temperatures in the 
strain-hardening range, the latter flow is charac- 
teristic of the slow creep of metals at elevated 
temperatures in the creep range. This latter 
group of the permanent deformations of solids 

* Presented at the Symposium on Metals jointly spon- 
sored by the American Institute of Physics and the 


Massachusetts Institute of Technology, Cambridge, Janu- 
ary 29, 1937. 
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rate relation is a simple power function. Various types of 
flow including purely plastic flow and viscous flow can be 
discussed by changing the value of the exponent in the 
power function relation. Several specific examples which 
include the case of the thick walled cylinder with closed 
ends loaded by internal pressure and the case of a disk 
having a circular hole and loaded in its plane are worked 
out. Equations are developed which show the steady state 
stresses in such cases. Expressions for the stress concen- 
tration factors in creep are also obtained. 


has been experimentally studied by engineers 
during the past years intensively for the ductile 
metals at elevated temperatures in connection 
with creep. While the present paper primarily 
was concerned with this latter group of phenom- 
ena and with the creep of metals under combined 
stresses, it was later found that the formulae 
which were derived for the steady state of 
distributions of stress assuming certain speed 
relations valid for creep would express also the 
distributions of stress in the strain-hardening 
range of metals under the same loading con- 
ditions for a corresponding strain-hardening law. 
This will be explained later in more detail. 

Further restrictions refer to the order of 
magnitude of the permanent strains (whether 
infinitesimal or finite) and to the type of strains, 
which have to be considered in certain cases 
(whether elastic or permanent or both). In view 
of the extremely small deformations occurring 
during creep in the following only infinitesimal 
strains and strain rates need to be considered. 
Comparatively little is known even under these 
simplifying assumptions about the redistribution 
of the stresses around holes or notches or in 
heavily stressed, thick walled cylinders or rotating 
disks when the elastic distribution of stress due 
to the exposure to high temperatures and to 
creep is gradually changing with time.! 





‘A considerable amount of valuable information on the 
creep of metals of an experimental nature became available 
during the recent years, which cannot be quoted here in 
detail. Among the papers which dealt recently with 
problems of a mechanical nature in creep, the following 
should particularly be mentioned: R. W. Bailey, ‘‘The 
Utilization of Creep Data in Engineering Design.”’ Insti- 
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Existence of a steady state of creep was 
assumed. After a certain time has elapsed, it is 
assumed that the creep rates do not change with 
time more under constant stress. The formulae 
permit to compute perhaps for the first time a 
few exact expressions for factors of stress concen- 
tration in simple cases corresponding to the 
special creep speed law or strain-hardening curves 
which were used. 

The elastic portions of the strains were neg- 
lected. It is expected that by this assumption, no 
serious error in the values of the stresses in the 
steady state of creep resulted. It should be 
expressly noted, that the formulae do not con- 
verge in the neighborhood of zero stress to the 
right expressions.” For the purpose of finding if 
possible simple expressions for some distributions 
of stress in the steady state of creep, use of 
certain mean quantities dependent on the six 
components of stress and of strain and of the 
strain rates is made. They are the measures of the 
intensities of stress, plastic distorsion and rate of 
distorsion. 

One of the simplest and at the same time most 
symmetrical expressions which offers itself for 
this purpose is the sum of the squares of the three 
principal stress differences. In the homogeneous 
stress field in the neighborhood of a point, there 
are eight symmetrically situated plane sections 
or four sets of parallel planes in which the shearing 
stress is equal to 


tn= 4[(o1—02)? 


+ (¢2—03)?+(03—0;) }} (1) 


and in which at the same time the normal stress 
a, becomes equal to the mean stress o 


on= (01 +02+03) 3=<0, (2) 


where oj, o2, 3 designate the principal stresses. 
tution of Mech. Eng., London (1935). Cf. also Engineering 
140, 595 (1935). F. K. G. Odqvist, ‘‘Theory of Creep Under 
Action of Combined Stresses with Applications to High 
Temperature Machinery.” Proc. Royal Swedish Institute 
Eng. Research, No. 141, Stockholm (1936). (A part of a 
paper communicated also to the Fourth Congress for 
Applied Mechanics, Cambridge, England, 1934. Cf. Proc. 
p. 228.) C. R. Soderberg, ‘‘The Interpretation of Creep 
Tests for Machine Design.” Trans. A. S. M. E. R. P. 58- 
15 (1936). Also, ‘‘Working Stresses,’’ same journal (1933). 
Paper APM 55-16 p. 131, and Schweizernche Bauzeitung 
104, 127 (1934). G. M. MacCullough, ‘‘Applications of 
Creep Tests” Trans. A. S. M. E. 55, 87 (1933). H. Hencky, 
‘‘The New Theory of Plasticity, Strain Hardening and 
Creep,”’ Trans. A. S. M. E. 55, 151 (1933). 

2For example the differential quotient 00/de when 
o =0, e=0 should become equal to the modulus of elasticity 
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These planes are the octahedral planes with 
reference to the sides of a cube oriented with its 
faces normally to the principal stress directions. 
If one selects a small regular octahedron as an 
element of the plastic material oriented with its 
axes parallel to the principal stress directions, the 
normal pressure o, on each of the eight faces of 
this element is equal to the mean pressure co. 
Since this latter has no influence on the yielding, 
it seems reasonable to consider the other stress 
component, namely, 7, as the measure of the 
stress intensity producing plastic flow or yielding. 
This stress r, defined by Eq. (1) is called the 
octahedral shearing stress. 

Similarly, the unit shear in the octahedral 
plane 


Yn = 3L(e— 2)? + (2-63)? + (6-41)? ]!, (3) 


where «1, €, €; are the principal strains is called 
the octahedral shear. Among the various planes 
in the neighborhood of a point in which the unit 
shear is equal to y, the eight octahedral planes 
have the additional remarkable property that the 
strain in the normal direction to them is 


€n=(ertetes)/3 (4) 


and consequently equal to the mean strain. 
Since for small strains the sum €,+«+ 6; repre- 
sents the volume dilatation and during a plastic 
deformation the volume does not change appreci- 
ably, so that this sum remains equal to zero, 
we see that 

€, = 0. 


Therefore, the distorsion of the octahedron con- 
sists only in the unit shears y, and the intensity 
of a small plastic deformation is measured by the 
octahedral shear y,. 

When a ductile metal is deformed at low 
temperatures, stress is dependent on the strains 
and must be increased for maintaining flow. 
This is expressed by assuming that 


Tr =f(Yn) (3) 


is a monotonously increasing function. (5) ex- 
presses in the case of the metals what is known as 


“cold working” or “‘strain-hardening.”’ In the 


statement (5) is implicitly contained the other 


E under a pure tensile stress o independently of the 
strain rate. The assumptions which were introduced disre- 
gard this fact. 
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observation that the speed with which the plastic 
deformations were produced was not essential for 
the magnitude of the stress. In the special case of 
(5) when a substance yields under a constant 
octahedral shearing stress, or when 


tT, = V20)/3=const. (6) 


the substance is called perfectly plastic. The 
normal stress oo is the plastic limit in pure 
tension. o9 is the normal stress under which the 
substance would yield under tension. 

At elevated temperatures, on the contrary, 
stress is more dependent on the speed of deforma- 
tion, less on the strains and one may assume as a 
second fundamental form of flow a speed de- 
pendent plastic flow defined by a relation 


Ta=f(¥n); (7) 


where y,=dy,/dt designates the rate of change 
of the octahedral infinitesimal shears y,- 7, is the 
octahedral rate of shear. 

When the stress r, is increasing in the same 
proportion as the shear rates 7, are increasing, 
the substance is called perfectly viscous: 


Ta=7'° Yn) (8) 


where 7 is a material constant. 

It is admitted that in the preceding Eqs. 
(5)-(8) and in the applications to follow, empha- 
sis is laid on ideal conditions of flow. Eq. (7) and 
some of its consequences will particularly be 
investigated in the hope of finding simple 
distributions of stress in the steady state of creep 
for a few cases of practical interest. This is done 
also in the hope of verifying some of the formulae 
derived after assuming a more general creep 
speedlaw through the two limiting special cases 
of flow, which have already been quoted referring 
to Eqs. (6) and (8), respectively. Inspection 
of Eqs. (5) and (7) expressing conditions of flow 
has shown that these two more general types of 
flow must include as special cases the perfectly 
plastic substance characterized by Eq. (6) and 
the perfectly viscous substance characterized by 
Eq. (8). If solutions based on the two more 
general Eqs. (5) and (7) are later established, 
they must contain as special cases the solutions 
which would have been obtained if they would 
have been based on the conditions (6) or (8), 
respectively. 
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It will become apparent that further con- 
nections exist between the two Eqs. (5) and 
(7) defining strain and speed dependent plastic 
flow respectively and between the well-estab- 
lished laws for a perfectly elastic substance on 
one side with the law for viscous flow on the 
other side and finally with the law for perfectly 
plastic flow. For such a comparison it is necessary 
for quoting here briefly some partially well- 
known experimental facts. 


THE STRESS STRAIN AND THE STRESS-STRAIN 
RATE RELATIONS FOR PLASTIC FLOW 


When a material element is permanently 
deformed, the rates of shear dy/0t at any instant 
are proportional to the corresponding shearing 
stresses 7, independently of the orientation of the 
planes of shear or of slip. This is expressed by 


Oy /dt=T/n, (9) 


valid for any direction of the planes of shear. 
The corresponding is true for the unit shears + 
in an elastic material : 


y¥=1/G (10) 


with the difference that the proportionality 
factor G (modulus of rigidity) is a material 
constant. Applying (10) to the case of pure 
tension in an elastic substance, for example when 
a normal stress o, acts parallel to the axis and the 
unit strains e, in the same direction and the 
strains ¢, and e, in two perpendicular directions 
y and z are 


¢,=0,/E, e=e 


= — ve, (11) 


‘(E modulus of elasticity, v Poisson's ratio), since 


(10) must hold also for the maximum shear 
Ymax=€2—€, and the corresponding shearing 
stress 1S Tmax =0,/2 


€,—€y= (1+ »)0,/E=Tmax/G=0,/2G 
the relation for the elastic constants follow: 
E=2(1+ )G. (12) 


Consequently three normal stresses, oz, oy, ¢: 
acting in an elastic material in the directions of 
x, y, 2 produce in the x direction the strain 


€.=(¢,—vo,—vo,)/E (13) 
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and adding the three similar expressions for 
Ex, €y, €z ZiVeS 


e=3(1—2y)c/E, (14) 


where e stands for the elastic volume dilatation 
e=e,+e,+e, and o is given by (1). When v=}, 
e=0 and 
E= 3G. (15) 
Similarly, for establishing the stress rate of 
strain relations valid for a steady plastic flow or 
for creep, we note that when a normal stress ¢, 
acts, the material flows with a rate of strain 
u,=0e,/dt in the direction of the x axis and 
contracts in the perpendicular directions with 
the strain rates 


u,=de,/0t and u,=de,/dt. 


Since for small permanent strains e=e,+¢,+e«,=0 
also the sum of the three permanent strain rates 


Uztuy+u,=0 (16) 


must vanish. Therefore, a normal stress o, must 
produce the following strain rates: 


Uz=0;,/>, Uy=u,= —Uu,/2=—o,/2¢. (17) 


Since the maximum rate of shear corresponding 
to this simple distorsion is u,—u,=3u,/2 and 
Tmax =02/2= ou;/2 of Eq. (9): 
Tinax =(Uz—Uy,) and we see that for the perma- 
nent strain rates: 


on account 


= 3n. (18) 


Consequently the plastic stress-strain rate rela- 
tions when the normal stresses o,, o, and a, act 
simultaneously, must be 


u,=[o,—(o,+0,.)/2 ]/¢, 
uy =[o,—(o,+0;) 2] , 
uz=[o.—(o,+0,)/2 ]/¢. 


To these come the relations for the shearing 
stress components expressed by Eq. (9). We will 
designate in the following shear rates with the 
letter v, so that v stands for dy/dt. For the 
shearing stress therefore 


(19) 


(20) 


We note that the Eqs. (19) satisfy already the 
condition of incompressibility Eq. (16), so that 
(16) needs not be introduced more, and that the 
quantities ¢ and » which have been introduced in 
the plastic stress-strain rate relations (19) and 


v=r/n. 
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(20) are analogous to the familiar moduli of 
elasticity and rigidity E and G, respectively, 
appearing in the stress-strain relations (10) and 
(13) of an elastic substance. Since both sub- 
stances are assumed as incompressible we have 


o=3n E=3G. 


¢@ and » could be called the plastic normal and 
shearing moduli, but for other reasons be- 
coming soon apparent they have received other 
designations. 

The permanent strains ¢e,, --- and strain rates 
u,,*** to be considered in this article remain 
small. Products of these latter quantities where 
they would be necessary in dynamic equations 
can be neglected. The forces producing flow or 
creep can also be assumed to be in equilibrium. 

Under these circumstances the correspondence 
between the group of equations valid for a 
perfectly elastic substance and that for the 
steady flow of a perfectly viscous substance ¢, 
n=const becomes apparent. All that is necessary 
for establishing full correspondence is to assume 
incompressibility for the elastic material. The 
following variables and material constants corre- 
spond to each other: 


and 


Elastic substance: Strains: ¢,, €,, €:, y; Constants: E=3G. 
Viscous substance: Strain rates: uz, Uy, Uz, V; 

Constants: ¢=3n. 
The quantities ¢ and 7 in the case of a perfectly 
viscous substance are known as the coefficients of 
viscosity for tension (or compression) and for 
shear, respectively. During the slow steady flow 
of a perfectly viscous solid substance the distri- 
bution of the stresses is described by the same 
functions of the coordinates x, y, z as that which 
would be obtained in an incompressible elastic 
solid under the same boundary conditions.* In 
many cases this correspondence exists even 
though an ordinary (compressible) elastic solid is 
substituted for the incompressible one, namely in 
all those numerous cases of stress more, in which 
the functions representing the stresses are inde- 
pendent of Poisson’s ratio v. For example, the 
equation of the elastic line of a beam in bending 
represents the shape of the bent line of a beam 
of a viscous material subjected to the same 
forcefield in the steady state. Or the formulae 


3See also, T. N. Goodier, ‘‘Slow Viscous Flow and 
Elastic Deformation,”’ Phil. Mag. 22, 678 (1936). 
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for the radial and tangential stress o, and a, in 


function of the radius 7 in an elastic cylinder 
subjected to internal pressure express also the 
distribution of o, and ¢, in a similar cylinder of a 
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Fic. 1 (left). Strain-hardening curve. Fic. 2 (right). Creep 


speed law. 


steady state. In the 
disk of 


revolving with an angular speed w 


viscous material in the 


formulae for an _ elastic radius ‘‘a”’ 


o,= (3+ 7)yow"(a* —r*) / 8g, 
o,=[(3+v)a?—(1+3,r)r* }yow"/ 8g.4 


However, v would have to be taken equal to } for 
obtaining the corresponding stress formula for a 
revolving disk of a viscous material : 


o,= 7yow"(a*?—r*) / 16g, 
o,= (7a® — 5r*) you" 16g. 


From the stress-strain rate relations for speed 
dependent flow given by Eqs. (19) and (20) the 
stress-strain relations for strain dependent flow 
are obtained after substituting in (19) and (20) 
for the letters u and v the letters « and y: 


e,=[0,—(0,+¢,) 2 

¢,=[o,—(0.+0,)/2 ]/9, 
e,=[o0,—(0,+0,)/2 ] | 
y=T/n. 


When the strains, respectively, the strain rates 
are infinitesimal and equilibrium prevails, a 
strain dependent plastic distortion bears the 
same relation to incom- 
pressible elastic material than a slow steady 
speed dependent flow bears to the steady flow of a 


the distorsion of an 


viscous substance. This is perhaps illustrated bv 
comparing the graph of the strain-hardening 
curve (Eq. 5) indicated in Fig. 1 with the graph 
of the creep speed law (Eq. 7) shown in Fig. 2. 
Let us assume that both of these curves would 
correspond to the same mathematical law, i.e., 
r=fi(y) would be the same function of y as 


‘ yo is the specific weight. 
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t=f,(v) is of v. Then the solution expressing the 
stresses in function of the coordinates is the same 
for a plastic strain-hardening material as it is in 
the steady state of creep. If, for example, the 
strain-hardening law is expressed by 


T= T1(¥ ¥;)™ 
and the creep speed law by a power function 
T = T2(V/ V2)” 


after a proper selection of the constants and 
using the same exponent in both functions, the 
stresses under given boundary conditions would 
be expressed by the same functions of the coordi- 
nates. When for the exponent m is taken =1 the 
curves in Figs. 1 and 2 are straight lines, corre- 
sponding to the special cases of the elastic, 
incompressible and the perfectly viscous sub- 
stance respectively. 


Two—DIMENSIONAL CREEP FLOW 


For the permanent deformations the following 
equations are available: The strain rates u,, u,, u. 
and the rate of shear v expressed in terms of the 
components w, and w, of the velocity vector w 
of a point (x, y) in the solid are 


Ow, OW, OW, OW, 
lz=—, ty=—, v=—+ (22) 
Ox oy oy Ox 
the condition of compatibility of these three 
equations is 
Cur, Ou, dv 
EE, BaF he ETRE (23) 
Ox" Oy’ Oxdy 
the equilibrium conditions are 
Oo, OT Oo, OFT 
— +—=0, —+—=0, (24) 
Ox oy dy ax 
the plastic stress-strain rate relations are 
u,=(20,—0,—¢,)/6n, | 
z z y z 1, | 
u,=(20,—0.—¢,)/6n, ™ 
4 4 - n q (25) 


u,=(20,—0,—<0,) /6n, 


Vv=T/ YN. 


We noted already that this last group of relations 
takes care of the condition of incompressibility 
expressed by Eq. (16). Similarly as in the theory 
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of elasticity a few special cases are of importance 
in which one additional condition is imposed on 
the cinematic relations (25) which shall be 
mentioned now: 


a. Plane flow 


The normal strain rate u, in the direction of the 
z axis is constant. The simplest case is u,=0. 
From Eq. (16) follows then: 


(26) 


Uz=—Uy 


and from the third of the Eq. (25) that in this 
case the axial stress ¢, must become equal 


o,=(0,+0,) 2. (27) 


Instead of the group of Eqs. (25) the stress- 
strain rate equations are condensed to: 


Uz= —Uy=(o,—a,)/4n, v=T/n. (28) 


b. Thin plate or disk o.=0 


After inserting this value for ¢, in (25), the 
following is obtained : 


u,=(20,—0,)/6n, Uy=(20,—¢0,)/6n, v=T7/7. (29) 


If we count now the variables appearing in 
the Eqs. (23, 24, and 28 or 29, respectively), we 
find seven: o;, dy, T, 0, Uz, Uy, v, for which only six 
equations are available. The one additional equa- 
tion needed is derived from what was called the 
‘‘speed”’ or “creep law’’ or the material stating 
the physical relation which must exist between 
the speed of deformation (the creep rates) and the 
stresses under which they are produced. 


c. Perfectly viscous substance 


One special case has already been quoted, 
namely that of a perfectly viscous substance: The 
viscosity coefficient 7 in the cinematic relations is 
constant. To Eq. (22) or r= nv must be attributed 
two simultaneous meanings: it expresses the 
cinematic relation of flow and also the speed law 
in this particular case known as the law of 
viscous flow. 


d. Perfectly plastic substance. inde- 


pendent flow 


Speed 


In this case n is considered as an unknown 
function of the coordinates x, and y. One 
additional equation is still needed which is 
obtained by assuming a pure stress condition of 
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yielding. For this latter the condition is taken 
that the shearing stress 7, in the octahedral 
planes is equal to a constant or that 


? 


(01 — o2)*+ (o2—03)*+(03— 0)? = 97,7 = 20;". 


(30) 


As experiments have shown it represents fairly 
well the yielding of the ductile metals at com- 
paratively low temperatures. (01, o2, a3 are the 
principal stresses, a» is the yield stress in pure 
tension.) 

In the case of two dimensional flow of a 
perfectly plastic substance stresses o,, o,, 7 can 
now directly be computed from (24) and (30) and 
the remaining equations serve for the evaluation 
of uz, Uy, v, n. The compatibility equation shows, 
however, that in this case the u,u,v can be 
multiplied with an arbitrary factor. These 
quantities thus remain undertermined. With 
other words the velocities with which the sub- 
stance flows can be increased or diminished in the 
same proportion in all points, without altering 
the values of stresses during yielding. 






> 


CREEP 
STRAIN € 





TIME T 


Fic. 3. Steady states of creep. 


e. Steady creep obeying given law. Speed de- 
pendent flow 


In the steady stage of creep the existence of a 
function 
(31) 


was assumed in which 7, and v, are the shearing 
stress and the rate of shear in the octahedral 
planes which were defined by the expresstons 


tn? = (1/9)[(01—02)?+ (62-03)? +(¢3—01)"], (32) 
Vn? = (4/9)[ (4; — U2)? + (ue — U3)? + (u3— 1)? ]. (33) 


Eq. (31) is based on the fact that the creep time 
curves e=/,(t) derived from observation in tensile 
tests under a constant stress o=const. seem to 
flatten out (Fig. 3) so that after sufficient time 
has elapsed they might in the first approximation 
be considered to consist of a family of inclined 
straight lines e=a+dt, the slopes u=de/dt of 
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Fic. 4 (top). Fic. 5 (bottom). Stresses in thick walled 
cylinders under internal pressure. 


which when plotted against stress furnish a 
relation o=f2(u) similar to Eq. (31). This there- 
fore refers only to ‘the steady state of creep’’ and 
the following theory is not capable of predictin: 
the conditions in the earlier transition period. 


ROTATIONAL SYMMETRICAL CREEP 


a. Creep in cylinder with closed ends under 
internal pressure 

We shall consider first the case that in axial 
direction creep does not occur, i.e., that the 
strain rate in the direction z parallel to the axis of 
u,=0. the radius r as 
coordinate, o, and o,; for the radial and tangential 
stress (Figs. 4 and 5), wu, and u, for the radial and 
tangential strain rate, the cinematic relations 
according to Eq. (28) are 


rotation Introducing 


uy= —U,=(o,—<,)/4n, (34) 

the compatibility condition is 
(d/dr)(ru,) —u,=0, (35) 

and the equilibrium condition 
(d/dr)(ro,) —o,=0. (36) 


To these four equations for the five unknowns 
u,Uu.c,om Must now be added a statement about 
the speed law. We have first to compute the 
octahedral shearing stress 7, and the shear rate 
v,. Omitting the subscript when these latter 
quantities are expressed and noting that o,, 01, ¢ 
are principal stresses, “,, u:, (u,=0) principal 
strain rates and finally that on account of u,=0, 


7 
2 
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up=—uU,, and o,=(0,+0,)/2 the octahedral 


shearing stress from Eq. (32) is found equal to 


tT=(0,—<0,)/+/6 (37) 


and the octahedral shear v from Eq. (33) is found 


equal to 


v=2[2(u?+uu,+u/)/3}'=2(2/3)'-u,. (38) 


For the creep law 


°)= 


(39) 


T= r*(v U 


is introduced where r*, v* and m are constants. If 
- are corresponding values, then (39) 
can also be written as 


71V1, T2Ve, * 


(40) 


T= 17;(v/v;)" = T2(v/v2)™. 


Either of the pairs of constants: 7r*, v*, 7%, 


ae 


To2, ++ can be considered as the “‘material con- 
stants’’ of the creep law. The exponent is 
O0O=m=1. We note that when (37) is divided by 


(38) and use is made of (34) between the octa- 
hedral shearing stress tr and shear v the second 
equation 

(41) 


T= NV 


must be obtained because the cinematic relations 
must also refer to the octahedral planes of shear 
in the stress tensor. Combining now the two 
expressions (40) and (41) 7 can be expressed in 
terms of r alone 


9=1(71/7)!-», (42) 


Here +7 is replaced using (37) by the 


difference o,—<¢, 


stress 


(43) 


n, designates the value of »;=7,/v;, and the new 
exponent ” introduced in (43) is 


n=1/m—1. (44) 


When m varies between 0 and 1 the new exponent 
n varies between ~ and 0. 
On account of u,=—u, the compatibility 


Eq. (35) can immediately be integrated giving 
(45) 


u,=(o,—0,)/4n=Cc/r’, 


the integration constant c is found from the 


boundary condition 
r=7}, 


(46) 
(47) 


Or = =, 
c=(01+ p)ri?/4m. 


o1=91, 


which gives 
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Computing the value of 7; from (37) 
T1=(01+ p)//6 

and inserting 7; in (43) and c from (47) in (45) we 

obtain 


(48) 


Lo (71+?) r; 


n  fa-et 


and (49) 


These last two expressions allow to compute the 
stress difference. ¢,—o, and the viscosity coeff- 
cient in terms of r: 


o,—0;= (oitp)(r 7)~2/ (n+) | 


2n/(n+1) 
° 


a=: (r/nri) 


After using the value for the stress difference 
a:.—o, Which has just been computed, the equi- 
librium condition (36) can be integrated and 
furnishes with the second boundary condition 
r=feo, o,-=0 the unknown radial stress o,. With 
the abbreviations 


p=r/To, pri=i/fe (52) 


the following exact expressions for the steady 
radial and tangential creep stresses in a heavy 
walled cylinder subjected to an internal pressure 
p are obtained 


_ —c[p 2m 1 |, 
0 OM ii | (53) 


where C=ppr m (1 — p1 2m) | 


The tangential stresses o, at the inner and at the 
outer surface of the cylinder are found equal to: 


-o,=(2m—1+ );"")p (1 —pi"),| 


p=pi-° 5 
: 54 
‘o2=2mpi"p (1—p,;""). { (S 


p=1-- 


Some of these stress distributions are plotted in 
Fig. 6 exe for the exponents m the values: 
m=1, %, 3, } and 0. The corresponding curves 
reguenaiian the creep speed laws are given in 
Fig. 7. The stress formulae (53) when m is chosen 
equal to 1 are: 


bpr*(p* — 1) pe:*(p* +1) 
CS ae, eae, ( 
1—p; 1—p;’ 


wn 
wn 
~— 


These coincide with the well-known expressions 
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for the stresses in a viscous or an elastic tube.This 
was to be expected since when m is taken equal to 
1 according to the speed or creep law (39) 


stress r is proportional to the shear rate v and the 
case when m=1 evidently corresponds to flow of 
a tube of a perfectly viscous substance. 

When on the other hand the exponent m is 
approaching the value zero the diagram repre- 
senting the creep laws Fig. 7 shows that r tends to 
become a constant r=7,, the shearing stress 7 is 
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Fic. 6. Steady stress distribution in heavy walled cylinder 
after creep. 

becoming independent of the shear rate v-m=0 
therefore corresponds to flow in a tube of 
perfectly plastic substance. 

When m=0 the stresses appear in the form 0/0 
but their expressions can easily be derived from 
(36) by direct integration : 


In p 





o;=— 
In pi 


We see furthermore that an exponent m exists for 
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Fic. 7. Shearing stress r versus rate of shear v for various 
speed laws r-7;(v/0,;)"; m=0 perfectly plastic substance ; 
m=1 perfectly viscous substance. 


which the steady tangential creep stresses in the 
tube wall are uniformly distributed. This occurs 
when 


m=, o.= pp. (1—p;) =const. (57) 


(o, in this case can even become smaller than ), 
when p; <3). 

These three special cases of the formulae (53) 
have also been plotted in Fig. 4, which gives an 
idea of the o, curves to be expected in a heavy 
walled cylinder for the various creep laws when 
the steady state has been attained. The figure 
discloses an interesting result, namely that the 
maximum stress o; max Which in an elastic tube is 
found at the inner surface after creep took place 
may shift to the outer surface of the tube. The 
stress o; at the inner surface becomes larger than 
the stress o2 at the outer surface as long as the 
exponent in the creep law (39) is }<m<1. 
Fig. 6 for the stresses in a tube contains the two 
extreme cases of flow for a perfectly viscous 
(m=1) and for a _ perfectly plastic (m=0) 
substance. 

With regard to the axial stress o, one remark is 
still in order. We treated the case making the 
assumption that the tube yields or flows without 
changing lengths in the axial direction, i.e., when 
u,=e,=(. In this case the axial stress according 
to Eq. (27) is o.=(oe,+¢,) /2. 

Using (53) this gives 

pmp,*" 


¢,=- ’ (58) 
( 1 — pi” ) p>” 
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If now using this o, the mean value of the axial 
stress ¢, over the tube wall or 


6.=2S,,'0.pdp/(1—p,°) (59) 


is computed, it turns out to be independent of the 
exponent m, namely equal to: 


&.=pp:’/(1— 1’). (60) 


A tube with closed ends subjected to internal 
pressure p has the same mean axial stress 
6.= pr;’/ (ro —ry) = ppy’/(1—p:*). We must con- 
clude therefore that our condition u,=0 imposed 
at the beginning on the stress distribution ex- 
pressed by the Eqs. (53) is identical with the 
other condition that the mean axial stress ¢, in 
the tube is so large that it just balances the axial 
tensions which would be produced when closing 
both ends of the tube. It can indeed easily be 
shown that this is not only true for the resultant 
of the axial stresses o, but that the condition 
u,=0 is satisfied in any tube with closed ends 
subjected to internal pressure p. Indeed when 


o,=(d/dr)(ro,) and when o,=(o¢,+0;,)/2 


the integral 


forar=3 fordr+3 f d(ro.) dr \-rdr. 


Using partial integration and r=r,; and r=fr2 as 
limits this is found equal to 


J erdr= pri 2. 


The mean stress ¢, is therefore always pr,’ 
(ro2—r,*). This is the same value as obtained ina 
closed tube. In such a tube therefore the condition 
u.=0 prevails. According to what was said 
before, we note finally that the same distribution 
of stresses in a cylinder as shown in Fig. 6 and as 
found in the steady state of creep using as speed 
laws the power function law, Eq. (39), Fig. 2 
and expressed by (53) would have been derived 
in a tube of a plastic material flowing under 
pressure the strain-hardening curve of which for 
cold working would be expressed by a power 
function 
r=T*(y/y*)", O<m<il 


of the shears y or a set of curves similar to those 
shown in Fig. 7. 
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Stresses around cylindrical cavity 


Let 7; be the radius of the cavity and let an 
equal stress o act in larger distance from it in the 
directions perpendicular to the axis of the cavity. 
Let u,=0. The boundary conditions to be 
satisfied by the radial stress are r=r, ¢,-=0 and 
r= x, ¢,=o0=const. The integration constant c 
in Eq. (45) is found from the first of these con- 


ditions. When r=r7, let o-,=01, n=; then 
c=o,r;"/4n. Eqs. (46)—(51) give 
o,—0,=a1p 2m n= 7p?" m) (61) 


where p stands now as abbreviation for the ratio 
p=7\/f. 


The radial and tangential stress are now given by 
the formulae : 


o-=o[1—p "], | 


(62) 
o,=o[1+(2m—1)p"]. | 
At the surface of the cavity r=ri, p=1 the 
tangential stress is 
o,=2m<e. (63) 
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Fic. 8. Steady stresses around cylindrical cavity 
for various speed laws. No axial rate of flow: Uz=0 
stress concentration factor K =2m. 
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The 


usually is the ratio: 


stress concentration factor k defined as 


k=0o,/0=2m. (64) 


k is thus found to be equal to twice the exponent m 
appearing in the creep speed law: 


T=7T (0/0,)™. 


For a few values of m the steady creep stresses 
according to (62) are plotted in Fig. 8 around the 
cylindrical cavity. For m was taken 

—_e 2 1 1 1 
m=1, 3, 2) 3 6 


The value of m=1 corresponds again to a 
perfectly viscous substance, the distribution for 
o, and o, when m=1 is identical with the one 
around a cavity in an elastic material. The 
factor k=2. 

An interesting case occurs when m is taken 


1 


equal to 3. The tangential stresses are redis- 


tributed in a uniform manner and 
(65) 


o,-=—a(1—p)/p, o.=o=Cconst. 


When the exponent m in the creep speed law is 
smaller than 3 the stress o; at the cavity becomes 
smaller than o2.° 


5 In order to check the formulae (62) further the octa- 
hedral shearing stress 7; was computed at the bore. 1; is 
found from (37) taking ¢,=0, ¢,=0,;=2me which gives 


71 =01/\/6=(3)!- mo. (66) 


When m tends to approach the value zero, the octahedral 
shearing stress according to the definition of a perfectly 
plastic substance Eq. (8) must become equal to v2-00/3; 
oo is the yield stress in pure tension. Making 7; = V2-00/3 
for o the value 


o=a9/\/3-m (67a) 


and after inserting this in the expression for a, (Eq. 62) 
o,=a9(1—p-")/V/3-m (67b) 


is obtained. If the exponent m tends now to become zero, on 


account of 
1 —— 2m 0 s 
lim - ) =-=2Inp 
mae m 0 


the formula (67b) for the radial stress o, converges into 
the function: 


o,=(2e0//3) In p. (68) 
Similarly when m=0 o,; converges into the function: 
o1=(200/+/3)(In p+1). (69) 


These are the known formulae for the stresses around a 
cylindrical cavity pierced through a_ perfectly plastic 
substance. As can be seen from (67a) the stress o in this 
last case (m=) increases beyond all limits, so that cases 
when m tends to approach the value zero are of no practical 
interest. 
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Disk HAVING CIRCULAR HOLE AND STRETCHED 
IN Its PLANE. STRESS CONCENTRATION 
FACTOR IN CREEP 


Let the radius of the hole be 7; and the stress 
for large values of r be equal to o,=c. It is 
required to find o, and o; as functions of r and 
the value of the tangential stress o; at the bore 
after a redistribution of the stresses in the disk 
took place and the steady state of creep is 
reached. As usually in such cases a factor of 
stress concentration k=o,/o can be defined and 
it is particularly required to find the value of k 
for a given exponent m of the creep speed law: 

7T=7,(v/v,)™. (70) 
The normal stress in the perpendicular direction 
to the plane of the disk is zero: o,=0. For the 
radial and tangential stress and strain rate the 
following equations are available: 


u,=(20,—0,) /6n, (71) 
u,=(20,—<0,)/6n. (72) 
The compatibility condition 
d(ru,)/dr=u,, (73) 
the equilibrium condition 
d(ra,)/dr=a;. (74) 


The octahedral shearing stress + defined by 
Eq. (30), on account of ¢,=0 is given by the 
expression 


r=}-[2(0?—0,0,+<¢,’) }! (75) 
and the shear rate v by 

v=2-[F(u?+uu,+u,) }}. (76) 
From these 7 equations by gradual elimination 
two simultaneous equations for the two unknown 
stresses o, and o, can be obtained. One more 
relation already implicitly contained in the 

stress-strain rate relations, namely 
T= nv (//) 
allows the elimination of v. From (70) and (77): 
n=m(ti/7T)", (n=1/m-—1). (78) 
General integration of the preceding group of 
equations becomes possible after making a slight 
change in the expression of 7 in Eq. (75). When 
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r=r,, o-=0 and when r=, o,=0,=0. We 
observe that at both ends of the interval 
ry<r<@ the right side of Eq. (75) becomes the 
same expression, so that there 


r=Vv2-<a; 3. (79) 


This suggests to make 7+ in the whole disk equal 
to this last expression instead of using for r+ 
Eq. (75).® 


Substituting 7 from (79) into (78) gives 
n=171(01/0;)". (80) 


ao, and 7; are here the values of o; and 7 at the 
bore r=r,. Eliminating the strain rates from (71), 
(72) and (73) gives 


d {2e.-—<, 3 
r ( )+ (o,—0,)=), (81) 
dr n n 


and after substituting in this equation the value 
for » from (80) the first of the following two 
differential equations for o, and o, is obtained 
while the second is Eq. (74) 


r(d/dr)[o,"(20,—¢,) |+3e,"(o.—0,)=0, (82) 
r(do,/dr)+o,—0¢0,=0. (83) 

Using instead of r the ratio 
p=r/r; (84) 


and after changing the variables o,o, and p 
according to: 


x=20,—0,, y=o.—0,, 2=Inp. (85) 


Eqs. (82) and (83) are transformed in the two 
differential equations for x and y in function of 
the new independent variable z: 


nx(x’ —y’)+(x’+3y)(x—y) =0, (86) 
x’ —2y’—y=0. (87) 


The radial and the tangential stress o, and o; are 
given by 
o,=x—2y, o.=x-Y. (88) 


6 Examples show that in a thin disk stretched in its 
plane by tensile stresses the values of + computed in both 
ways (Eq. (75) or (79)) do not differ much. For a perfectly 
plastic material the quantity r= V20 /3 is const. (a9 yield 
stress in tension) and (75) represents an ellipse in the 
coordinates a,, ¢;. The points corresponding to the states of 
stress which are possible in the disk under the given 
assumptions are found in a comparatively short arc of the 
ellipse. According to Eq. (79) this arc is replaced by a 
secant situated near to the ellipse. 
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A set of particular integrals of (86) is easily 
found assuming that 


(89) 


y=ax, 


where a does not contain z. Substituting (89) in 
(86) and replacing 1+ by 1/m 


(1—a)x[x’+3max ]=0, (90) 


(1—2a)x’—ax=0 (91) 


is obtained. Eq. (90) is satisfied when either a=1 
or when x=0 or when the bracket is equal to 
zero. When a=1, 


and x=y=ce*=ce"" *’=c/p 


is a set of particular integrals. The corresponding 
stresses according to (88) are 


G,-=—C/f, 2. 


This is the radial distribution of stress in a sector 
which has been cut out of the disk. The sector is 
pulled in the radial directions only, the straight 
sides being free of stress. The second case (when 
x=y=0) is trivial. In the third case however 
when the bracket in (90) vanishes, we have the 
two homogeneous linear equations for x’ and x: 
x’ +3amx=0, 


(1—2a)x’ —ax=0. 


With 


vanishes: 


condition that their determinant 
a(1i+3m—6ma)=0. When 
g,=0,=C, which is again a trivial case. However, 


the 


a=(, 


when the second factor vanishes, we have 
a=(1+3m)/6m and x’ = —3max so that 
x= y a=ce 3maz =C grne (92) 


is a solution. In this case the boundary condition 


that for p=1, o,-=—p can be satisfied. The 

corresponding expressions for the stresses : 
(3m—1) 

o,=—p-p (L+3m)/2. o,.=———_— pp (1+3m) /2 (93) 


2 


show that this is the distribution of stresses in a 
large disk having a hole with the radius 7;. The 
inside of the hole is subjected to the pressure p 
in radial direction.’ 

7 On account of the practical impossibility to apply large 
pressures in the plane of a thin disk without producing an 
instability and buckling, this case need not to be dealt with 
further. Taking m= 1 the solution is identical with one in an 
elastic or a viscous material 
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Returning to the integration of (86) and (87) 





we eliminate x’ from these two equations, 
obtaining 
2y(y'+2y) 
x=— (94) 
(n+2)y’+(n+4)y 
From (87) we see that also 
x’ =dx/dz=(dx/dy)-y' =2y'+y. (95) 
After substituting now for 
y=y/w. (96) 


Eqs. (94) and (95) can be rewritten in the form: 





2y(1+2w) 
x= (97) 
n+2+(n+4)w 
and dx/dy=2+w. (98) 


If x from (97) is now inserted in (98) a differential 
equation connecting the variables y and w: 


d (1+2w)y w 
pal PEE 2 = poe 
dylLn+2+(n+4)w 2 





(99) 


is obtained. Eq. (99) can further be simplified by 
making use of the new variables w and v: 


Or y 


- 


w= 1 -—-—_ =——_ 
oO; x—-y 


and v=Iny. (100) 


From (100) 
1+w 


x= y. 
Ww 


(101) 


w is expressed by means of w after equating (101) 
and (97): 


n+2+nw 
—-— ——- (102) 
n+4+ nw - 


and this enables us to rewrite finally (99) using 
w and v as variables: 

dw n+4+(n—2)w 
a, (103) 
dv n+4+nw 
Integration of (103) gives finally : 


v=In y=Inw 


+ (2/n—2) In [n+4+(n—2)w ]+const. 
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Fic. 9. Stresses around hole in disk in 
viscous (m=1), in creeping (1 >m>0 
and in plastic (m=0) substance. Stress 
concentration tactor 


K =(2/(14+3m) )?™/1-3m 


This equation, when u is written as abbreviation 
for the bracket can be rewritten as 


y=cwu?!* (104) 


and we see at the same time from (101) that 


x=c(1+w)u?!("-?), (105) 
The formulae (104) and (105) express the distri- 
bution of the stresses in the disk as the stresses 
o, and a, linearly depend on the variables x and y. 
Using (88) we obtain the stress expressed in 
parametric form: 

o,=cC(1—w)u?™, | 
‘ (106) 

a, =cu 
where u=n+4+(n—2)w. It remains finally to 
express the radius r or the variable p=r_r; 
as a function of the parameter w. From (96) 
dz/dy=w_ y. Remembering that v=In y, dv=dy/y 
we see that 

dz dv=w (107) 

and that 


dz dw=(dz/dv)-(dv/dw)=w-dv/dw. (108) 


Since the right-hand terms w and dv, dw in (108) 
are known functions of w given by Eqs. (102) and 
(103), respectively, (108) is the differential equa- 
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tion for z=In p as a function of the parameter w. 
We obtain 
dz n+2+nw 


=— ; (109) 
dw wl n+4+(n—2)w | 


Since the boundary conditions for a disk having a 
hole r=r, and stretched by a uniform stress o are 
p=1, o-=0 and p=, o,=o0,=0, the corre- 
sponding values of the parameter w=1—<a,/ 0, 
must be: when p=1, w=1 and when p=“, w=0. 
Eq. (109) is therefore to be integrated under the 
condition: w=1, z=In p=In1=0. The integral 
of (109) satisfying these latter conditions is: 


r 


9 


u 4(n+) [(m—2)(n+4)] 
x| | . (110) 
2(n+1) 


The integration constant cin (106) is found from 
the second boundary condition that when w=0, 
a,=a. This is the case when c=(n+4)o?/"™, 
The two Eqs. (106) together with (110) solve the 
problem completely. 

The steady creep distribution of the stresses 
around the hole in the disk is given in parametric 
form by the following three equations : 


n—2 \2/(n-2 
rad. stress = (1—w)o( 14 ) . £eEn) 
n+4 


n—2 ora 
tang. stress imo 1 + «) : (112) 
n-+4 


nt4+(n—2)w yr" n2—4) 
| . (113) 
2(n+1) 


The factor of stress concentration k is found 
from (112) by taking for w=1 


k=a,/0=([2(n+1)/(n+4) P/® (114) 


or after introducing here the exponent m ap- 
pearing in the creep speed law, this factor is 
found to be equal to 


2 2m/(1—3m 
b=( ) , (115) 
1+ 3m 
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Exponent m=1| 2/3 1/2 
n=(1/m)—1=0} 1/2 1 
Factor of stress concentration k=2 | (2/3)~*3 (4/5)-? 
=2)| 1.72 1.56 


Substance: elastic or viscous 


Formula (115) is validfor0<m<1.Asn=1/m—1 
Eqs. (111-114) are valid for 0<n<a. The 
preceding equations must contain as special 
limiting cases the distributions of stress in a disk 
around a hole for a perfectly viscous and for a 
perfectly plastic substance. The first limiting case 
is obtained taking for the exponent m the value 
one. When m=1, n=0 we obtain from (113) 


and from (111) and (112) 
2(1—w)e, 20 
F,= ’ ¢,;=-—~ ’ 
2-—w 2—w 
after eliminating w we obtain 


o,-=a(1—r;/r*), o,=o(1+r;?/r*) (116) 


or the well-known formula for the stresses in a 
disk having a hole, valid in a viscous as well as in 
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Fic. 10. Factor of stress concentra- 
tion K = (2/(1+3m))?"/@-38™ plotted 
versus exponent m of speed law 

=17,(v/v,)™ for disk with hole. 
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TABLE I. Factor of stress concentration for creep and plastic substance. 


1/3 | 1/4 1/5 1/6 0 
2 3 4 5 % 
el/3 (8/7)? 5/4 (4/3)2/3 | 1 
140 | 1.31 ae a 
after steady creep | plastic 


an elastic substance. When m=1 Eq. (115) gives 
the factor of stress concentration k=2 known 
for the elastic substance. 

The second limiting case is given when 
m=0, n=. From the Eqs. (111-115) one 
finds w=7r;/r, o,=o0=const. and o,=(1—w)e 
=o(1—r,/r). These are the well-known formulae 
for the stresses in a disk having a hole for a 
perfectly plastic substance assuming that the 
disk material is yielding uniformly under a 
constant maximum shearing stress r=0/2. The 
factor k is in this case equal one. 

Table I contains the values of the factor of 
stress concentration k computed with the formula 
(115) for a few values of the exponent m of the 
creep speed law. 

The values of & are also plotted in Fig. 10 in 
function of the exponent m. The curve repre- 
senting the values of k=f(m) is not different 
much from the straight line connecting the points 
m=0, k=1 and m=1, k=2. A stress concen- 
tration factor in the steady state of creep can 
therefore be taken approximately as 


k=1+m (117) 


under the assumptions as to the creep speed law. 
A few distributions of the stresses are shown in 
Fig. 9. One interesting case deserves perhaps to 
be mentioned especially. This latter occurs when 
the exponent approaches the value m = 3 or n the 
value 2. The formula for k in this case becomes 
indefinite and can be written as 


k=lim [1+2(m+4)/2(n—2) PO 
or after writing here for 2/(m—2) =a also as 


1 a 
k=tim (1+) =e!=1.396. (118) 


a= 3a 
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The formulae (111-113) fail in this case, direct 
computation however furnishes when m=2 the 
rather simple functions 


a(1—w)e*’’, r; 
=wi-e@-D/3, (119) 
o,=ae*'*, r 


a, 


again valid for 0<w<1. Taking w=1 the tan- 
gential stress at the bore is obtained equal to 


o,=e'!-¢=1.3960 (120) 


in agreement with the value of k just computed 
above in (118). 

We note again that the formulae expressing 
the stresses o, and o, Eqs. (112), (113), (114) and 
the factor of stress concentration k Eq. (115) 


or (117) are also valid in the case of the strain 
dependent plastic flow in a thin disk having a 
hole, the material of which is strain-hardening 
according to the power function law 


r=7*(y/y*)™. 





Notes on the Wetting of Filaments by Molten Metals 


M. ALDEN COUNTRYMAN 
Physics Department, Iowa State College, Ames, Iowa 


(Received April 2, 1937) 


HE following results of some experiments on 

the wetting of hot filaments by molten 

metals are presented in the hope that they may 

prove helpful to other experimenters who, like 

the author, may have searched the literature in 
vain for such information. 

The object of the experiments was to deter- 
mine which filament material was most desirable 
for use in evaporating a metal directly from the 
filament. In general if the molten metal adhered 
to the filament it would evaporate successfully 
unless its vapor pressure was too low at the 
highest safe temperature for the filament or 
unless it united with the filament to form a 
low melting point alloy. 

The apparatus used consisted of a bell jar 
placed on a flat plate of heavy glass ground to 
fit; the joint sealed with mercury. The lead-in 
wires were introduced through a seal in the 
evacuating tube which in turn was sealed into 
the lower side of the flat glass plate by means of 
a ground joint sealed with de Khotinsky cement. 
Included in the system near the bell jar side of 
the trap was a discharge tube having two alumi- 
num disk electrodes 24 mm in diameter spaced 
12.8 cm apart. By applying the potential from a 
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hand Tesla coil to one electrode, the other being 
grounded, and observing the form of the dis- 
charge it was possible to make a rough determi- 
nation of the pressure in the system. 

The metal to be evaporated was secured in 
the form of a wire of convenient size and wound 
around the filament to be tested. The wound 
filament was then formed into a spiral with 
short straight ends. These were inserted into 
connectors provided for them; pieces of glass as 
test strips were then placed above and below 
the filament at a distance of about 5 cm and 
the bell-jar lowered into place. When the po- 
tential from the Tesla coil failed to produce a 
discharge across the test gap the pressure was 
considered low enough for satisfactory evapo- 
ration. 

A potential was applied to the filament and 
the current gradually increased until the metal 
melted and either adhered or dropped off. If it 
adhered, the current was further increased until 
the metal evaporated (as indicated by the 
coating of the strips of glass included in the 
system) or the filament burned out. 

In Table I three results are indicated: (1) 
whether or not the metal when molten ad- 
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Aluminum 
658.7°—1800° 


Silver 
960.5 


(x ld 
1063°—2600 


Copper 
1083°—2310 


Ni kel 
1452°-2900 


Chromium 
1615 


Platinum 
7 


1755°-4300 


Aluminum 
658.7°- 1800 


Silver 
960.5°—1950 

Gold 
1963°- 2600 


Copper 


1083°—2310 
Nickel 
1452°- 2900 


Chromium 


1615°-2200 
Platinum 
1755°-4300 


hered to the filament, (2) whether or not 
metal evaporated successfully, and (3) whether 


or not 


1950° 


2200° 


Cc 


Cc 


~ 


NICKEL 
1452°-2900°C 


No experimental evidence. 


| 1. Wets perfectly. 


TABLE I. 





PLATINUM 
1755°-4300°C 


2. No success—filament burned out in 
spite of great care. 


3. Metal probably united with filament 


1. Wets perfectly. 

2. Evaporates well. 

3. Probably alloys—deposited 
may contain some nickel. 


metal 


No experimental evidence. 


1. Wets perfectly. 

2. No success. 

3. Copper probably alloys with nickel to 
form low melting point alloy. 


Requires care but evaporates very suc- 
cessfully. 


Obviously impossible. 


Obviously impossible. 


TANTALUM 
2850°-5300°C 


1. Wets perfectly. 
2. Evaporates very well. 
3. 


No obvious effect of metal on filament. 


Melts and drops off. 


1. Melts and clings in globules. 
2. No success 


Melts and clings in droplets. 
terial. 


1. Wets perfectly. 
2. No success—filament burns out 
3 


. Nickel seems to unite with filament to form low melting 


point alloy. 


No experimental evidence 
(Metal not available in wire form.) 


. Melted, then sputtered and seemed to burn up. 


1 
2. No success—filament burned out 
3 


. Probably due to interaction between metal and filament. 


there was obvious interaction between 
the molten metal and the filament material. 
The filament material used is given at the 
head of each column, and the metals which were 
to be evaporated are listed in the column at the 
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gold appears coated with a slag or scum. 
3. Gold may have united with filament to form alloy. 


to form low melting point alloy. 


No experimental evidence. 


No experimental evidence. 


. Evaporates fairly well. 
No apparent interaction 
metal and filament. 


1. Wets perfectly. 
? 
3. 


between 


No experimental evidence. 


No experimental evidence. 


No success as yet. 


. Wets perfectly. 


MOLYBDENUM 
2620°—3700°C 
| 1. Wets perfectly. 

2. Evaporates but requires careful con- 
trol of temperature. 

| 3. Not much apparent alloying. 





| Melts and drops off. 


1 
| 2. Evaporates well. 
3. No obvious effect of metal on filament 
1. Melts and clings in droplets. 
2. Evaporates well—though drops occa- 
sionally fall off. 
3. No obvious interaction between copper 
and filament material. 
1. Wets perfectly. 
2. No success—filament burns out. 
| 3. Nickel seems to unite with filament to 
form low melting point alloy. 


No experimental evidence. 

(Metal not available in wire form.) 

1. Wets perfectly. 

2. Not much success—vapor pressure of 
metal too low at safe maximum fila- 
ment temperature. 

3. No obvious interaction between metal 
and filament. 


TUNGSTEN 
3370°—5300°C 


1 
2. Evaporates but requires some care. 
3 


. Filament appear 


s to have been affected slightly by metal. 


| Melts and drops off. 


1. Wets but hangs 
2. Evaporates well. 
3. No apparent inte 


in gk ybules. 


-‘raction between metal and filament. 


Melts and drops off. 


Evaporates fairly well, though drops occasionally drop off. | 
No obvious interaction between copper and filament ma- 


1. Wets perfectly. 
2. Not successful 
3 


Tried chromium plz 


Wets perfectly. 


wr 


the 
melting point and 


The writer ackn 


department who su 


filament burns out easily. 


. Probably some interaction between metal and filament. 


ating filament with little success. 


Requires care but is fairly successful. 
No obvious interaction. 


left. Under the name of each metal is given the 


the boiling point in degrees 


centigrade at atmospheric pressure. 


owledges with gratitude his 


indebtedness to Dr. P. H. Carr of the physics 


ggested the research problem 


and did much to encourage its continuation. 
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A Note on a Problem in Potential Theory 


Morris MusKAT 
Gulf Research Laboratory, Pittsburgh, Pennsylvania 


(Received March 22, 1937) 


Two methods are given for treating a new type of poten- 
tial problem arising in a recent study of the author on the 
encroachment of water into an oil sand. The problem is that 
of finding the potential distributions in two regions of 


different ‘‘constants”’ (“‘conductivities”), separated by a 
moveable surface of unknown shape, this motion and shape 
to be determined by the conditions that each point of the 
interface should have a velocity proportional to the vector 
gradient of the potential at the point, and that the area 
swept out by the moving surface shall assume the ‘‘con- 
stant" appropriate to that of the encroaching side of the 
interface. The first method is a perturbation method, valid 
when the difference between the ‘‘constants”’ is small, and 
reduces the problem of determining the shape and motion 
of the interface and the potential distributions simul- 
taneously, to a series of problems in which potential dis- 
tributions are to be found on the sides of known interfaces, 
and the motion of the surface is to be found in a known 
(time varying) potential field. This same type of reduction 


N a recent paper in this journal,' the author, 

in discussing the encroachment of water into 
an oil sand, outlined a new type of potential 
problem. This problem was formulated as fol- 
lows: Find the potential functions , between a 
fixed surface S, and the moveable surface 
F(x, y,2,t)=0, and #2 between F(x, y, 2, t)=0 
and the fixed surface S,, such that? 


>,=V,,(x, y, 2, 4) on Sy, 


b,=V,(x, y, 2, t) on S,, 
?, = Po, | 
on F(x, y, 2, t)=0, 


4 





C,0PD,/OnN=C2d0P2/On 
where 
OF /at Cs. VP), eVF= 0; Fi X,Y; 2, Q) = Sie 


n being the normal to F. The system is diagram- 


1M. Muskat, Physics 5, 250 (1934). 

* The time is to be considered as a parameter when it 
occurs in the potential distributions, as the changes in the 
potential distributions due to changes in the boundary 
values or positions of the interface are to be considered 
as taking place instantaneously. 
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is effected also in the second method, which is one of direct 
successive approximations, but which may be used for any 


values of the potential ‘‘constants.’”” The perturbation 


theory is carried through explicitly to the third order for 
the case of a one dimensional system, and agrees to that 
order with the exact results which were derived for this 
case in a previous paper. The method of successive approxi- 
mations is carried out for this same case through the second 
approximation for a ratio of the “conductivity” of the 
interior side of the interface to that on the exterior side 
equal to 5, and through the third approximation for the 
ratio 1/5. In the former, the calculations give a rapid 
monotonic convergence to the correct solution, and in the 
latter the convergence is oscillatory, and only slightly less 
rapid. These behaviors are shown to be generally true for 
the linear problem, in the course of the construction of an 
analytic convergence proof for the one dimensional problem. 
The convergence proof of the method of successive ap- 
proximations is also given for the case of radial symmetry. 


matically indicated in Fig. 1, where the moveable 
surface F is denoted by Sj, its initial form being 
taken for convenience as S,. 

This problem was new in the sense that 
instead of the normal potential problem in 
which the potentials are to be found on the two 
sides of a fixed interface separating two regions 
of different “‘constants,”’ that in question involves 
a moving interface, the motion and shape of 
which must be determined simultaneously with 
the potential distributions. It corresponds physi- 
cally, in the problem where it arose, to a system 
in which oil, of viscosity uw, originally occupying 
the sand between S, and S,, is displaced by the 
water, of viscosity ue, encroaching into the oil 
sand across the surface S,, the pressures at S, 
and S, being considered as known. The surface 
S; is then the continuously changing water-oil 
interface, and the constants c are the reciprocals 
of the liquid viscosity. The potential functions 
® are the velocity potentials kp of the liquid 
(neglecting effects of gravity), where k is the 
permeability of the sand, taken as uniform, and 
p is the fluid pressure. 
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No general treatment of this problem was 
given in the above mentioned paper. Rather, 
rigorous solutions were presented only for the 
three simple cases of linear, radial and spherical 
flow, where the shape of the interface S; was 
obvious a priori from the geometry of the 
systems. For the more general cases, the discus- 
sion was based upon the “zeroth” approximation 
to the correct solution, in which the differences 
between c, and cz were neglected, and the shape 
and motion of S; were taken to be those of a line 
or surface of fluid particles in a homogeneous 
system of the same geometry, leaving S, at the 
same instant as the real passage of the water 
across S, was supposed to take place (cf. Eq. 
(4) below). 

In this note two methods will be presented for 
reducing the problem of the unknown moving 
interface to a series of potential problems, to be 
solved successively, in each of which the shape 
and motion of the interface are to be considered 
as known from a previous approximation. The 
time will enter in the analysis of the potential 
problems only as a parameter defining the shape 
and position of the interface. The first method 
will consist of a formal perturbation theory, 
applicable to cases where the difference between 
the constants ¢,, C2 may be taken as small. The 
second will be a method of direct successive 
approximations valid for any values of c; and ¢¢. 
Although no general proofs have been developed 
regarding the convergence of the methods, the 
correctness of the perturbation theory to the 
third order and the rapid convergence of the 
successive approximation method in the simple 
case of linear flow, where the exact solution was 
derived in the paper cited above, will be demon- 
strated by the actual calculations. 


PERTURBATION THEORY FOR |¢2—¢1, = || Ki, 2 


Taking for convenience ¢;= 1, and c.=c=1+n, 
the perturbation theory consists simply in con- 
sidering the functions #,, #2, and F to be ex- 
panded in power series in the small parameter 7. 
Thus: 


P, =Pjo+Pii1+77Pi2+ ---,) 
Do = Doo +P21+° Poet * ++, (2) 
F= Fo+nF\+n?Fe+ Pe 
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where all the %’s are potential functions. Ap- 
plying these to Eq. (1) and denoting by 4, 2°”, 
F™ the resultant values of ®,,2, F to the mth 
approximation, the general problem in the mth 
approximation may be expressed as that of 
finding the solutions to the equations: 


V76, =y7°6,( =0, ) 


p,(™ =W,,on S, - @,\=V, on Sn 


P,™ =H, 
| 
OP, (m) OP, ' m) IP» (m) OP» | m—1) > , 3) 
—— =(1+n)—— = +n : 
on on on on 


on F("—)(x, y, 2, t) =0, 


a F' m) 





: =7?,™ -yF™. FO) (x, y, z, 0) = 5, 
Ot 


or OF,,/dt= > VHi;-V Fi. 


j+k=m 


Carrying this formulation through explicitly 
for the first two approximations, we have: 


Zeroth approximation : 
V*bi0=V"*P20= 0, 


Vv, on Sy, 


Y 


Pio = Popo as 
V,on S,, 


OF)/dt=VPi0:'V Fo; Folx, y, 2, 0) = a 





Fic. 1. 
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C2(u2)=C Ci(y)= 1 -$=Y(t)=0 





/ PeeogP 622 228020 ceeee8 aZ--- a ane 4. 


}=,(t)=! 














Fic. 2. Diagrammatic representation of a linear 
encroachment system. 


First approximation : 
V°?1 1 =V "bo = 0, 
©:;,=00n S,; ®2,=00n S,, 
Pi, = Po) | 
OPD;; On = IAPDe;/ ON+ APa9 an) 
OF, dt=TP io VPitVeuV Fo, 


F‘) = Fy+nF,—S,, at t=0. 








Second approximation : 
V*Pi2=V Poe = 0, 


P jo =(Q on ‘a : Po. = 0 on a 


Pio Poo + (1 n) (Po, — P14) 
OPi2 IPo2 = IPe; | 
= } | 

on on on fon FO, 

(6) 

1/dOD2o IOD2,; IP), 
} +——--— }| 
n\on on an J } 


OF 2/0t=TPi0° TV Fot VP VPs t+VPia VV Fo, 
FO = Fo+nFi+n°F.—-S,, at t=0, 


Ste. 





THe LINEAR PROBLEM 


For the one dimensional case, diagrammati- 
cally represented by Fig. 2, where the interface 
in the whole system, of unit length, is defined 
by its distance x» from the inflow face S,, with 
a unit potential difference between ,S, and S,, 
i.e.: V,~=0, V,=1, the above equations lead to 
the following analysis: 
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on Fo, (5) 


Zeroth approximation : 


Pio = Poo — 1 um Ze 


OF,/dt= —OF)/dx; Fo=t—xo=0. 


First approximation: 


);=a1;(x—1); Po; = 21%, 


where 


41(X9 — 1) =Ae1X0; Q\;=a2,—1, 


so that a),;= —X9=—l; deo =1—xX)=1-1. 


Further 


0 F, ot= —O F, ox—a S2> F, =t—x-+t i Z. | 
Hence | 
(t—x™)(1+7)+nt?/2=0; 





x) =t+nl?/2(14+ 7). J 
Second approximation : 
Pio=Ay2(xX—1); Po2=do22x, 
Ayo(x6? —1) =deox’? +(1/n) | adoix"! 
— ay (x6? —1)} =deox? +£/2(14+7n), 
19 =d22+de,4+ (1 n)} —1+da1—-d11; 
=doet+1-t. 
Hence 
t* 
dy2=t(1-—t)— ji—n(1—2)}, 
2(1+ 7) 
t- 


ao22>= —(1-—t)*-— 
2(1+7) 


\1—n(1—2)}. 


Further, 
0 F, ot= —d F, Ox—a 11 @12, 
so that 
{3 £3(1 —n) nt' 
F.= t—x+-+ + . 
3 6(1+n) 8(1+7) 


Finally 
nt? nt? 


(¢—x)(1++7n*)+ : + + 


n’?(1—n)t 


6(1+ 7) (9) 


8(1+ 7) 
and 





x?) =t+ nl? /2—n7t?/2+77t?/2+0(n*). 
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Third approximation : 
Here the extensions of Eq. (6), reduced to the 
problem of linear flow, are: 


Di3=Q13(X—1); Po3=Ao3x, 
P13 = Po3+ (1/7) (P22— Piz) 


+(1 n*) (Po; — P41) on x), 


OP}; OPo3 IPoe 1 IOP oo OPo, OP jo 
=o pt $a ) 
Ox Ox ox n\X Ox Ox Ox (10) 
1 | OPo; IPo OP), | 
+—| —+—-—} 


n* | Ox Ox Ox | 
OF; OF; 
=— —@11—@12—4}3, 
ot Ox 





where, as F; will enter into x“ with the coeffi- 
cient 7°, only the zero-order terms need be 
retained in dj. and d,3. Solving these in the same 
manner as those for the second approximation, 
it is found that to the zero order 


Q43= —t+2f—28, 
so that 


F,;=t—x+/2—/6+14/2, 
and finally 


2 


am 242 wt? 
x= 
2 


“ 9 
+ (t—1)+ 
2 8 


(4—8t+52*). (11) 


That these approximations of Eqs. (8), (9), 
and (11) are correct to the orders indicated, 
may be readily verified by expanding the exact 
solution for the one dimensional problem, 


namely : 
a+ly 2nt \*| 
Xo= 1 —_ 1 —_ ) ° (12)3 
n | nt+1 


Although the actual determination of the shape 
of the interface F has been unnecessary in the 
above perturbation calculations, just as in the 
case of the rigorous analysis of the linear prob- 
lem, the analytical content of the perturbation 
method is fully illustrated by the treatment of 
the one dimensional case. The only essential 
difference between the one and the two and 
three dimensional problems is that in the latter 


3 Reference 1, Eq. (11), where «€ is to be replaced by 
1/(1+ 7) and xo by 1—x. 





VOLUME 8, JUNE, 1937 


the partial differential equations for F,,, the last 
of Eqs. (3), will involve two or three of the 
coordinates (x,y,z) rather than just x. In all 
cases, however, the program contained in Eq. 
(3) gives an explicit reduction of the original 
problem of determining the shape and motion 
of the interface and the potential distributions 
simultaneously, to the series of problems in 
which these parts of the complete solution may 
be derived by independent analytic processes. 


THE METHOD OF SUCCESSIVE APPROXIMATIONS 
The method of successive approximations 
proposed here begins with the same zero approxi- 
mation as the above perturbation theory. It 
may be formulated as follows, where the func- 
tions ®,,, Pea, and F,, represent here the complete 
solutions to the mth approximation. 
V7Pin = V Pom == (), ) 
d,,=¥,0n S,; 9,2, On Se, 
| on | (13) 


? es Pom, 


COP ym, on= C2dOPom /on 





OF, ot — Dm “UF. Fault, y, Z, 0) _ Ses j 


where 3,, is the velocity vector —C,V®im at Fr. 

Thus the first approximation involves the 
finding of the strict potential distribution for a 
system with constants ¢, (2, separated by the 
known moving interface Fy given by the zeroth 
approximation. The motion of a line of fluid 
particles in a homogeneous medium of constant 
¢,, but with the velocity field corresponding to the 
new first approximation potential distribution, 
gives the new interface F,. This same process is 
repeated with F, taken as the known moving 
interface, to give the second approximations 
Pio, Poe, and then F2, etc. Although the physical 
interpretation of this approximation method 
would lead one to expect that the process will 
actually converge to the correct solution, an 
analytical convergence proof has not as yet been 
constructed except for the cases which had 
previously been solved rigorously. Nevertheless 
the application of the method to these simple 
cases does show that at least there the numerical 
convergence is quite satisfactory, as well as 
attainable in an analytical sense. 


437 





THE LINEAR PROBLEM 


Using the same notation as in Fig. 1, except 
that the constant co=c for region (2), which 
really represents the ratio ¢2/c, is no longer 1+, 
we have: 


Zeroth approximation : 


P19 =Px_p9=1-—x; Fo=t—x=0, 





(14) 
cr =f 
First approximation : 
Pi, =41,(xX—1); Oo,=—adax+1, 
(15) 
Ayi\xX ° —1)=do,x)+1; a\; > C€daQ\, 
so that 
—1 
on —- = ’ 
1+ at 
—_— i a +1) l—c 
ade; = ° a= 9 
1+at Cc (16) 
OF, —1 OF, of dt 
. , Fy=x- =0, 
al 1 + al Ox 7/9 1+ at 
x“? =(1/a) log (1+ eat). 
Second approximation : 
Pj\2=A;2(x—1); Poo =doaex+1, 
(17) 
Qyo(x"? —1) =daox'? +1; 6u™ Con, | 
so that 
d\2= —ve= —1/(1+log (1+ <a?)). 
Hence, 
0 F, —1| ra F, 
ot 1+log (1+at) dx 
which gives 
vf dt 
x= 
79 1+log (1+at) 
1 (18) 


II 


| ka(—u)— Fi(1)}, 
ae 
where 
—u=1+log (1+ at). 


Carrying through the same procedure so as to 
obtain the third approximation, it is found that 
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‘ dt 
xD) = [ on ; (19) 
Yo 1+ax 
and in fact, for the mth approximation, 
vf dt 
xV= | ; (20) 
J 1+ax 1) 


which may also be derived directly on noting 
that the velocity in a linear system separated 
at x“) into regions of ‘‘constants” 1 and ¢ is 
simply the integrand of Eq. (20). 

These formulas have been applied numerically 
to the second approximation for c=5, (a= —0.8) 
and the third approximation—this last by graphi- 
cally integrating Eq. (19)—for c=2, (a=4). 
The results are shown in Fig. 3. It will be seen 
that for c=5 the second approximation already 
gives a very close approximation to the correct 
motion of the interface. In the case c=}, the 
convergence is not quite so rapid. However, in 
view of the large difference between the correct 
and zeroth approximation solutions, the con- 
vergence may still be considered to be quite 
satisfactory. 

The monotonic convergence for c=5(>1) and 
oscillatory convergence for c= 4( <1) are general 
consequences of Eq. (20). In fact, noting that 
the correct solution x(t) obeys Eq. (20) in the 
limit where x") =x =x, namely : 


t dt 
x)= ff . (21) 
0 1+ax 


it readily follows that for 
a<@Q, then x” (t)=x(t), 

as x®-)(t)Sx(), 
while for (22) 
a>d, x™(t)h=x(t), 

as x—)(¢)=x(?t). 

Observing, furthermore, that 
a<0, x™(t)Sx%-) (2), ) 
as) x—)(t)Sx"-?2)(2), | 

a>0, x(HjSx%—- (2), 


as x ®-)(t)=x-?)(2), 
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it follows that since for a<0, x(t) >x(t) 
>x(t), x(t) must increase monotonically 
and converge to x(t) (or at least to values less 
than x(t)). Similarly, since the inequalities of 
Eq. (23) can also be applied to the elements of 
the sequence x, x, x@ ‘4, it follows for 
a>0 that since x >x®: 


XO) Sx) D>z4M D> -- x(t) 
(24) 
KX xO CHO) << -+- x2), 


thus proving the convergence to x(t) of both 
the odd and even sequences of the approxi- 
mations. 


THE RADIAL PROBLEM (CF. FIG. 4) 


The equation for the interface in the mth 
approximation is here evidently 


OF ,,/dt= —v,0F,,/ 9dr, ( 


i) 
wn 


and from elementary potential theory it is 
readily found that: 


A® 


Z —_—_— — 


r{ log (re/%w) tea log (7/179 )] 


a=1/c-—1, (26) 


ro Y(t) giving the position of the interface to 
the (7—1)st approximation. 
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“> 14 16 16 20 22 24 26 26 30 
t 


Taking the function F, in the form 
F,=r—ro™ (t) =0, (27) 


Eq. (25) becomes 





Ory! mr _ A® 
ot ro [log (re/fw) ta log (r./ro"—) J 
so that 
ry\™ 2 (28) 
r. 
4A f dt 
r.2 Jo log r.7/fw? +a log (7,/ro"—" (t))? 


Furthermore, the exact solution 7o(/) satisfies 
the equation 


To : 4A ‘ dt 
(“) -1- f . (29)4 
r, r.2 Yo logr2/n,°+a log (r2/1r9") 


From Eqs. (28) and (29) it follows that for 


ew |} 


a<O:ro™(t)Sro(t), as ro"—?(t)Sro(t), 
while for (30) 


a>O: ro (t)Sro(t), as ro"? (H=ro(t). 


4M. Muskat, reference 1, Eq. (18). 


439 











Fic. 4. Diagrammatic representation of a radial 
encroachment system. 


And from Eq. (28), one may readily show that 
for: 


a<O : 79°? (t)Sro""— (2), 

| 

as) ro"—)(t)Sr9""—*) (8), | 

, | 
while for ~(31) 

( | 

a>O : ro (t)Sr_o—" (8), | 

| 

as 7ro\"—)(th=ro‘"—-?? (2). 


The convergence of the sequences of approxi- 
mations then follows by choosing® for a<0, 
zero approximation functions for which either 
ro) <ro™ if ro >ro(t) or ro > ro if ro <ro(t), 
and for a>0O, zero approximations for which 
if ro >ro(t), and rmO>ro if ro 
<r(t). For the zero approximation as given by 
Eq. (28) (a=0), it follows from Eqs. (28) and 
(29) that ro >ro(t) and ro <ro for a<0, and 
ro <ro(t) and ro >ro for a>O0, thus giving 


ro?) <r, 


* Insofar as the problem is reduced to the successive 
application of Eq. (28), one may choose for the zero 
approximation any function ro (t), the choice affecting 
only the rapidity of the convergence once the proper 
inequality with respect to ro” or ro is satisfied. 
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the conditions for the actual convergence of the 
resulting sequences. 

In an exactly similar manner one may demon- 
strate the convergence of the method given here 
in the problem possessing spherical symmetry, 
hydrodynamically to that of 
spherical encroachment. 


corresponding 


Although here, too, the method of successive 
approximations has been illustrated by examples 
in which the shapes of the interface have been 
known a priori, the determination of this shape 
will, in the general problem, involve no analytical 
processes other than those already occurring in 
the treatment of the linear and radial cases. 
The determination of the potential distributions 
in the various approximations will again reduce 
to potential problems with known surfaces of 
discontinuity, the time entering only as param- 
eters. The shape and motion of the interface 
will be given by the solution of the same partial 
differential equation, such as the last of Eq. (15), 
except that it will involve 3 or 4 of the inde- 
pendent variables (x, y, 2, ¢) rather than just the 
two (x,t). Indeed, the solutions of these indi- 
vidual problems will in general be very difficult 
and perhaps even be intractable, from a practical 
point of view,® but at least they will all be of a 
kind that may be considered as part of known 
analytic methods. 

The writer is indebted to Dr. L. J. Peters for 
helpful discussions of the problem treated here 
and to the executives of the Gulf Research & 
Development Company for permission to publish 
this paper. 

6 In fact, if one is willing to use graphical or numerical 
methods, it will probably be easier to attack the original 
problem directly by the following step-wise integration 
process. First, ®; and #2. would be computed (graphically 
or numerically) for the initial position of S;. The velocity 
vectors —cy® would then be constructed over S;, the 
terminae of which, with a proper scale, will give the position 
of S; at a short time interval after the initial instant. For 


this new position of S;, ®: and #2. would again be com- 
puted, the velocity vectors drawn on Sj, etc. 
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Heat Flow in an Infinite Solid Bounded Internally by a Cylinder 


Lioyp P. SMITH 
Cornell University, Ithaca, New York 
(Received April 1, 1937) 


The problem of the flow of heat through a long tunnel wall in an infinite solid has been 
treated and the results for the rate of flow and total heat flow have been expressed in a con- 
venient form for numerical computations. This general problem has become of interest in 


connection with the cooling of deep mines. 


N deep mining operations it has become in- 

creasingly important to cool the tunnels to a 
point where the temperature is such that men 
can work efficiently. In newly excavated tunnels 
the rock temperatures sometimes run from 120°F 
to 140°F and in order to reduce the temperatures 
so that work can be carried on efficiently in a 
reasonable time it is necessary to know the rate 
at which heat flows through the tunnel walls in 
the course of time. 

In order to ascertain the heat flow under actual 
conditions it is necessary to have the solution of 
the following heat flow problem: Given an 
infinite solid of uniform characteristics bounded 
internally by a very long cylinder of radius R. 
At time t=0, the temperature through the solid 
will be a function f(r) of the radial coordinate 
alone, i.e. w=f(r) for r>R and t=0. In practice 
deep borings have shown that the initial tempera- 
ture distribution is essentially constant. Since 
the wall temperature may vary in the course of 
time it will be considered as an arbitrary func- 
tion of the time for the present, i.e. w=u,(t) for 
r=R. A heat flow problem of this type has been 
considered by J. W. Nicholson! but the for- 
mulae are not very convenient for numerical 
computations. 

The mathematical problem reduces to that of 
finding a solution of the heat equation 


0u/dt=a(d*u/dx*+0°u/dy*?+0*u/dz"), (1) 


which satisfies the boundary and initial condi- 
tions specified above, where a is the diffusivity 
of the solid material. Since the temperature 
does not vary along the length of the cylinder, 
it may be shown? that a solution of Eq. (1) in 


! Proc. Roy. Soc. 100, 226 (1921). 
2 See Carslaw, The Conduction of Heat, p. 172. 
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cylindrical coordinates which satisfies the bound- 
ary conditions laid down above is given by 


eo) 


u(r, 0,1) = | ae’ | g(r, 0;7', 0: t)f(r')r'dr’ 


0 A 
at 2r 
+a arf u,(R, 6’, 7) 
0 “o 
| 98 , | ’ 
Xi—(7, 0;7°;t-—7)} Rdé’ (2) 
lar’ Deas 


provided the proper Green's function g(r, 6; 
r’, 0; t—r) can be found. The proper Green’s 
function is that solution of the heat equation 
which gives the temperature at any point 1’, 6’, 
at time ¢, due to a line source of unit strength 
located at the point 7, @ at time rt, the solid being 
initially at zero temperature and the surface 
being kept at zero temperature. Therefore, to 
obtain the Green’s function for this case, it is 
convenient to start with the solution of the 
heat equation representing the temperature from 
a line source of unit strength located at the 
point 7, 6, namely 


1 —r?+r’—2rr’' cos (0—6') 
exp ; (3) 
4rat 4at 


In order to satisfy the condition that Green's 
function must be zero at r’ = R it is convenient to 
expand the above function in the Fourier 


series 
1 9 
——exp —7r°+r’°—2rr’ cos (@—6’) 
4rat 
1 ” ao , 
= > cos n(o—6") | se—*"" J,,(sr’) J,,(sr)ds, 
2a nm Dp a |] (4) 
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where J,(x) is Bessel’s function of order n. For 
future application of this series it will be found 
desirable to express: the above integral as a 
contour integral in the complex z plane. This 
may be done by expressing one of the Bessel 
functions in terms of Hankel functions of the 
first kind. As is well known the Hankel functions 
HH, (x) and H,® (x) are also fundamental 
solutions of Bessel’s equation. The reiation® 
between J,(x) and //,,“ (x) is 


2J (x) =H1,? (x) —e"* HT,“ (—x). (5) 


With this relation the integral in Eq. (4) becomes 


ao 
e~*" J (sr’) J, (sr)sds 
e 


0 
= 


= | et? J (sr) HH, (sr)sds. 
e x 
Consider next the integral 
fe ate? 7 (sr HH, (er)adz 


taken over the contour indicated in Fig. 1. Since 
the integral has no poles in the positive half- 
plane the integral around the complete contour 
is zero, hence 


v . 


| et? J (sr’)H,, (sr) sds+ e~at? J (zr’) 


p Cc 
X IT, (sr)zdz+ integral over arcs=0. 


It may readily be shown that if r>r’ and the 
argument of z on the right-hand end of C is such 
that O=arg z<7/4 and at the left-hand (37/4) 
<arg z=7 then the integral over the circular 
arcs vanish as p—* so that understanding C to 
mean the path with the above characteristics 


et? J (sr) HH, (sr)sds 


« 


— -fe atz? J (sr) HT, (2r)edz, r>r’. (6) 
. 


For the case when r <r’ we need only interchange 
r and r’ in the integral on the right. 
7 


* Watson, Bessel’s Functions, p. 75. 
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The series in Eq. (4) thus reduces to 


t +. . 
—~ > cos n(d—86’) e~ ot? J,.(2r’) HH, (er)adz, 


Ar n x . Cc 


A Green’s function which vanishes for r’=R 
can be obtained by adding the solution 


1 


> cos n(@—86') A(z)e—** J, (2r') HH, (sr)adz 
4r n 0 “Cc (8) 


to solution (7), the A being determined so as to 
make the sum vanish at r’=R. It is readily seen 





Fic. 1. The path C in the z plane. 
that if A=+J,(2R)/H, (sR) then a solution 
which does vanish at 7’ =R is 
1, w+ . | 
3 cos n(@— 6") " oad 
Ag 2=—« a P 


HT, (gr) 


g(r,0;7r,0,)=- 


x iJ ,(2r) HH, (sR) 


IT, (sR) 
— HT, (2r’)J,(2R)}2dz, r>r’; (9) 


when r<r’ it is only necessary to interchange 
rand r’. 

In order that Eq. (9) be the correct Green's 
function. it must reduce to zero for all 7 and 7’ 
when ¢=0. That this condition is satisfied by the 
function (9) is easy to show by means of a 
suitable contour. The derivative of g is required 
in Eq. (2) and this is 


0g | 1 $a » L 
=— > cos n(@—@’) Sei 
ar’) vor An n x “Cc 
H,,°? (ar) d 
x H,,?(¢R) J,(2R) 
H,(2R) | d(zR) 
d 
—J,(zR) HT, (2R) } 27dz. 
d(zR) 


Making use of the relation 


d d 21 
J,(2)—H,, “9 (sz) — H,, ©? (z)—J,,(z) = 


dz dz WZ 
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the above equation simplifies to Substituting (9) and (10) in (2), the integration 
do a . H.. (er) with respect to 6’ can be carried out at once 
_= > cos n(@—6’) | e~ ate? zdz. since neither f(r’) nor u, depend on the angle 6’ 
Or 2n?R =o J ¢ H,,°? (sR) 


in our problem. This yields 
(10) : 


er 


. HH)“ (zr) 
u(r, t)= -3{ r' f(r')dr’ } e— ete" | Jo(sr’) Ho? (2R) — Ho (2r’) Jo(zR) | edz 
all “Cc TT, (¢gR) 


. . Hy? (2r’) 
—}3 r' f(r')dr’ re | Jo(sr) Ho (2R) — Ho"? (sr) Jo(2R) } adz 
er “Cc IT, : (sR) 


ai . _ Ho" (2r) 
a us(r)dr fe aie—t) zdz. (11) 
Te’9 ~“C Hy‘? (2R) 


Since we are interested in the heat flow the derivative of the above expression is required. Making 
use of the relations 


d d 
Jo(z)=—JSJi(z) and —H,“)(s) = —H, (sz) 
dz dz 
the derivative is 
Ou ” 7 IT, (sr) 
=} r'fir')dr" | gets? | Jo(sr’) Ho“? (2R) — Ho (2r’) Jo(2R) } 2*dz 
se Or op J ¢ Hy‘? (2R) 
? 
. . Ho‘ (zr’) 
+} r' f(r’ )dr’ e— ets? | Ji (sr) Hy (2R) — HW,“ (sr) Jo(2R) | 2*dz 
. J, Jo HT (sR) 


at 7 . 1 (ar) 
_ us(r)dr | ours z*dz. (12) 
Te’ 9 ~“C TT)? (sR) 


It is not profitable to proceed further with the general formulae without introducing the particular 
temperature variation at ‘=0, namely f(r), and the way in which the wall temperature changes with 
the time. 


APPLICATION TO A SPECIAL CASE 


The value of the derivative given by Eq. (12) will be computed for the special case where f(r) = uo 
=const. and where the function u,(¢), i.e. the temperature at the wall, is a step function of the time, 
that is, 


=U 5} O<t<hy 

U,=Uss th<t<te 

| —— he 
=Usn a, «4 


Substituting these special functions in Eq. (12) gives 


VOLUME 8, JUNE, 1937 443 





Ou Us . IT,“ (2r) 
; 2 ; o 
. rdr' | eat: ( To(sr’) Ho (¢R) — Ho (sr’)Jo(2R) | 22d2 


Or 2¢p Jo II)? (2R) 
Uo (* . HT)? (2r’) 
+ dr | e~st IJ (sr) Hy) (sR) — Hy? (sr) Jo(2R) | 2*dz 
2¢, J [To (eR) 
at afl . IT, l (or) af? . al . 
— we | dr |e lil 2°dz tu. | dr | -+-dzg+---4 ten | dr | ds} 
7 0 o < IT, ; (cR) v ti “Cc tat” é 





For r>R the integrals with respect to r’ and 7 can be evaluated thus: 


| or’ J (2r' dr’ =rJ\(2r) —RJ,(2R), 


R 
or’ Hy (ar) dr’ =rH(sr) — RH, (sR), 
. R 
. 
or’ Hy (2r' dr’ = —rH1, (ar) for values of z on path C 
e r 
alk 
a| s*e* dr=e—2t — ee a\t—tk-1 
Y th] 


Substituting these expressions in the above equation for du/dr one obtains: 


Ou Ue 11? (zr) : 
= Je ote {{rJ (sr) —RJ\(2R) JH (2R) — (rH, (er) — RH, (zR) JJ o(2R) | edz 
Or 2WJe I, (2R) 





Uy (* rll (2r) 
- e~ ate? IJ (er)? (2R) — Hy (2r) Jo(2R) | edz 
2 5 HT)? (2R) 
1 , IH, (2r) . - — Hy (ar) 
+ jaf fe ett —e~ ate") dz+ ue | § ea (tba) 2° galt ty) 2") dz 
rl J Ho? (sR) J ¢ HT? (2R) 
, _ (zr) 
+++++inf {1-e e(t—ta-1)2"} dz |. 
Jc HT)‘ (2R) 
> [1,"(2r) ; 
If it is noted that dz=0 for all r>R and then the limit of the result taken as r-R, the 
J ¢ I1o\? (sR) 
above expression becomes 
Ou uoR 7 H, (zr) 
lim = — Je ate? \Ji(2R)Ho?(2R) —1H1,? (2R)Jo(sR) | edz 
rok OF 2 we TI)? (cR) 
1 . 11"? (2R) , 1," (2R) 
+ —un fe aad dz+(wa—u.2) | e ili dz+-:- 
71 Jo Hy‘ (zR) ¢ Ho (2R) 
° HH," (2R) ° 11? (2R) 
(We Mare) | a(t—th) 2° dz+-*-+(u,, ten) fe o(t—tn—1)2° dz }. 
"J , IH) (2R) J Ho? (cR) 
Making use of the relation 
HH, (sR) Jo(sR) —J1(2R) Ho? (2R) = (—21/r2R) 
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ene yg 
- > 








and letting 


1 HH, (sR) 
F(t, R)= fe ate? dz 


11s C 
the above equation for (du/dr),.r becomes 


(0u/Or),—r = (Uo— Us) F(t, R) 


= (Ug—U51) F(t, R) + (tg: — U2) F(t—t, R) 


Hy‘? (2R) 


=(Uy—Ug) F(t, R) + (si — Use) F(t—y, RY He + +g n-1— Msn) F(t—tn_i, R) 


For R very large and uw,=u.2=-+-+ =U, the 
above derivative should approach that for the 
case of an infinite solid bounded by a plane. 
To check this, it is only necessary to take the 
limit as R approaches ~, i.e. 


lim (%o—u,,) F(t, R). 


R- 


For R very large the Hankel functions in 
F(t, R) may be replaced by their asymptotic 
expressions 


HH, (2R) = (2/2R) te! @8—-34/4) 


IT o I (<R) (2 rzR) ei Gk QT 4) 


II 


in which case 


1 
F(t, R)=- e~2'*dz, 


Te? CO 


Since the integrand has no poles the integral 
over C is the negative of that along the real axis 


or 
1 pet 1 
F(t, R)= - e atz* dy = 


T “ (rat)? 
Thus 


Ou Ug—tUs1 

lim —} = 

R-« 2 | 5 
OT) ma (mat) 


, 


which is exactly the value of du/dr for an in- 
finite solid bounded by a plane. 

In order to evaluate F(t, R) it is convenient 
to let 


F(t—t,, R) =(1/R)G(a—an), 
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0<t<t; 
t,<t<le 
(13) 
| aS oe 
where 
1 H(z) 
Gla) = fe aa*__ dz (14) 
Tid © IT)? (z) 


and a=at/R? and a,=at,/R?. 

The contour integral representing G(a@) can 
be transformed into an integral along the real 
axis as follows: 


1 H(z) 
G(a) = fe as" dz 
rie © TTy(? (sz) 


e7@ lz 


ads 


1e _(d/dz) H(z) 

—§ IT‘? (z) 
le d 

— e~*°— log TH“) (z)dz. 


Tle C dz 


Integrating by parts, 


2a 


ze~*" log Ho) (z)dz. 


Tle Cc 


G(a)= — 


Since the integrand has no poles in the positive 
half-plane and the integral over an indentation 
at the origin which is a branch point of HZ)“ (z) 
vanishes as the indentation shrinks to the origin 
the above integral becomes 


2a 7*@ , a 
G(a) =— } xe—** log Ho“ (x)dx 
11? _ 2 
2a 7 ' 9 ! . 
= -| xe~**" }log | Ho? (x) | +4P(x) } dx, 
Tid _.,, 


where y¥(x) is the argument of fH (x), ice. 
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tan (x) = Yo(x)/Jo(x) for x=0. Since | HW)“ (x)! 
is an even function the above integral reduces to 


) 


a¢ 


» ae 
G(a) = | xe 2" (x) dx. 
. 


7 . 


It may be shown that $(x)=yY(x)+7/2 is an 
odd function of x. Consequently 


xe? (x) —2/2!dx= xe~*"Y(x)dx 


‘ : e r 


= xe?" (x) dx 


and the expression for G(@) may finally be 
written 
. 


4a - 
G(a) = xe~*y(x)dx. (15) 
T +0 


APPROXIMATIONS FOR G(a) 


The integral G(a@) cannot be evaluated in 
terms of simple functions but a rather good 
approximation of it can be given. This can be 
done by noting that xt(x) is very nearly a 
linear function of x for values of x less than one. 
The approximation 


xiy(x) =Ax, O=x=1 (16) 


will be made for values of x <1 where A = 1.660. 
An idea of the agreement between xiy¥(x) and 
Ax can be had by consulting Table I. 

In the interval x>1, ¥(x) can be represented 
very closely by the first three terms of its asymp- 
totic expansion, namely 


¥(x) =x+7/4—1/8x. (17) 


A reference to Table I will indicate how well the 
asymptotic expression agrees with (x). Sub- 
stituting the approximate expressions in (15), 
G(a) becomes 


- fox ! - 
G(a)= A xle~*"dx 
T | ee’) 


7 (x?-+7/4x— fe “ds. (18) 


salt 


Carrying out the integration yields 
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TABLE I. 
7 1 
x ¥(x xdy (x) 1x +43, 
0 0 | O 
0.02 0.3719 0.0526 0.0332 
04 4416 | 0883 .0664 
.06 4923 | .1206 .0996 
08 5367 .1518 .1328 
10 .5765 .1823 .1660 
.20 7415 3316 .3320 
30 8806 .4823 4980 
40 1.0079 .6374 .6640 
50 1.1285 .7980 .8300 
.60 1.2446 .9641 .9960 
10 1.3577 1.1359 1.1620 
80 1.4686 1.3135 1.3280 
.90 1.5778 1.4968 1.4940 
1.00 1.6856 1.6856 1.6600 1.6604 
1.10 1.7924 1.7717 
1.20 1.8984 1.8812 
1.30 2.0037 1.9892 
1.40 2.1083 2.0960 
1.50 2.2125 2.2020 
1.60 2.3163 2.3072 
1.70 2.4197 2.4119 
1.80 2.5227 2.5159 
1.90 2.6256 | 2.6242 
2.00 2.7282 2.7229 


_ 8aA o (— 1)"a” 4\e-* 
G(a)= y -+- 1+ ) 
r n=0(4n+5)n! r/ 2 


1 a 
+ (:- ) (Ext (a)! (19) 
(7a)} a 


2 pla 
el Laie ods. 


(qr) % 


where 


0 


The series is convergent for all values of a but 
is not particularly convenient for computations 
for large values of a so that another expression 
will be given for large a. It may be noted that 
for very small a the above expression becomes 
1/(ra)+ which it should do in order to give the 
flow of heat from an infinite solid bounded by a 
plane when R-~, i.e. a—0. The above expres- 
sion has been checked with G(a) for several 
values of a by integrating Eq. (15) numerically. 
The agreement is quite good for values of @ less 
than 50. For values of a> 50 it is probably best to 
integrate expression (15) numerically. 

For large values of a, i.e. a greater than 3 or 4, 
it may be shown that 
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vl 1 
| 3 2e az“dy ms 
“0 


| (5/4) 
i 
1 3 
—a'‘e “(14 = +) 
4a 16a" 


: e 1 3 
(1—Erf (@)*)= Wen iam —) 
(ra)? 2a (2a)? 


Substituting these expressions in Eq. (18) one 
finds that 


AT(5/4) | r 7 
—e* (4- ~ ) 
| a | 4 8 
1/4A —9 1 /11—6A | 
+-( }- ( )r- | (20) 
a 16 a’ 32 


where I'(5/4) is the gamma-function of 5/4 or 
(5/4) =0.9064. 


2a° 


and 


2 


Gia)= 





7 


TotraL HEAT FLow 


The total amount of heat which has flowed 
through the cylinder wall during the time inter- 
val (=0 to {=¢ is given by 


» OU | 
O=k dt, 
vn ar|, R 
where k is the heat conductivity. From Eq. (13) 
this becomes 


at 
Q=k(uo—ua) | Flt, Rode O<t<t, 


0 


t 


F(t, R)dt 


-A| (uo— ter) | 
“0 


t 
(amma) f Fut, Rid Iy<t<lo 


e/ 1} 


=H (wy) f F(t, R)dt 


0 


+(ua—wes) | F(t—t, R)dt+--- 
e/ tt} 


vf 
+ (sn tn) | F(t—t,-1, R)dt 
¢Y tn—1 


bn—1<b<bn. 
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Consequently an integral of the type 


*t 
F(t—tn, R)dt 


tm 
must be evaluated. From Eq. (14) it is seen that 


1 vf 
| G(a—a, )dt 


| ‘F(t—tn, R)dt= 
R © tm 


tm 


R poem R 
= | G(a)= K(a—a,). (22) 


ade 0 a 


With the approximate expression (19) the inte- 
grations indicated above can be carried out 
giving 

(—1)"(a—a,,)* 


? 2 
K(a—a,,) = 4.1 (a—aQ»)? > 
n—0(4n+5)(n+2)n! 


nr 23 (a— Aa!) 
T + (1 —e— e-em) ) + gq enee 


a— mn 
+(r(a—a,,))?f 1—- (1—Erf (a—a,,)*) 
12 


(23) 


Although the series in the above expression con- 
verges for all values of (a—a,,) it is not con- 
venient to use for large values of (a—am) so 
that an asymptotic series is desirable for large 
values of (a—a,,). This may be obtained most 
conveniently by going back to the first integral 
in Eq. (18) from which the above series origi- 
nated. It is then seen that a value of the integral 


2 a— am l 
[= | ada f xie—2dx 
bel | “0 


is required for large values of (a—a,,). By means 
of an obvious transformation J may be written 


8 — a 


[= f ida | sie ‘ds, 


0 0 


(24) 


where B=a—a,,. The integral 


2a 


ste—*ds 


0 


is just the integral representation of the in- 
complete gamma-function which with the help 
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of the confluent hypergeometric function’ may 
be written as follows: 


sie*ds = 1'(5/4) —a’'e~* 


(1+ £/a)'e‘*dé. 
0 
Substituting this in Eq. (24) and inverting the 
order of integration and integrating with respect 


to a one finally obtains 


4 4 
[=-1(5/4)B3!4 


é | | 
3 3 p e 


Expanding the integrands in powers of y/8 and 


integrating term by term, the required asymp- 


(- 7\1 

2— + 

ahi 

‘ 11\ 1 a 15\ 1 | 
4\4 ae 44\4 “& 


Using this expression the complete expression 


for Kia Am) is 


totic series is 


4 4 €° 
I = 1*( 5 4 )B° on + 
3 





w 
~~ 


8A -. ea 





‘Whittaker and Watson, p. 341. 
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a Om 
i i ~ - 1 
+ (7(\ a — Am) )? 1- (1—Erf (a—a,,)?) 
12 


B=(a—am) 3. 


For large values of a, i.e. 4<a< 50, the above 
expression is very convenient and quite accurate 
since many of the terms yield only small con- 
tributions. 

the dimensionless 


am) the expression (21) for the total heat 


In terms of 
K(a- 


which has flowed up to the time ¢ per unit area 


quantity 


of tunnel wall is 


RR 
Q= Ug —Us,)K(a—O) 0<t<t, 
a 
kR 
= [ (uo — U1) K(a—0) + (51 — U2) K(a—a) | 
a 
t) < t < lo 
kR 
= [ (to — U1) K(a—O0) + (5) — U2) K(a@— a) 
a 
7° °° -(08,, 1—Us,)K(a-—a, 1) | Fans <5 <8, 


It does not seem practical to carry the analysis 
any further for the special case where the tunnel 
wall temperature is a step function of the time. 
The above formula is quite useful in all those 
which the wall 


cases for temperature can be 


approximated by a step function. 
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